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ABSTRACT 


Josephson  Parametric  Amplification  for 
Circuit  Quantum  Electrodynamics: 

Theory  and  Implementation 

Christopher  Michael  Quintana 
Department  of  Physics 
Bachelor  of  Arts 

Circuit  Quantum  Electrodynamics  (cQED)  is  a  highly  promising  platform  for  quantum 
control,  quantum  simulation,  and  quantum  optics  in  a  solid-state  environment.  Such  sys¬ 
tems  consist  of  artificial  atoms  coupled  to  superconducting  resonators,  and  can  be  probed 
with  microwave  signals  on  a  chip.  This  thesis  theoretically  and  experimentally  studies  a 
method  of  quantum-limited  Josephson  parametric  amplification  for  the  improvement  of  such 
microwave  measurements.  A  theoretical  analysis  of  parametric  amplification  with  a  nonlin¬ 
ear  quantum  oscillator  is  given,  along  with  an  analysis  of  a  particular  implementation  with 
a  superconducting  coplanar  waveguide  resonator  shunted  by  an  array  of  SQUIDs.  A  para¬ 
metric  amplifier  chip  is  designed  and  a  fabrication  process  is  engineered.  Measurements  and 
analysis  of  a  successfully  fabricated  device  are  presented,  and  future  outlook  is  discussed. 
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Chapter  1 


Introduction  and  Motivation 


1.1  Quantum  Computation 

Building  a  quantum  information  processor  is  currently  an  active  area  of  experimental  re¬ 
search  in  physics  and  engineering.  In  the  standard  model  of  quantum  computation,  the 
familiar  logical  bit  (0  or  1)  of  a  classical  computer  is  replaced  by  a  two- level  quantum  sys¬ 
tem  to  form  a  quantum  bit,  or  qubit ,  with  computational  basis  states  denoted  by  |0)  and 
|1).  Algorithms  that  take  advantage  of  uniquely  quantum  behavior,  such  as  superposition 
of  basis  states  and  entanglement  of  qubits,  allow  certain  computational  tasks  to  be  com¬ 
pleted  in  an  inherently  more  efficient  way  than  any  known  (and  in  some  cases  any  possible) 
classical  algorithm.  Two  now-famous  examples  are  Shor’s  algorithm,  which  uses  a  quantum 
mechanical  speedup  in  taking  a  discrete  Fourier  transform  to  factor  numbers  exponentially 
faster  than  we  know  how  to  classically,  and  Grover’s  algorithm,  which  yields  a  quadratic 
speedup  in  the  brute-force  solution  of  NP  problems  (essentially  any  problem  whose  answer 
can  be  efficiently  checked)  [1,2].  The  techniques  of  quantum  computation  are  also  promising 
as  a  means  of  studying  the  behavior  of  computationally  intractable  quantum  systems. 

Since  any  two-level  quantum  system  will  in  principle  function  as  a  qubit,  many  physical 
architectures  have  been  and  are  still  being  pursued.  There  are  certain  obvious  requirements 
that  any  physical  implementation  should  satisfy,  widely  known  as  the  DiVincenzo  criteria  [3]: 

(Dl)  Scalable,  well-characterized  qubits. 

(D2)  Initialization  of  qubits;  say,  to  1 000. . .) . 

(D3)  Long  coherence  times  relative  to  gate  speeds. 

(D4)  A  universal  set  of  quantum  gates. 

(D5)  Qubit  state  readout. 
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Criterion  (Dl)  means  that  we  should  have  qubits  in  the  first  place.  If  we  don’t  have  a  nat¬ 
ural  two- level  system  like  a  spin-);  particle,  we  require  a  quantum  system  with  anharmonic, 
isolatable  energy  levels  so  that  a  two-level  subspace  approximation  can  be  made.  This  way, 
our  computations  will  only  drive  transitions  between  two  chosen  levels  (although,  the  clever 
use  of  higher  auxiliary  levels  can  aid  in  computations).  Requirement  (Dl)  also  means  that 
we  should  know  the  qubit  parameters  necessary  for  the  performance  of  the  quantum  gates 
of  (D4).  In  the  prominent  gate  model  of  quantum  computation,  “gates”  are  simply  pre¬ 
determined  unitary  time-evolution  operations  on  quantum  states  resulting  from  carefully 
engineered  Hamiltonians.  “Universal”  in  (D4)  refers  to  the  fact  that  any  quantum  gate  can 
be  written  as  a  composition  of  single-qubit  gates  and  a  single  type  of  entangling  two-qubit 
gate  [3]  (i.e.,  a  gate  which  cannot  be  written  as  the  tensor  product  of  two  one-qubit  gates). 
This  requirement  poses  the  necessary  challenge  of  coupling  different  physical  qubits  together 
for  the  generation  of  entanglement.  The  controllable  coupling  of  more  than  a  few  qubits 
has  shown  to  be  an  extremely  difficult  experimental  task.  Decoherence  of  fragile  quantum 
states  due  to  internal  noise  and  unwanted  interaction  with  the  environment  puts  a  limit  on 
how  many  quantum  gates  can  be  performed  before  the  quantum  states  change  undesirably 
by  decaying  to  their  ground  states  (relaxation)  or  losing  their  relative  complex  phase  infor¬ 
mation  (dephasing).  In  fact,  it  is  surprising  that  useful  quantum  computation  is  thought 
to  be  possible  at  all  in  the  face  of  decoherence  -  remarkably,  feedback-based  quantum  error 
correction  schemes  have  been  developed  which  show  that  below  certain  error  and  fidelity 
thresholds,  computations  can  succeed  with  arbitrarily  high  probability  [4],  In  particular, 
the  use  of  so-called  “surface  codes,”  in  which  fault  tolerance  emerges  out  of  a  fabric  of  logical 
qubits  formed  by  many  physical  qubits  with  nearest-neighbor  coupling,  is  a  highly  promising 
(though  highly  ambitious)  route  to  take  with  superconducting  qubits  [5,6].  The  ability  to 
quickly,  reliably,  and  nondestructively  read  out  quantum  states  is  therefore  an  essential  part 
of  (D5).  Even  one  working  quantum  computer  would  be  an  incredibly  marvelous  accom¬ 
plishment.  Due  to  these  difficulties  and  the  specific  nature  of  known  quantum  algorithms 
(not  to  mention  the  expensive  laboratory  equipment  needed),  one  shouldn’t  expect  to  see 
personal  quantum  computers  in  the  future,  but  rather  servers  allotting  (valuable)  quantum 
computation  time  to  users,  not  unlike  in  todays  supercomputing  clusters. 

1.2  Circuit  QED:  Quantum  Optics  on  a  Chip 

This  thesis  focuses  on  a  method,  to  be  described  shortly,  of  improving  (D5)  in  arguably  the 
most  promising  implementation  of  quantum  computation:  circuit  quantum  electrodynamics 
(cQED),  a  field  which  has  seen  incredible  progress  over  the  past  decade.  As  illustrated 
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(a)  Cavity  QED  (CQED) 


(b)  Circuit  QED  (cQED) 


Figure  1.1  (a)  Cavity  Quantum  Electrodynamics  (CQED):  A  two-level  atom  interacting 
with  a  single  mode  of  the  electromagnetic  field  with  coupling  strength  g.  n  is  the  rate 
at  which  photons  decay  out  of  the  cavity,  set  by  the  transparency  of  mirror-like  walls. 
7_l  is  the  rate  at  which  the  atom  decays  into  modes  not  trapped  by  the  cavity  via  other 
mechanisms.  Xtransit  is  the  amount  of  time  that  the  atom  spends  in  the  cavity,  (b)  Circuit 
QED  (cQED):  Modern  solid-state  offspring  of  CQED  in  which  a  spatially  fixed  artificial 
atom  interacts  with  a  resonator  formed  by  a  coplanar  waveguide  circuit  on  a  chip.  Figure 
adapted  from  [7]. 


in  Figure  1.1,  cQED  stems  from  cavity  quantum  electrodynamics  (CQED),  the  study  of 
atoms  interacting  with  photons  in  cavities,  which  was  well-studied  before  it  was  considered 
for  quantum  computation  [8].  In  CQED,  the  interaction  of  an  atom  with  a  single  mode  of 
the  electromagnetic  field  can  cause  a  quantum  of  energy  to  coherently  oscillate  between  an 
atomic  excitation  and  a  cavity  photon,  a  phenomenon  known  as  vacuum  Rabi  oscillation. 
This  process,  depicted  in  Figure  1.1(a),  occurs  at  a  rate  related  to  the  coupling  strength  g 
until  photons  leak  out  of  the  cavity.  This  phenomenon  has  no  obvious  classical  counterpart 
and  requires  a  quantum  theory  of  light  to  explain.  To  a  good  approximation,  such  a  system 
is  described  by  the  Jaynes- Cummings  Hamiltonian 1 


#JC 


hujr 


a) a  +  -  j  +  -fujOacrz  +  hg((T+a  +  cr_ 


a 


(1.1) 


where  a  and  <r_  are  lowering  operators  for  the  cavity  mode  and  the  atom,  respectively. 
In  circuit  QED,  the  atom  and  cavity  are  replaced  by  a  solid-state  “artificial  atom”  and 
(usually)  a  coplanar  waveguide  resonator  on  a  chip.  This  modern  development  has  among 
other  advantages  the  capability  of  being  conveniently  fabricated  on  a  chip  with  standard 
clean  room  techniques,  allowing  for  a  significant  amount  of  man-made  control  over  design 
parameters  and  for  artificial  atoms  that  are  spatially  stationary. 

It  can  be  shown  (see  Appendix  B.2)  that  when  the  detuning  |A|  =  \uja  —  ujr  \  of  the  atom 

^Tor  a  discussion  of  electromagnetic  field  quantization  and  an  in-depth  derivation  of  (1-1),  the 
interested  reader  may  consult  Appendices  A  and  B. 
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from  the  resonator  mode  is  large  compared  to  the  coupling  g ,  under  a  change  of  basis  the 
“dispersive”  Hamiltonian  can  be  well-approximated  by 
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(1.2) 


Comparing  the  first  term  with  the  form  of  a  standard  harmonic  oscillator  Hamiltonian,  we 
can  interpret  the  effect  of  the  dispersive  coupling  as  an  effective  qubit  state-dependent  shift 
in  the  cavity’s  resonance  frequency.  By  probing  the  cavity  near  resonance  with  a  signal 
weak  enough  not  to  disturb  the  qubit,  one  can  reliably  read  out  the  qubit’s  state  [8].  In 
fact,  the  above  Hamiltonian  generalizes  to  the  case  of  several  qubits  coupled  to  a  single 
cavity  mode,  and  the  dispersive  readout  method  can  be  used  to  infer  the  joint  state  of 
multiple  qubits  at  once  [9,10].  The  dispersive  regime  also  provides  an  effective  coupling 
between  two  qubits  coupled  to  the  same  resonator,  allowing  the  performance  of  two-qubit 
gates  [8,9,11].  “Building  a  quantum  computer”  is  a  phrase  that  has  been  thrown  around  a 
lot  over  the  past  decade  or  so,  but  in  just  the  past  couple  of  years,  very  impressive  progress 
in  superconducting  qubit  coherence  times  and  scalability  (for  example,  [12-14])  have  made 
this  goal  no  longer  insane  (for  lack  of  a  better  word),  but  rather  only  insanely  (for  lack  of 
better  word)  difficult.  If  steady  progress  continues,  it  is  not  unreasonable  to  speculate  that 
a  large-scale  practical  quantum  computer  could  be  built  within  ~  20  years  from  now. 

Besides  its  potential  use  for  quantum  information  processing,  circuit  QED  offers  an 
exquisite  solid-state  test  bed  for  the  rich  field  of  quantum  optics  in  the  microwave  regime. 
A  special  case  of  this  is  the  dispersive  readout  of  qubits  mentioned  above.  Because  cQED 
systems  have  characteristic  energies  in  the  microwave  regime  (much  lower  than  optical  ener¬ 
gies),  it  is  a  difficult  task  to  detect  signals  at  the  level  of  a  few  photons.  For  nondestructive 
quantum  feedback  and  other  interesting  quantum  optics  experiments,  it  is  necessary  to 
probe  systems  using  signals  of  low  average  photon  number.  Commercial  amplifiers  can 
work  for  many  of  these  experiments,  at  the  expense  of  adding  significant  noise  and  thereby 
requiring  long  averaging  times  and  potentially  disturbing  the  quantum  system  being  probed. 
This  in  general  arises  from  the  lack  of  control  over  many  microscopic  degrees  of  freedom 
(such  as  electrons  in  a  high-electron  mobility  transistor  [HEMT])  that  need  to  be  coupled 
to  the  input  signal  for  energy  transfer  to  occur.  It  would  be  nice  to  have  a  linear  amplifier 
with  a  noise  temperature  below  1  K  at  microwave  frequencies.  In  accordance  with  Frii’s 
formula  for  a  chain  of  amplifiers  [15, 16],  an  initial  amplification  with  a  low  noise  figure  and 
high  gain  will  lead  to  an  effective  SNR  at  later  stages  of  amplification  determined  by  the 
added  noise  figure  of  the  first  stage  (and  the  SNR  of  the  original  input  signal). 

To  date,  the  most  successful  microwave  amplifiers  have  been  based  on  so-called  para- 
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metric  amplification  [17],  which  can  achieve  the  best  added-noise  performance  allowed  by 
quantum  mechanics.  In  general,  parametric  amplification  is  a  phenomenon  which  occurs 
when  some  reactive  parameter  oscillates  near  twice  the  frequency  of  an  oscillating  system, 
allowing  energy  to  be  pumped  into  or  out  of  an  input/output  port  to  create  an  amplifier. 
Such  an  oscillation  can  naturally  be  achieved  in  systems  consisting  of  a  nonlinear  oscilla¬ 
tor.  In  the  context  of  superconducting  electrical  circuits,  such  microwave  detectors  were 
first  explored  in  the  late  1980s  by  Yurke  et.  al.  [18,19].  In  recent  years,  interest  in  such 
quantum-limited  amplifiers  based  on  Josephson  junctions  has  seen  a  great  resurgence  due  to 
applications  to  quantum  information  and  astronomical  metrology  [20-26] ,  with  each  variant 
of  Josephson  parametric  amplifier  having  its  own  strengths  and  weaknesses.  As  a  result,  in 
very  recent  years,  impressive  experiments  have  been  enabled  by  such  amplifiers,  including 
the  realtime  monitoring  of  quantum  jumps  in  a  superconducting  qubit  [27],  feedback  con¬ 
trol  and  state  preparation  of  superconducting  qubits  [28,29],  the  measurement  of  nontrivial 
microwave  photon  correlation  statistics  [30,31],  significant  advances  in  quantum  measure¬ 
ments  [25],  and  the  monitoring  of  a  qubit  state  via  a  partial  back-action  measurement  [32]. 
Quantum-limited  microwave  amplifiers,  then,  allow  interesting  quantum  optics  experiments 
to  be  performed  that  would  otherwise  be  too  difficult  due  to  a  loss  of  detection  efficiency.  In 
general,  better  amplifiers  lead  to  more  efficient  experiments.  In  particular,  the  implemen¬ 
tation  of  quantum-limited  amplifiers  for  quantum  feedback  is  an  extremely  important  and 
active  area  of  research  in  the  current  effort  to  achieve  scalable  quantum  error  correction. 

1.3  Thesis  Overview 

The  goal  of  this  thesis  is  to  thoroughly  study  the  theory  behind  Josephson  parametric  am¬ 
plification,  engineer  a  physical  design,  and  experimentally  make  and  test  such  an  amplifier 
in  the  laboratory.  In  Chapter  2,  an  in-depth  treatment  of  the  theory  of  parametric  am¬ 
plification  with  a  quantum  Duffing  oscillator  will  be  given.  Unlike  in  this  introduction,  an 
attempt  to  be  thorough  and  guide  the  reader  through  all  of  the  important  physics  will  be 
made.  We  will  explore  the  origins  of  parametric  amplification  in  nonlinear  quantum  sys¬ 
tems  and  analyze  important  practical  considerations  and  tradeoffs  such  as  gain,  bandwidth, 
dynamic  range,  phase-sensitivity,  and  added  noise.  In  particular,  it  will  be  shown  that  this 
method  of  parametric  amplification  can  in  principle  yield  an  amplifier  that  adds  the  least 
amount  of  noise  allowed  by  quantum  theory.  In  the  following  chapter  (Chapter  3),  an  in- 
depth  theoretical  treatment  of  a  particular  implementation  with  superconducting  circuits 
is  discussed;  namely,  how  a  quantum  Duffing  oscillator  naturally  arises  out  of  a  system 
composed  of  an  array  of  SQUIDs  terminating  the  centerpin  of  a  A/4  coplanar  waveguide 
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resonator.  Such  an  amplifier  can  be  driven  with  a  single-mode  pump  field  of  reasonable 
power  and  also  have  its  amplified  band  tuned  in  frequency  via  the  SQUIDs.  Finally,  in 
Chapter  4,  we  will  first  describe  in  detail  the  development  of  a  process  to  physically  make 
such  an  amplifier  from  scratch  using  clean  room  techniques.  Various  practical  design  fac¬ 
tors  are  considered  based  on  the  theory  developed  in  the  previous  two  chapters.  Finally, 
measurements  characterizing  a  successfully  fabricated  Josephson  parametric  amplifier  are 
presented.  The  behavior  of  the  parametric  amplifier  is  analyzed  and  future  prospects  for 
using  it  as  a  tool  in  other  experiments  are  discussed. 


Chapter  2 


Parametric  Amplification 


2.1  A  Classical  Intuition:  Parametric  Resonance 

Before  delving  into  a  thorough  quantum  mechanical  treatment,  in  order  to  obtain  an  in¬ 
tuition  behind  parametric  amplification  we  will  briefly  discuss  the  classical  phenomenon  of 
parametric  resonance,  to  which  we  will  eventually  compare  the  quantum  case.  A  paramet¬ 
ric  oscillator  is  one  whose  parameters ,  such  as  a  spring  constant  or  damping  term,  vary  in 
time  [33,34],  A  simple  example  of  such  a  system  is  a  small-angle  pendulum  with  a  vertically 
oscillating  support  point,  or  a  child  (partially)  pumping  a  swing  already  in  motion  by  peri¬ 
odically  standing  and  squatting.  Here,  the  time-varying  length  of  the  pendulum  multiplies 
the  dynamical  variable  0  in  its  equation  of  motion.  In  the  context  of  the  classical  spring 
oscillator  studied  in  freshman  physics,  a  parametric  oscillator  has  the  second-order  linear 
differential  equation  of  motion 

x  -f-  2y(t)x  -(-  ———x  =  0.  (2.1) 

m 

Without  loss  of  generality  we  can  assume  7  (f)  =  0  by  making  the  change  of  variables 
q(t)  =  e^o  7('r)drx(f)  [33,35],  effectively  moving  into  an  “undamped”  frame  in  which  the 
equation  of  motion  becomes 

q  +  vi(t)q  =  0,  (2.2) 

where  =  k(t)/m  —  y2(f)  —  "y(t)  is  the  effective  natural  frequency  of  the  (/-oscillator. 

When  ui(t)  is  a  periodic  function  of  time,  this  is  known  in  mathematics  as  a  Hill  equation. 
Let  us  suppose,  as  is  often  the  case,  that  uq(f)  is  varied  via  a  sinusoidal  perturbation.  We 
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write  this  in  the  canonical  form 

wf(t)  =  ljq  [1  +  h cos(2tOpf)] ,  (2.3) 

where  h  <C  1.  With  this  form  for  t c2(t),  equation  (2.2)  is  known  in  mathematics  as  the 
Mathieu  equation.  To  gain  an  intuitive  understanding  of  parametric  resonance,  we  substitute 
(2.3)  into  (2.2)  and  bring  the  oscillating  parameter  to  the  right  hand  side, 

q  +  oo^q  =  —Hujq  cos(2 uipt)q.  (2-4) 

This  represents  a  driven  simple  harmonic  oscillator,  with  an  effective  driving  term  propor¬ 
tional  to  the  response  q.  Starting  for  the  moment  with  h  =  0,  the  solution  is  simple  harmonic 
motion,  q(t)  =  acos(wot)  (choosing  the  appropriate  origin  of  time).  As  an  approximation, 
we  can  substitute  q  =  acos(cuo)t  on  the  right  hand  side, 

q  +  ui^q  =  —huj^a  cos(2cupt)  cos(u;ot) 

=  —Hujq-  [cos([2tUp  —  uJo]t)  +  cos([2cjp  +  wo]£)]  .  (2-5) 

Resonance  will  occur  when  the  effective  driving  frequency  is  equal  to  the  effective  natural 
frequency  u>o.  We  see  that  this  happens  in  the  first  term1  when  2up  =  2loq  and  the  sec¬ 
ond  term  then  becomes  an  effective  off-resonant  drive  at  frequency  3<uo  which,  as  studied 
in  freshman  physics  [36],  is  heavily  filtered  out  and  can  be  ignored  to  a  good  approxima¬ 
tion.  Therefore,  if  (t)  is  perturbed  with  a  frequency  close  to  2wo  (i.e. ,  twice  the  natural 
frequency  of  the  oscillator),  we  effectively  have  a  simple  harmonic  oscillator  driven  near 
resonance,  with  the  driving  force  proportional  to  the  amplitude  of  existing  oscillations.  But 
this  suggests  that  the  amplitude  of  oscillation  will  grow  exponentially  (unless  the  ampli¬ 
tude  is  initially  zero)!  Roughly,  this  idea  is  known  as  parametric  resonance.  Of  course,  our 
analysis  was  by  no  means  rigorous  and  is  not  exact,  but  this  is  in  fact  what  happens  -  it 
can  be  shown  through  a  more  thorough  analysis  [33]  that  for  a  parameter  modulation  of 
frequency  of  2cuo  +  e  (where  |e|  <  o  is  small),  the  amplitude  of  the  non-transient  sinu¬ 
soidal  response2  at  frequency  ojq  +  \e  is  proportional  to  est,  where  s  =  \  \J (^/tcuo)2  —  e2. 
When  we  add  damping  [as  in  (2.1)],  the  proportionality  constant  with  respect  to  the  lab 

1  Alternatively,  the  first  term  would  satisfy  a  degenerate  “resonance”  condition  if  and  only  if 
Wp  =  0,  which  means  no  parametric  variation  in  the  first  place. 

2  Consistent  with  Floquet  theory  [37],  the  response  will  also  have  Fourier  components  at  frequen¬ 
cies  (2 n  +  l)(wp  +  \e)  for  nonzero  n,  but  it  can  be  shown  that  the  amplitudes  of  these  components 
are  very  small  in  the  steady  state.  For  example,  the  leading  third-harmonic  correction  at  frequency 
3(wo  +  |e)  is  smaller  in  amplitude  by  a  factor  of  h/ 16  from  the  fundamental  harmonic  response  [33]. 
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frame  dynamical  variable  x  becomes  e^s~'Y',t.  We  therefore  see  that  to  a  good  approximation 
the  parametric  oscillation  manifests  itself  as  an  effective  negative  damping  constant  that 
competes  with  (and  if  large  enough  can  overcome)  the  natural  damping  of  the  system. 

To  roughly  see  how  parametric  resonance  can  be  utilized  for  amplification,  consider 
adding  an  external  driving  term  to  the  system 


x  +  2-fx  +  u2(t)x 


Fq 

m 


cos  (Udt), 


(2.6) 


and  then  (naively)  approximate  the  effect  of  the  parametric  oscillation  as  a  modification  of 
the  damping  term  of  a  non-modulated  resonator, 


X  +  2(7  —  s)x  +  CJqX 


F°  (  A 
-  COS \UJdt). 

m 


(2.7) 


If  the  parametric  oscillation  strength  h  is  such  that  the  effective  damping  constant  7  —  s 
is  still  positive  but  small,  then  from  the  study  of  the  forced  damped  harmonic  oscillator  of 
freshman  physics  [36] ,  the  resulting  steady-state  amplitude  of  oscillation  at  frequency  uq  is 


a  = 


Fq _ 1 _ 

m  \J i^d  ~  ^o)2  +  (2 [7  -  s\ujd)2 


(2.8) 


showing  a  Lorentzian  resonance  near  u)d  =  wo  with  a  maximum  amplitude  proportional  to 
1/(7  —  s) .  For  strong  enough  parametric  oscillation,  s  can  approach  7  and  yield  a  large 
parametrically  amplified  steady-state  response,  and  if  s  is  large  enough  to  overcome  7,  the 
system  response  will  continue  to  increase  in  amplitude  (even  without  an  explicit  driving 
term,  as  discussed  above).  This  is  by  no  means  an  exact  solution,  and  is  only  meant  to  give 
an  intuition  behind  the  origins  of  parametric  gain  (a  rigorous  treatment  of  the  quantum 
case  will  be  given  later). 

Such  an  amplification  scheme  has  several  advantages.  It  can  easily  be  made  very  sensi¬ 
tive  while  simultaneously  having  a  low  level  of  noise,  assuming  one  can  pump  the  parameter 
u(t)  in  a  low-noise  manner,  and  that  a  reactive  rather  than  resistive  parameter  is  made 
to  oscillate.  This  can  for  example  be  done  inductively  or  with  an  electromagnetic  drive. 
This  pump  tone  (i.e.,  the  signal  that  modulates  the  system  parameters)  is  ultimately  what 
provides  the  energy  for  the  amplification  of  signals.  Besides  forming  a  natural  cryogenically- 
operable  parametric  amplifier,  we  will  see  through  a  quantum  mechanical  treatment  that 
the  parametric  amplifier  implemented  in  this  thesis  in  fact  theoretically  adds  the  minimum 
amount  of  noise  allowed  by  quantum  mechanics,  allowing  for  example  signals  at  the  one- 
photon  level  and  for  signal-to-noise  ratios  to  overcome  the  limitations  of  commercial  HEMT 


2.2  Parametric  Amplification  with  a  Quantum  Duffing  Oscillator 


10 


amplifiers.  In  addition,  when  operated  in  the  quantum  regime,  such  an  amplifier  can  let  us 
access  some  very  interesting  quantum  phenomena,  such  as  squeezing. 

2.2  Parametric  Amplification  with  a  Quantum  Duff¬ 
ing  Oscillator 

How  can  we  make  a  parametric  amplifier  (paramp3  for  short)  for  microwave  signals  in 
cQED?  A  practical  way  to  achieve  parametric  amplification  is  by  employing  a  nonlinear  os¬ 
cillator,  which  can  have  a  bifurcation  point  near  which  the  resonator  response  depends  very 
strongly  on  the  input  [21].  By  perturbing  a  strong  pump  signal  sent  into  such  a  nonlinear 
cavity,  one  can  achieve  parametric  amplification.  If  we  pursue  this  practical  method,  the 
question  then  becomes  how  to  create  such  a  nonlinear  quantum  oscillator.  As  we  will  see  in 
detail  in  Chapter  3,  the  effect  of  embedding  Josephson  junctions  in  a  coplanar  waveguide 
resonator  is  to  give  rise  to  a  nonlinear  resonator  that  can  be  treated  quantum-mechanically 
under  experimental  conditions.  Before  delving  into  this  specific  implementation,  we  study 
the  general  theory  of  parametric  amplification  with  a  quantum  Duffing  oscillator,  which 
will  allow  us  to  use  the  powerful  technique  of  input-output  theory  from  quantum  optics  to 
describe  how  we  can  send  input  microwave  signals  to  such  a  resonator  and  receive  amplified 
output  signals.  Input-output  theory  also  provides  a  straightforward  way  to  include  both  lin¬ 
ear  and  nonlinear  loss,  as  well  as  other  general  nonlinearities4.  Finally,  as  already  mentioned, 
a  fully  quantum  treatment  will  allow  us  to  see  that  such  amplifiers  are  quantum-limited, 
in  principle  yielding  the  least  amount  of  noise  allowable  by  the  quantum-mechanical  un¬ 
certainty  principle.  For  these  reasons,  we  will  study  parametric  amplification  in  its  general 
quantum  context,  allowing  us  to  derive  among  other  things  expressions  for  the  bandwidth 
and  gain  of  a  parametric  amplifier  derived  from  a  general  class  of  nonlinear  effects. 

The  Hamiltonian  of  the  quantum  Duffing  oscillator  follows  naturally  from  its  classical 
counterpart.  We  start  with  the  classical  Duffing  oscillator  Hamiltonian 

P(q,p)  =  +  ^rnulq2  +  ^ffq4,  (2.9) 

and  convert  the  first  two  terms  into  a  quantum  Hamiltonian  using  the  standard  method 
of  annihilation  and  creation  operators  (which  we  will  denote  by  A  and  AQ  studied  in 

3We  will  often  use  this  abbreviation  for  brevity,  even  though  “para-amp”  may  be  more  logical. 

4In  particular,  it  allows  us  to  include  the  nonlinear  kinetic  inductance  of  the  superconducting 
transmission  line  itself,  although  we  will  see  that  such  nonlinearities  are  negligible  in  our  parametric 
amplifier. 
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introductory  quantum  mechanics  in  the  context  of  the  simple  harmonic  oscillator.  For  the 
third  term,  we  use  the  fact  that  q  oc  {A  +  AQ  to  obtain  the  quantum  Duffing  Hamiltonian 
operator 

-fflCerr  =  frhJoA^  A  +  K—{A  +  H^)4 

PS  hcc0HtH  +  hy(Ht)2H2,  (2.10) 

where  in  the  second  equality  we  have  dropped  rapidly  oscillating  terms  whose  net  time- 
averaged  effect  is  zero  in  the  rotating  wave  approximation.  Seen  in  another  way,  we  can 
drop  these  terms  provided  K{A^  A)  isn’t  large  compared  with  coq,  because  they  represent 
non-energy-conserving  processes  with  an  unbalanced  number  of  photon  creations  and  an¬ 
nihilations.  With  the  normalization  chosen  as  above,  K  is  usually  referred  to  as  the  “Kerr 
constant,”  by  analogy  to  the  classical  nonlinear  Kerr  effect  in  which  a  material  has  an 
electric  field-dependent  index  of  refraction  n  =  no  +  ri2E2  [38]. 

2.2.1  Input-Output  Theory 

Suppose  we  have  a  quantum  Kerr  oscillator.  To  utilize  such  an  oscillator  (which  we  will 
often  refer  to  as  the  cavity)  as  an  amplifier,  we  need  to  couple  it  to  the  outside  world  -  in 
particular  via  an  open  port  coupling  the  oscillator  to  a  microwave  transmission  line.  But 
before  specializing  to  the  case  of  a  coupling  capacitor,  we  consider  the  general  case  of  a  single 
cavity  mode  interacting  with  an  external  bosonic  field,  and  derive  the  relation  between  input 
and  output  signals  to  and  from  the  cavity  mediated  by  this  field.  To  do  this,  we  will  use 
the  method  of  input-output  theory  developed  in  chapter  7  of  [39].  Namely,  we  will  model 
the  transmission  line  used  to  drive  the  cavity  by  a  collection  of  harmonic  oscillators  6(w), 
and  describe  the  signals  we  send  into  and  out  of  the  cavity  by  traveling  wave  combinations 
ain(a;)  and  aout(w)  of  these  transmission  line  modes.  We  will  then  extend  this  method  to 
account  for  the  leakage  of  energy  from  the  cavity  through  other  ports  as  well  to  account 
for  possible  damping  of  the  system.  In  particular,  our  model  will  yield  a  coupling  rate  71 
between  the  transmission  line  and  cavity,  and  will  assume  a  linear  loss  port  with  coupling 
72  (what  this  means  will  be  described  shortly).  If  desired,  the  reader  may  skip  to  the  main 
results  of  this  section  [(2.23)  and  (2.25),  which  make  sense  in  the  context  of  the  schematic 
in  Figure  2.1]. 

Following  [39],  we  break  up  the  total  Hamiltonian  into  the  sum 

H  =  Hk 

err  +  -f^bath  +  Hint, 


(2.11) 
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where  Tfbath  is  the  Hamiltonian  of  the  external  field  and  Hm t  describes  the  coupling  between 
the  Kerr  cavity  and  external  fields.  We  decompose  this  field  into  its  component  bosonic 
modes,  whose  annihilation  operators  we  denote  by  b(uj).  Technically,  when  an  open  port 
is  present,  the  harmonic  modes  of  the  joint  system  are  not  strictly  A  and  b(u),  but  this 
model  is  accurate  when  the  quality  factor  of  the  system  is  much  greater  than  unity  [17] 
(our  system  will  have  Q  >  100).  We  then  have  an  infinite  bath  of  harmonic  oscillators  with 
(dropping  the  zero-point  energy  for  now) 

#bath  =  J  frwb\w)b(w)du.  (2.12) 

This  continuum  decomposition  is  valid  in  the  limit  of  an  infinitely  long  transmission  line 
bath5.  We  posit  that  the  interaction  Hamiltonian  consists  of  terms  representing  the  swap¬ 
ping  of  excitations  between  the  cavity  mode  and  each  external  mode  with  a  rate  proportional 
to  a  coupling  factor6  /( w), 

/°o  r  -| 

f(u))b(u)A^  —  f*(uo)Ab\uj)  duj.  (2-13) 

-OO  L  J 

Physically,  the  lower  limit  of  integration  should  be  u  =  0,  but  since  most  of  the  weight  of 
the  integral  will  occur  near  the  resonance  frequency  of  our  system  (which  is  large  compared 
to  the  couplings  and  therefore  the  characteristic  bandwidth)  it  is  a  good  approximation  to 
extend  the  lower  limit  to  w  =  — oo.  We  proceed  first  by  solving  for  the  time  evolution  of 
b(w).  Using  the  standard  bosonic  commutation  relations 

[6(a»),  b\w')\  =  5(lu  —  J)  (2-14) 

and  the  basic  operator  identity 

[AB,  C\  =  A[B,  C\  +  [A,  C\B,  (2.15) 

we  obtain  the  Heisenberg  equation  of  motion 

bt(w)  =  j\H,  bt(w)\  =  —iubt(uj)  +  f*(u)A(t),  (2-16) 

where  the  subscript  on  bt(oj)  explicitly  denotes  time  dependence.  Following  the  approach 
of  Collett  and  Gardiner  [40,41],  we  integrate  this  ordinary  differential  equation  to  obtain 

5The  spectral  density  of  modes  is  implicit  in  our  definition  rather  than  explicitly  kept  track  of. 
6We  make  the  rotating  wave  approximation  here  to  ignore  terms  proportional  to  Ab(ui )  and 
Atbt{ui),  which  are  present  in  the  most  general  linear  coupling  between  cavity  and  bath  modes. 
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the  solution 

bt{w)  =  e-^t~to)b0(u)  +  f*(uj)  [  e-M*-0 A{t')dt'\  (2.17) 

Jt0 

where  the  first  term  represents  free  evolution  of  the  bath  and  the  second  term  represents 
waves  emitted  from  the  cavity  into  the  bath.  We  choose  to  <  t  to  be  an  arbitrary  time  in 
the  distant  past  when  the  input  signals  we  care  about  have  yet  to  be  sent  to  the  cavity.  We 
use  this  solution  to  solve  for  the  time  evolution  of  our  cavity  mode  A: 

d  A  i  i 

—  =  -  [17,  A(t)}  =  -  [HKerr,  A(t)]  ~  J  f{w)bt(u)du 

•  roo  rt  /*oo 

=  [HKeri,A(t)}~  /  /  A(t')  /  If^e-^-^dudt', 

"J  J  —  OO  j  to  J  —  oo 

(2.18) 


where  we  have  interchanged  the  time  and  frequency  integrations  in  the  final  term.  Following 
[39],  we  make  a  partial  Markov  approximation  by  treating  f(cj)  to  be  a  constant  /  over  a 
band  of  frequencies  close  to  the  characteristic  resonance  frequency  of  the  system,  which  as 
mentioned  before  is  where  the  frequency  integrals  carry  the  most  weight  and  is  valid  for 
our  system  containing  a  high-Q  coupling  capacitance.  This  allows  the  concept  of  an  ‘input’ 
and  ‘output’  mode  to  be  of  great  utility.  Namely,  we  define  the  input  mode 

Oin(t)  =  (2.19) 

which  has  the  interpretation  of  a  wave  propagating  towards  the  cavity7,  for  example  in 
a  transmission  line,  as  it  is  simply  the  free  (pre-interaction)  evolution  of  the  bath  modes 
relative  to  the  initial  conditions  at  to-  Similarly,  we  can  define  the  output  mode 

1  t"00 

flout  (t)  =  — ==  /  (2.20) 

where  ti  >  t  is  a  time  in  the  future  when  our  input  signals  are  no  longer  interacting  with  the 
cavity.  This  operator  is  simply  the  (backwards)  free  evolution  of  the  bath  from  the  distant 
future,  associated  with  a  solution  of  the  form  (2.17)  but  with  the  integration  performed 
from  t  to  t\,  and  can  therefore  be  interpreted  as  the  outgoing  signal.  We  note,  importantly, 

7Often  in  the  literature,  this  mode  is  defined  differing  by  a  minus  sign.  The  sign  is  inconsequential 
as  long  as  it  is  kept  track  of,  but  can  lead  to  errors  otherwise  (for  example,  some  of  the  typographical 
errors  in  [42]).  We  attempt  to  give  a  careful  treatment  and  avoid  sign  errors  in  this  thesis,  although 
some  mistakes  are  likely  to  have  been  made  here  as  well. 
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that  these  operators  obey  the  necessary  bosonic  algebra.  For  example, 

r  ,  i  roc  roo 

J  Z7T  J  _oo  J  —  oo 
|  roc  roo 

=  —  I  /  -  u/) 

J  —oo  J  —  oo 

=  —  /  e-^-^du  =  5(t  -  t'),  (2.21) 

2vr  J.oo 

where  we  have  used  (2.14)  and  the  inverse  Fourier  transform  of  the  Dirac  delta  distribution. 
Defining8 

7  =  7T/2,  (2.22) 

using  our  Markov  approximation  and  Fourier  analysis  we  can  rewrite  (2.18)  as 

J  A  A  roo  ft  roo 

—  =  V  [^Kerr,  A(t)]  -  f  /  e^-^b^du  -  f2  /  /  A{t')dudt' 

dt  n  J  -  oo  j-oo 

=  T  [#Kerr,  ^(t)]  “  \/27ain  ~  f  /  A(t')2Tl5(t  -  t')dt' 

n  Jt0 

=  ^[HKen,A(t)]  -  yQr/ain(t)  -  7^4 (i),  (2.23) 

a  quantum  stochastic  Langevin  equation  for  the  intra-cavity  field  in  which  the  Markovian 
approximation  has  manifested  itself  as  a  linear  damping  term  that  depends  only  on  the 
state  of  the  system  at  time  t  and  not  at  previous  times  910.  We  can  similarly  obtain  an 

8For  future  convenience,  we  uniformly  redefine  the  phase  of  the  relevant  modes  so  that  f  =  \  f\ 
is  real  within  the  band  of  the  Markov  approximation.  Also  note  that  the  definition  of  7  in  the 
literature  is  sometimes  defined  to  be  twice  as  large.  In  this  thesis,  we  will  sometimes  use  k  =  2j 
to  denote  this  (we  will  see  later  that  this  is  the  difference  between  the  full  width  at  half  maximum 
[FWHM]  and  half  width  at  half  maximum  [HWHM]  of  the  associated  resonance  peak). 

9Note  the  factor  of  1/2  that  arises  in  the  final  term  due  to  the  one-sided  integral  with  an  upper 
limit  exactly  at  the  spike  of  the  6  kernel.  This  can  be  intuitively  understood  by  thinking  of  the  d' 
kernel  as  a  limit  of  integrals  of  half  of  a  Gaussian  curve.  Although  we  have  been  somewhat  haphazard 
with  our  manipulation  of  the  Dirac  delta  distribution.  For  a  rigorous  yet  accessible  introduction  to 
the  theory  of  distributions,  see  [43] . 

10We  should  pause  and  note  that  despite  the  apparent  discrepancy,  the  units  in  this  equation  of 
motion  are  consistent.  From  the  fact  that  the  Hamiltonian  must  have  units  of  energy,  using  “[w]”  as 
a  shorthand  for  “angular  frequency”  we  see  from  (2.12)  that  b(u>)  has  units  of  [w]-1/2,  so  that  since 
A  is  unitless,  (2.13)  implies  that  /  has  units  of  [w]1/2  and  therefore  7  has  units  of  [w].  Finally,  from 
our  definition  (2.19),  ain  has  units  of  [w]1^2.  The  asymmetry  of  units  between  A  and  a;n  derives  both 
from  convention  and  the  fact  that  our  cavity  has  discrete  modes  whereas  we  model  the  transmission 
line  bath  as  a  continuum  of  modes. 
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equation  of  motion  involving  aout, 

rl  A  i 

=  ^[#Kerr,  A(t)]  -  y/&yaout(t)  +  7 A(t )•  (2-24) 

Subtracting  (2.24)  from  (2.23)  then  yields  the  main  result  of  this  section, 

a0ut(t)  -  am(t )  =  s/Zy A(t),  (2.25) 

which  together  with  (2.23)  will  allow  us  to  calculate  the  cavity  output  as  a  function  of  the 
cavity  input  and  cavity  parameters. 


2.2.2  Pumping  the  Oscillator 

We  now  have  the  tools  necessary  to  understand  the  parametric  amplification  of  inputs  to  a 
Kerr  cavity.  We  model  the  cavity  as  being  coupled  to  two  ports,  one  capacitavely  coupled  to 
an  input/output  port  and  one  to  a  linear  dissipation  port,  which  models  the  loss  of  energy 
from  the  imperfect  cavity  to  the  environment11.  We  label  the  traveling  wave  modes  of  the 
input/output  port  as  a±  and  those  of  the  linear  loss  port  as  02,  with  couplings  to  the  cavity 
of  71  and  72  respectively  as  schematically  illustrated  in  Figure  2.1. 
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j\N\r> 


7i 


Transmission  Line 
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<(t) 


Nonlinear  Resonator 


HKeiI  =  hw0A^A+h-{Al)2A2 


72 


a. 


2  W 

ANW 


Linear  Loss 


iVW> 


a 


out 


(*) 


Figure  2.1  Pumping  a  nonlinear  Kerr  cavity  with  a  single  “pump”  tone.  71  is  the  cou¬ 
pling  strength  associated  with  the  port  connecting  to  the  main  input/output  transmission 
line,  while  72  represents  a  coupling  strength  associated  with  an  effective  port  representing 
linear  loss.  The  intra-cavity  held  obeys  the  bare  Kerr  Hamiltonian  in  the  absence  of  these 
couplings. 

11  More  generally,  there  could  be  a  nonlinear  dissipation  port  through  which  a  bath  couples  to  A2 
rather  than  A.  Such  ports  are  usually  negligible,  but  we  will  later  quote  a  modification  of  our  result 
to  account  for  this. 
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Upon  evaluating  the  commutator 

[HKe„,A\  =  hu0[A^A,A\  +  £y[(At)2A2,A]  =  hu0[A\  A\A  +  A\  A\AA 

=  -hu0A  -  hKA]AA  (2.26) 

the  Langevin  equation  (2.23)  with  two  independent  sets  of  coupled  bath  modes  yields 

(j  A 

—  =  -iu0A  -  iKA^AA  -  -  (71  +  72 )A,  (2.27) 

and  (2.25)  yields  the  boundary  conditions 

a°ut(t)  -  a™(t)  =  1/271  A(t),  (2.28a) 

a°2 ut(t)  -  a%(t)  =  y/^A(t).  (2.28b) 

To  utilize  the  cavity  as  a  parametric  amplifier,  we  need  to  input  a  strong  classical  pump 
tone,  which  will  provide  the  energy  (and  oscillating  parameter)  for  the  parametric  gain  of 
small  input  signals.  That  is,  the  amplification  scheme  involves  inputting  a  pump  plus  small 
signal  perturbation, 

in  in,p  ,  in,s  /r>  nn\ 

a,  =  al  H  +  al  .  (2.29) 

We  will  follow  the  analysis  of  [21]  by  first  calculating  the  cavity  response  to  the  classical 
pump,  and  afterwards  treating  the  signal  as  a  perturbation.  We  write  the  incoming  pump 
(assumed  to  be  in  the  coherent  limit12)  as 

=  bfe~i^pt+^\  (2.30) 

where  b™  is  a  real  constant,  and  the  reflected  signal  as 

aout,P  =  6out  (2.31) 

where  b°ut  is  in  general  complex.  Analogously  to  the  case  of  a  classical  Duffing  oscillator  [34], 

12As  studied  in  quantum  optics,  a  coherent  state  is  the  classical  limit  of  a  sinusoidal  drive,  consist¬ 
ing  of  an  infinite  sum  of  all  Fock  number  states.  The  amplitude-phase  uncertainty  product  of  such  a 
state  is  the  minimum  possible  for  the  harmonic  oscillator.  The  expectation  value  rotates  clockwise 
about  the  origin  of  the  Re  [a],  Im[a]  plane  and  the  “wavepacket”  in  this  phase  space  is  very  localized 
for  mean  photon  number  n  1,  allowing  the  state  to  be  treated  as  a  complex  number  proportional 
to  e~lu,t. 
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we  posit13  a  steady-state  zeroth-order  response  of  A  in  the  form 


—  Pg-iOpt+fiB) 


(2.32) 


where  B  is  real  and  4>b  encodes  the  internal  cavity  phase  response.  Note  that  we  are 
treating  the  pump  mode  solution  operators  as  commuting  c-numbers,  when  in  reality  the 
steady  state  response  is  something  more  like  a  harmonic  oscillator  coherent  state.  Plugging 
these  three  equations  into  (2.27)  and  (2.28a)  in  the  absence  of  signal  and  noise,  we  obtain 
respectively 

[i{u o  -  up)  +  71+72 ]B  +  iKB 3  =  - (2.33) 

and 

6°ut  =  bf  +  .  (2.34) 

In  the  special  case  of  a  linear  cavity  with  no  Kerr  term  (or  for  sufficiently  low  drive 
powers  with  a  nonlinear  cavity,  such  as  the  black  curves  in  Figure  2.2),  we  can  analytically 
solve  for  the  output  in  terms  of  the  input  by  combining  (2.33)  and  (2.34)  with  K  =  0, 
obtaining  a  reflection  coefficient  of 


=  1 _ 271 

b f  i(u0  -  up)  +  7i  +72 

_  To  -  2 
" 1 ’ 


72  -  71  +  -  Up) 

7i  +  72  +  i(u0  -  up) 


(2.35) 


where  5u  =  up  —  uq,  Q  is  the  loaded  quality  factor 


1  _  1  |  1  _  27i  |  272 

Q  Qex  t  Q  int  Uq  Uq 


(2.36) 


13We  are  ignoring  the  possibility  exciting  higher  harmonics  characteristic  of  the  response  of  non¬ 
linear  oscillators  -  in  our  fundamental  assumption  of  the  Kerr  Hamiltonian  (2.10),  we  have  ignored 
potential  cross-Kerr  terms  (of  the  form  hKtj  A\  A% A)  A:j )  that  may  weakly  couple  to  higher  modes 
of  the  cavity,  and  also  in  our  rotating  wave  approximation  we  have  ignored  a  small  third  harmonic 
contribution  to  the  pump  solution.  This  is  ignored  in  [15,21],  although  it  is  usually  discussed  in 
treatments  of  the  classical  Duffing  oscillator  [34, 42] .  We  should  make  sure  that  our  single- harmonic 
zeroth-order  pump  response  is  accurate.  To  do  so,  note  that  adding  a  term  of  order  h^-A3  to  (2.27), 
plugging  in  the  zeroth-order  solution  (2.32),  and  then  integrating  to  obtain  a  more  accurate  solution 
yields  a  term  of  order  K^3e-*(3wpt+3<jiB)^  which  wiH  be  small  compared  to  the  zeroth-order  solution 
when  -j-B2  <C  1.  As  we  will  see,  in  our  implementation  K/up  ~  10~6  to  10— 5 ,  and  we  will  op¬ 
erate  at  photon  numbers  B 2  less  than  the  so-called  critical  number  (2.46)  of  order  103  in  a  cQED 
implementation.  We  can  therefore  treat  the  zeroth-order  classical  pump  solution  (2.32)  as  accurate. 
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and 


To  = 


72  -  7i 


(2.37) 


7i  +  72 

is  the  reflection  coefficient  on  resonance.  Note  that  being  on  resonance  yields  the  minimum 
amount  of  reflected  signal,  as 

7+4Q2(f|)2 


r  = 


1  +  4  Q2  S 


2  ’ 


(2.38) 


In  particular,  in  the  limit  71  =72,  we  have  To  =  0  indicating  full  leakage  of  the  input 
signal  through  the  linear  loss  port,  analogous  to  the  full  transmission  of  a  standard  c.QED 
two-port  A/2  resonator  on  resonance.  But  in  the  limit  71  3>  72  (which  is  the  case  for  our 
par  amp),  To  ~  1,  meaning  that  nearly  all  information  about  the  resonance  is  contained  in 
the  reflected  phase  rather  than  in  the  reflected  amplitude.  In  either  case,  T  — >  1  as  ujp  is 
detuned  far  from  ujq  relative  to  the  linewidth,  71  +72. 

Returning  to  the  general  case  of  K  A  0;  we  follow  [21]  by  taking  the  magnitude  of  both 
sides  of  (2.33)  and  defining  the  cavity  pump  photon  number 


TV  =  B2, 


(2.39) 


obtaining  the  following  cubic  equation  for  TV, 

K2N3  +  2(cn0  -  ujp)KN2  +  [(Wo  -  top)2  +  (71  +  72)2]TV  =  27l  (bf)2  .  (2.40) 

We  can  also  obtain  an  equation  for  the  phase  of  the  steady-state  internal  pump  field  relative 
to  the  input  pump  field  using  equation  (2.33)  and  Euler’s  formula, 

tan {(j>B  -  (f>p)  =  —  +  A^V'> — — .  (2-41) 

71+72 

Note  that  this  phase  response  is  nonlinear,  depending  the  internal  field  amplitude  TV.  We 
have  grouped  this  TV-dependent  term  together  with  loq  so  that  it  can  be  interpreted  as  a 
power-dependent  shift  KN  of  the  resonance  frequency  of  the  cavity.  For  negative  K  (which 
as  we  will  see  is  the  case  in  our  implementation),  the  cavity  resonance  is  pulled  downward 
in  frequency  with  increasing  pump  power.  The  expression  I\N  for  this  shift  also  coincides 
with  the  point  of  the  maximum  intracavity  response  TV,  as  is  easily  verified  by  solving  for 
dN/dujp  =  0  in  (2.40).  Equation  (2.40)  can  be  numerically  solved  for  the  steady-state 
amplitude  of  the  internal  pump  field.  There  are  two  possible  regimes  for  this  polynomial 
equation;  either  there  exist  three  real  solutions  for  TV,  or  there  exist  one  real  solution  and 
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Figure  2.2  MATLAB  plots  of  system  response  to  a  single  pump  input  tone  as  a  func¬ 
tion  of  frequency  (in  units  of  linewidths)  with  K  =  — 10_5wo  and  72  =  .I71,  for  drive 
amplitudes  b™  equal  to  ■01611“crit  (black),  .56I1ncrit  (red),  6^^  (green),  and  1.56jncrit  (blue). 
For  strong  drives  6"1  >  b\ncrit ,  dashed  lines  in  the  bistable  regime  indicate  the  unstable 
solution,  (a)  Intracavity  pump  response  Nnoim  [N  normalized  by  ^7l+72^  (bjn)2].  The 
actual  (unnormalized)  response  of  course  increases  in  maximum  amplitude  in  addition  to 
having  a  distorted  shape,  (b)  Intracavity  phase  response,  (c)  Magnitude  of  the  reflection 
coefficient,  (d)  Phase  of  the  reflection  coefficient. 


two  complex  solutions.  In  the  latter  case,  we  take  the  real  solution  to  be  the  physical  one. 
In  the  former  case,  it  can  be  shown  that  two  of  the  solutions  will  be  stable,  while  one 
will  be  unstable,  leading  to  the  feature  of  bistability  [21].  We  can  obtain  some  intuition 
behind  this  by  plotting  the  behavior  of  the  pump  solution  at  different  strengths  of  incident 
b™  as  a  function  of  frequency.  The  response  is  that  of  a  driven  Duffing  oscillator  (which 
we  have  arrived  at  in  the  Fourier  domain).  Figure  2.2  contains  MATLAB  plots14  of  the 
theoretical  pump  solutions  we  have  just  calculated.  The  left  column  depicts  the  steady-state 


14The  MATLAB  function  “roots”  is  quite  useful  in  solving  for  the  pump  response  in  all  regimes. 
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intracavity  amplitude  and  phase,  while  the  right  column  depicts  the  amplitude  and  phase 
of  the  reflected  signal.  Note  that  the  system  resonance  is  pulled  downward  in  frequency,  as 
expected  for  a  negative  value  of  K.  We  would  like  to  operate  our  parametric  amplifier  with 
a  pump  power  slightly  below  the  critical  input  power,  so  that  the  system  will  be  sensitive  to 
small  perturbations  but  not  yet  bistable.  The  bistable  regime  in  which  there  are  two  stable 
solutions  occurs  above  this  bifurcation  point.  Although  we  will  not  operate  at  such  powers, 
this  bistable  regime  has  been  successfully  used  to  create  sensitive  “Josephson  bifurcation 
amplifiers”  [20, 44-46] . 

We  can  calculate  the  critical  point  at  which  the  distorted  Lorentzian  response  transitions 
into  the  bistable  regime;  i.e.,  when  the  slope  of  the  response  first  becomes  vertical  as  the 
system  transitions  into  the  bistable  regime,  as  in  the  green  curves  in  Figure  2.2.  There  are 
two  conditions  that  are  satisfied  at  this  point.  One  is  that  the  slope  d N/dup  diverges,  or 
rather  dujp/dN  =  0.  Well,  differentiating  (2.40)  with  respect  to  N  and  collecting  the  terms 
involving  d up/dN  yields 

duj 

[2 KN2  +  2 N{0jo  -  Up)\  =  3K2N2  +  4KN{u0  -  u>p)  +  (w0  -  wp)2  +  (7i  +  72)2,  (2.42) 

which  means  (completing  the  square  on  the  right-hand  side)  that  the  critical  point  occurs 
when 

(w0  -  ujp  +  2 KN)2  -  I<2N 2  +  (7i  +  72)2  =  0.  (2.43) 

In  the  bistable  regime,  this  will  be  satisfied  at  two  pump  frequencies.  The  critical  point 
occurs  when  these  two  points  coincide,  which  will  happen  when  in  addition  d2uip/dN 2  =  0. 
Differentiating  both  sides  of  (2.42)  using  the  product  rule  and  setting  dup/dN  to  zero,  we 
obtain  the  additional  condition 


6K2N  +  4K(ujq  —  up)  =  0. 


(2.44) 


Combining  these  equations  then  yields  the  critical  values  of  ojp  and  N, 


I< 


u. 


',crit  =  W0  +  \/3  — (7l  +  72) 


K 


(2.45) 


and 


Amt  = 


2  71  +  72 


(2.46) 


v/3  \K\  • 

We  will  also  find  it  useful  to  calculate  [using  (2.40)  along  with  (2.45)  and  (2.46)]  the  input 
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power  corresponding  to  the  critical  intracavity  field, 

(2 

2.2.3  Parametric  Amplification  of  Small  Perturbations 

We  are  finally  ready  to  calculate  the  parametric  amplification  of  small  signals.  Intuitively, 
if  the  pump  is  operated  near  the  critical  point,  small  perturbations  will  lead  to  a  relatively 
large  change  in  the  system’s  response.  Suppose  a  small  signal  apb  =  c\(t)me~lUpt ,  with 
|c\n|  <C  |6f|.  Here,  c™  can  be  time-dependent  -  the  factor  of  e~lU}pt  is  there  simply  for 
calculational  convenience,  allowing  us  to  “peel  off’  the  time-dependence  of  the  pump  and 
define  c  in  a  rotating  frame  (as  we  will  explicitly  see  when  we  move  into  the  Fourier  domain). 
Our  input  mode  now  takes  the  form 


(hin  N2  _  4  (71  +  72)3 

1  1>crit;  3V3  7iA| 


ali(t)  =  e-i  Kt+^)  +  c^(t)e~iujpt, 

with  the  output  and  internal  resonator  modes  taking  the  form 

a°ut  (t)  =  +  cT\t)e-iuJpt 

A(t)  =  Be~i{uipt+ct,B)  +  c(t)e~iUpt, 


(2.48) 


(2.49) 


where  b°ut  and  B  are  the  solutions  to  the  classical  unperturbed  pump  as  derived  above. 
This  decomposition  of  input  and  output  signals  is  illustrated  schematically  in  Fig.  2.3.  We 
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Figure  2.3  Illustration  of  input,  output,  and  loss  for  a  parametric  amplifier  based  on  a 
nonlinear  Kerr  oscillator.  The  main  input-output  port  connecting  to  a  transmission  line 
is  shown  on  the  left,  while  an  extra  port  representing  linear  loss  to  the  environment  (as 
described  above)  is  shown  on  the  right.  The  decomposition  of  the  input/output  signals 
into  a  large  classical  pump  and  small  quantum  signal  is  also  illustrated. 
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assume  that  perturbations  to  the  pump  are  small  so  that  we  can  treat  c11n,c°ut,  and  c  as 
perturbations  to  linear  order,  which  can  include  both  a  small  signal  and  also  noise.  To 
derive  the  behavior  of  the  resonator  as  an  amplifier  of  input  signals,  we  will  assume  the 
input  from  the  loss  mode  a™  is  comparably  small,  so  that  plugging  (2.48)  and  (2.49)  into 
the  Langevin  equation  of  motion  (2.27)  and  subtracting  out  the  pump  response  from  both 
sides  [which  still  approximately  obeys  (2.33)]  yields,  after  dropping  all  quadratic  and  cubic 
appearances  of  c(t)  (which  is  justified  when  (cU)  <C  N),  the  linearized  equation 

-iupc(t)  +  ^  =  -iu0c(t)  -  iK  (cf(t)Ne~2l(/>B  +  2c(t)N^j  -  V^Ti cT(t)  -  (7i  +  72 )c{t). 

(2.50) 

Note  that  peeling  off  the  e~lUpt  term  from  the  signal  allowed  us  to  divide  out  a  common 
factor  of  e~lUpt  in  this  equation.  We  then  define  the  frequency-domain  signals 

c)n(w)  =  -=/  dtc^(t)elut 
i  r°° 

c?ut(w)  =  -f=  /  dtc°^\t)elu:t 
4-oo 

1  f°° 

C (w)  =  -w=  /  dtc(t)elujt  (2.51) 

V  J — oo 

in  the  rotating  frame  of  the  pump.  It  is  important  to  remember  this:  because  of  how 
we  have  defined  the  signal  term  in  (2.49),  c o  refers  to  the  frequency  relative  to  the  pump 
frequency,  that  is,  the  physical  signal  has  a  frequency  of  up  +  cu.  Using  the  basic  properties 
of  the  Fourier  transform  (in  particular,  how  it  acts  on  time-derivatives  and  under  complex 
conjugation),  we  obtain  after  Fourier-transforming  and  rearranging  (2.50) 

i  [(too  —  wp)  —  lo  —  i( 7i  +  72)  +  2 KN]  c(w)  +  iK Ne~2l^B c1  (—u)  =  —  y/27ic11n(u;).  (2.52) 

Defining 

W  =  i(uj0  -  Up)  +  (71  +  72)  +  2 iKN  (2.53) 

V  =  iKNe-2i4>B,  (2.54) 

(2.52)  can  be  written  as 


(W  —  iu)c{u)  +  Vc\—u)  =  —  y/27ic11n(w). 


(2.55) 
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As  is  often  the  case  with  working  in  the  frequency  domain,  a  little  algebraic  manipulation 
can  allow  us  to  solve  for  the  variable  under  question.  In  our  case,  we  can  solve  for  the 
perturbed  internal  cavity  field  c(u)  by  using  the  complex  conjugate  version  of  (2.55)  with 
oj  — y  — oj , 

(IF*  -  iu)c\-u )  +  V*c(u)  =  -v/2^c‘nt(— W),  (2.56) 

combined  with  (2.55)  to  yield  the  combined  equation 

-  V*c(u ))  _  . 

(W  -  iu)c(u)  +  F^ - —  _  — - =  -y/2^c?(u),  (2.57) 


which  after  collecting  terms  yields  the  solution 


-\/2t 7 

(IF*  -  iu)cf(u)  -  Vc'i\—u) 

\w\ 

2  —  2  Re[VF]wn  —  u2  —  F  2 

(2.58) 


We  will  find  it  instructive  to  factor  the  denominator,  which  is  simply  a  quadratic  polynomial 
in  the  variable  x  =  iu.  The  quadratic  formula  then  gives  the  roots 


A±  =  (71  +  72)  ±  \J I\2N2  -  (w0  -  wp  +  2KN)2, 


(2.59) 


where  we  have  used  the  definitions  (2.53)  and  (2.54).  Thus,  (2.58)  becomes 


c(w) 


?'w)c11n(w)  — 
\-)(iu  —  A+) 


(2.60) 


Finally,  using  this  intra-cavity  signal  solution,  we  can  calculate  the  cavity  output  c°ut  [as 
defined  in  (2.49)]  in  response  to  the  perturbation  via  the  input-output  boundary  condition 
(2.28a), 


c?utM 


=  1- 


271  (IF*  —  iu>) 


(iu  —  X-)(iu  —  A+) 

=  Q(u)&i(u)  +  M(u)c1^  (-u). 


c?(co)  + 


271 V 


dnt 


(iu  —  X-)(iu  —  A_(_) 


(-u) 


(2.61) 


This  result  shows  something  very  interesting:  if  the  signal  is  incident  on  a  pumped  Kerr 
resonator  with  detuning  u  from  the  pump,  (2.61)  implies  that  the  reflected  output  will 
contain  not  only  a  frequency  component  at  up+u,  but  also  at  up  —  u.  In  effect,  the  pumped 
nonlinear  resonator  couples  the  two  symmetrically  detuned  modes  via  a  mixing  process. 
Note  that  the  output  signal  at  frequency  up  —  u  only  occurs  when  a  Kerr  nonlinearity  is 
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present;  i.e.,  K  ^  0  — >■  V  A  0.  We  will  refer  to  this  tone  as  the  “image”  tone  (or  sometimes 
as  the  “intermodulation”  tone). 

Analogy  to  Classical  Parametric  Resonance 

Before  delving  into  the  properties  of  the  amplification  we  have  just  derived,  we  take  a 
brief  but  enlightening  aside  in  order  to  understand  the  correspondence  between  our  result 
and  conventional  parametric  resonance  in  the  classical  limit.  In  quantum  optics,  the  term 
“parametric”  generally  refers  to  processes  in  which  different  frequencies  of  light  are  mixed 
via  a  nonlinear  refractive  medium  [15,39].  But  if  we  put  in  a  little  effort,  we  can  in  fact  see 
how  our  parametric  amplifier  formally  maps  onto  a  conventional  classical  driven  parametric 
oscillator  (2.6);  i.e.,  an  oscillator  with  a  time-varying  resonance  frequency  of  the  form  (2.3) 
capable  of  exhibiting  parametric  resonance. 

We  can  do  this  by  temporarily  returning  to  the  time-domain  (2.27)  and  temporarily 
considering  the  intracavity  signal  perturbation  5(t)  =  c(t)e~lUpt  in  the  laboratory  frame. 
We  will  show  that  the  equation  of  motion  for  S(t )  maps  onto  equation  (2.6)  in  a  certain 
classical  limit.  As  in  our  derivation  of  the  rotating- wave  perturbation  equation  of  motion 
(2.50),  inserting  the  pump  plus  signal  ansatz  into  (2.27)  and  subtracting  out  the  pump 
response  from  both  sides  yields 

8(t)  =  -iu05{t)  -  iK[N5\t)e~2iu}pte~2i^B  +  2 NS(t)}  -  y/2 -  (71  +  72 )S(t).  (2.62) 

Differentiating  and  and  then  substituting  this  equation  yields  [using  the  shorthand  defini¬ 
tions  (2.53)  and  (2.54)] 

S(t)  +  (W  +  iu}p)S(t)  +  e~2i“ptV${t)  -  2iujpV5\t)  = 

=7< S(t )  +  (W  +  iup)S(t)  +  V[-y/2 rf(t)  -  W*5\t)  -  D*<S(f)]  =  -^2^8?(t) 

=> 5{t )  +  (W  +  iup)8(t )  +  e~2iUptV[-^ -  (W*  -  iup)tf(t)  -  e2iUptV*S(t)} 

-  2 iupe~2iUptV5\t)  =  -y/¥^8f(t) 

=> S(t )  +  (W  +  iojp)5(t)  -  e-2^V^;ntW  -  (W*  -  -  S(t ) 

-(W  +  iup)5{t)}  -  \V\25(t)  -  2 iupe~2iuJptV5\t)  =  -^[^{t) 
=>8(t)  +  2Re[W  +  iup]6(t)  +  \\W  +  iup\2  -  |D|2  +  2 KNupe~2iujpt-2^B]  S\t)  =  T(t), 

(2.63) 

where  T(t)  is  an  effective  driving  term  for  the  intracavity  perturbation  field.  If  we  make  the 
classical  approximation  that  5(t)  is  a  c- number  and  in  fact  real,  then  taking  the  real  part 
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of  this  equation  of  motion  and  using  the  definitions  of  V  and  W  yields 

^  +  2(71  +72)^  +uls{N,t)6  =  r  (t),  (2.64) 

where,  assuming  cop/uj 0  ~  1  and  K2 /wq  <  1  we  have 


wefr(^r  1 1)  —  w0  1  +  (71  +  72 )  + 


ANK  \ 
wo  / 


1  + 


2NK/u0 


1  +  (71  +  72) 2  +  ANK/uj0 


cos  (2  ojpt  +  2(f)  b) 


(2.65) 


This  is  precisely  the  form  (2.6),  with  Uq  (1  +  (71  +  72)2  +  ANK/ujq)  being  the  shifted  nat¬ 
ural  frequency  of  a  damped  nonlinear  oscillator  as  in  (2.3)!  [The  damping  arises  from  the 
coupling  of  the  resonator  through  the  input  and  loss  ports  71  and  72.]  We  can  therefore 
say  that  in  some  sense,  the  strong  pump  acts  to  modulate  the  effective  resonance  frequency 
seen  by  a  small  perturbation,  which  is  possible  because  of  the  nonlinear  Kerr  term.  This 
provides  energy  to  the  perturbation  when  it  is  near  parametric  resonance.  Along  with  our 
discussion  in  section  2.1,  this  gives  a  bit  of  intuition  (however  small)  behind  the  amplifier 
gains  we  will  now  study. 


Amplifier  Gain  and  Bandwidth 


Returning  to  our  quantum  mechanical  result  (2.61),  we  would  like  to  understand  the  prop¬ 
erties  of  our  parametric  amplifier  based  on  a  quantum  Kerr  oscillator.  The  bare  minimum 
way  of  characterizing  an  amplifier  is  of  course  calculating  its  gain  and  bandwidth,  which 
we  will  do  presently.  After  this,  we  will  characterize  other  important  properties  such  as 
dynamic  range  and  noise  performance.  We  define  the  signal  and  image  power  gains 


~out , 


G» 


7  = 


_  14 


w 


K»|5 


=  |£(w)|2  = 


\(ioj  —  A —  A_|_)  —  27i(W  —  iuj) 
(o;2  +  A2)(w2  +  A2_) 


(2.66) 


and 


Gi{oS) 


\M(u;)\2 


47i2m2 

(o;2  +  A2)(w2  +  A2_) 


(2.67) 


respectively.  We  see  that  the  gains  are  highest  when  the  signal  frequency  is  close  to  the 
pump  frequency,  00  ~  0,  and  also  when  the  cavity  is  driven  close  to  its  critical  bifurcation 
point  [at  which  point  the  denominator  diverges,  as  can  be  seen  from  equation  (2.43)].  As 
noted  in  [15,42],  the  3-dB  bandwidth  of  this  peak  is  of  the  order  of  the  linewidths  71,72- 
Following  [42],  we  can  straightforwardly  estimate  this  bandwidth,  assuming  the  cavity  is 
pumped  close  to  (but  not  at)  the  bifurcation  point  discussed  earlier.  For  oj  «  0,  we  can 
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approximate  (2.66)  and  (2.67)  as 

Gs( u  ^0)^  Gi(uj  «  0)  =  G(u  »  0) 


47?  1^1 2 


_  G(0) 


+  +  l  +  (^)2’ 

from  which  we  can  extract  the  3-dB  bandwidth 


(2.68) 


_  2  /  47l2|V|2 

\/G(0)  \l  x2_+x2+ 


(2.69) 


Using  equations  (2.45)  and  (2.46)  for  the  values  of  various  parameters  near  the  critical 
pumping  point,  we  have  A+  ~  2(71  +  272)  and  A_  ~  0,  so  that  (2.69)  becomes  in  the 
near-critical  limit 


2  /  hi  ■ 

(7(0)  V 
4  71 

\/3  y/G(0)' 


(2.70) 


We  see,  as  is  usually  the  case  with  amplifiers,  that  there  is  a  tradeoff  between  gain  and 
bandwidth  [in  fact,  we  also  see  this  in  the  classical  case  of  parametric  amplification  via 
equation  (2.8)].  This  key  result  states  more  specifically  that  the  gain-bandwidth  product 
y/GBu  is  proportional  to  71,  the  rate  at  which  photons  can  leak  out  of  the  output  port  of 
the  cavity  (it  is  interesting  to  note  that  although  72  will  broaden  the  resonance  peak  of  the 
pump,  it  doesn’t  affect  the  bandwidth  of  amplification  under  our  approximations).  In  other 
words,  the  bandwidth  is  roughly  given  by  that  of  a  Lorentzian  with  the  same  linewidth  as 
the  resonator,  divided  by  \[G.  This  is  an  important  reason  why  implementing  a  parametric 
amplifier  with  SQUIDs  is  so  appealing:  as  we  will  explicitly  see  later,  tuning  the  frequencies 
of  the  SQUIDs  embedded  in  the  centerpin  of  a  A/4  resonator  via  magnetic  flux  allows  one 
to  tune  the  resonance  frequency  of  the  system,  and  therefore  the  amplified  band.  This  is 
one  approach  to  overcoming  the  small  bandwidths  that  are  characteristic  of  a  large  class 
of  these  amplifiers.  Figure  2.4  shows  MATLAB  plots  of  the  calculated  gains,  in  which  the 
gain-bandwidth  tradeoff  is  visually  apparent. 


Dynamic  Range 

Even  if  we  increase  the  cavity  linewidth  and  compromise  by  decreasing  the  bandwidth,  we 
cannot  increase  the  gain  indefinitely.  The  dynamic  range  of  an  amplifier  refers  to  the  range 
of  amplified  signal  powers  for  which  the  device  functions  well  as  a  linear  amplifier.  In  our 
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Figure  2.4  Theoretical  signal  and  image  gains  (in  logarithmic  power  units  of  decibels 
[dB])  with  K  =  — 10_5wo  and  72  =  .I71,  for  pump  amplitudes  b “  equal  to  .•5611ncrit  (black), 
•76i"crit  (purple),  -9&incrit  (red),  .956^  (green),  and  .996^  (blue).  In  (a)  and  (b), 
a  signal  frequency  of  w  =  0  relative  to  the  pump  was  used.  In  (c)  and  (d),  the  pump 
frequency  for  each  drive  power  was  chosen  based  on  what  gave  the  maximum  gain  in  (a) 
and  (b). 


previous  derivation,  we  assumed  that  the  pump  tone  is  “stiff’  in  that  it  is  unchanged  by 
any  perturbing  signals  or  quantum  fluctuations,  whereas  in  reality  the  pump  is  ultimately 
what  supplies  the  energy  for  signal  amplification  in  the  first  place  and  could  in  principle  be 
depleted.  This  not  only  imposes  a  limit  on  the  power  of  the  input  signal,  but  also  imposes  a 
limit  on  the  power  of  the  output  signal  and  therefore  on  the  dynamic  range  of  the  amplifier. 
At  the  very  least,  the  amplifier  must  be  able  to  function  in  the  face  of  noise  due  to  both 
quantum  and  thermal  fluctuations  at  the  input  port.  To  be  safe,  heuristically  we  would 
want  this  amplified  noise  to  be  two  orders  of  magnitude  smaller  than  the  pump  power.  We 
can  roughly  estimate  the  amplified  vacuum  noise  from  frequencies  near  the  pump  frequency 
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(ignoring  72)  as 

Pnoise  «  klJlP  [  G(0)df,  (2.71) 

1  -h  Bu 

which  for  ujp/2it  =  6  GHz,  a  power  gain  of  100  (i.e.,  20  dB),  and  7i/27t  ~  15  MHz  yields  a 
rough  value  of  —120  dBm,  meaning  that  the  critical  power  of  the  pump  should  at  least  be 
on  the  order  of  —100  dBm  for  such  operating  parameters.  With  this  satisfied,  it  is  then  up 
to  the  user  of  the  device  to  ensure  that  the  signal  input  power  is  not  an  appreciable  fraction 
of  the  pump  power.  In  particular,  when  designing  the  amplifier,  one  should  ensure  that  the 
critical  pump  input  power  -Pinout  =  ficcp(b™CTit)2  is  greater  than  this  value.  With  these  same 
parameters,  by  (2.47)  this  would  require  the  Kerr  constant  \K\/2tt  to  be  less  than  on  the 
order  of  100  KHz. 

Quantum-Limited  Noise  Performance 

In  this  section,  we  will  arrive  at  a  very  important  property  of  a  parametric  amplifier  based 
on  a  quantum  Kerr  oscillator  -  at  zero  temperature,  its  added  noise  output  is  limited 
by  quantum  fluctuations  to  the  minimum  amount  allowed  by  quantum  mechanics.  More 
precisely,  besides  just  amplifying  quantum  (or  classical)  fluctuations  at  the  input  as  briefly 
discussed  above,  any  linear  amplifier  is  required  by  quantum  mechanics  to  add  a  certain 
minimum  amount  of  extra  noise  to  its  output.  We  will  show  that  the  parametric  amplifier 
theoretically  reaches  this  limit.  We  first  (non-rigorously)  introduce  the  concept  of  a  noise 
spectral  density,  which  characterizes  how  much  of  a  noise  signal  is  contained  in  a  region  of 
frequency  space.  Intuitively,  this  number  would  be  proportional  to  the  measured  output 
power  if  the  output  port  was  followed  by  a  narrow  bandpass  filter  centered  on  the  frequency 
in  question.  Classically,  the  energy  carried  by  a  signal  (such  as  a  voltage)  x(t)  in  the 
time  domain  is  proportional  to  f  \x(t)\2dt.  By  Parseval’s  theorem  [47],  this  is  equal  to 
(modulo  a  factor  of  27r  depending  on  the  convention  used)  1/27 r  f  \x(u)\2dt,  where  x(oj) 
denotes  the  Fourier  transform  of  the  signal.  By  the  orthogonality  of  Fourier  components,  it 
makes  sense  to  divide  the  total  power  into  independent  contributions  from  each  frequency, 
which  could  be  physically  measured  for  example  by  a  spectrum  analyzer.  Classically,  we 
therefore  identify  the  spectral  density  Sx,x  with  |x|2,  modulo  a  constant  factor.  Quantum 
mechanically,  we  can  deduce  a  similar  structure  due  to  the  orthogonality  of  (the  continuum 
of)  bath  modes  (although,  the  quantum  mechanical  spectral  density  is  not  in  general  real, 
making  its  interpretation  a  little  more  intricate).  Namely,  the  energy  of  a  given  mode  a(oj) 
is  a  multiple  of  h(up  +  u)  determined  by  the  thermal  occupation.  For  bosonic  modes  each 
obeying  a  harmonic  oscillator  Hamiltonian,  we  use  the  Bose-Einstein  distribution  [48]  to 
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deduce  the  symmetrized  spectral  density15  in  the  approximation  inp  +  in  ~  in. 


tunpSa,a{u)  =  tvjjp  la\in)a{in)  +  a(in)a\in)\  =  hinp  /  a\in)a( in)  +  - 


=  tun, 


1 _ I  =  cotu  (  \ 

P  \  ehu]p/kBT  _  1  +  2J  2  '  \2kBT )  ' 


(2.72) 


We  can  now  calculate  the  added  fluctuation  noise  at  the  output  of  our  amplifier.  Using 
(2.61),  we  have 


Sa,out,aout(Lu)  —  Gs(uj)Sain,ain(w)  +  Gj(ca)  S'ain,ain(  w), 


(2.73) 


which  shows  that  for  a  signal  at  in.  extra  noise  at  the  image  frequency  —  in  will  be  generated 
in  accordance  with  the  general  mixing  properties  of  the  parametric  amplifier.  Now,  it  can 
be  shown  that  the  signal  and  intermodulation  gains  of  (2.66-2.67)  satisfy  [45] 


Gs{in )  —  Gi{in)  +  1. 


(2.74) 


The  noise  spectrum  (2.73)  then  becomes 


•Saout^out  —  Gs(m)Sain,ain{u)  +  [Gs(w)  l]<S,ainj0in(  ^) 

1 


=  Gs(u>) 


Sain,ain(u)  +  1  - 


Gs(u>) 


'S'ain.ainC  w) 


(2.75) 


So,  interpreting  the  signal  at  in  as  the  signal  and  Saouttaout(in)  as  the  input  noise,  the  added 
output  noise  will  arise  in  the  image  channel,  so  that  using  (2.72)  we  obtain  an  added 
output  noise  spectral  density  for  frequencies  near  the  pump  (w  <  71),  in  units  of  tunp  and 
normalized  by  the  gain, 


An 


1 

GM 


-  coth 


(2.76) 


Note  that  in  the  limit  kBT  <C  tunp,  for  large  gains  this  approaches 

-l.Y  ^  2:  (2‘77) 

15We  can  intuitively  view  this  as  the  sum  of  the  absorption  and  emission  spectrum  (respectively) 
of  the  mode  of  frequency  in,  which  is  also  equal  to  to  the  expectation  value  of  fluctuations  of  the 
harmonic  oscillator  energy  as  shown  by  the  second  equality.  For  a  more  rigorous  treatment  of 
quantum  spectral  densities  based  on  correlation  functions  (which  in  quantum  mechanics  can  be 
imaginary),  see  [49]. 


2.2  Parametric  Amplification  with  a  Quantum  Duffing  Oscillator 


30 


meaning  that  for  large  gains  the  parametric  amplifier  will  always  add  half  a  quantum  of 
noise  (i.e.,  half  a  photon  of  energy  hujp)  at  its  output.  Conversely,  in  the  classical  limit 
ksT  3>  frwp,  we  obtain  the  classical  result  A n  oc  &bT. 

We  now  show  that  this  is  in  fact  the  best  we  can  do;  that  is,  this  is  the  minimum 
amount  of  added  noise  allowed  by  quantum  mechanics  for  any  kind  of  phase-preserving, 
phase- insensitive  amplifier  (these  terms  will  be  defined  shortly).  This  is  a  key  reason  why 
Josephson  parametric  amplifiers  are  so  appealing  in  circuit  QED  experiments  -  they  can 
be  used  to  amplify  signals  such  that  quantum  noise  becomes  even  more  important  than  the 
otherwise  significant  noise  at  later  stages  of  commercial  HEMT-based  amplification.  This 
fundamental  limit  is  partially  a  consequence  of  the  uncertainty  principle  for  Hermitian  op¬ 
erators  on  a  Hilbert  space  [47,50]  (along  with  the  identification  of  physical  observables  with 
Hermitian  operators  in  quantum  theory):  if  A  and  B  represent  two  self-adjoint  observables, 
then 

*a4>(^<[AB]>)  •  (2-78) 

To  utilize  this  relation,  we  must  connect  the  language  of  quantum  mechanics  to  familiar 
laboratory  microwave  voltage  signals.  If  the  frequency  of  the  input  signal  is  close  to  the 
fundamental  resonance  of  the  cavity  and  far  away  from  other  resonances,  we  can  model  an 
ac  drive  by  mapping  the  coherent  classical  signal  onto  a  quantum  operator  given  by16 

Vs(t)  =  Ho  Mi)  +  a8(tf)  ->•  V0(ase~iuJst  +  dje***)  (2.79) 

where  Vo  represents  the  amplitude  of  the  signal.  We  can  rewrite  this  as 

14(f)  =  2Vq  (Xi  cos(c ust)  +  X2  sin(o;sf)) ,  (2.80) 


where 


(2.81) 


are  identified  as  quadrature  operators.  These  two  operators  are  Hermitian  and  therefore 
observable,  and  represent  the  amplitude  of  each  of  the  two  quadratures  of  a  classical  signal 


oscillating  90°  out  of  phase  from  each  other,  just  like  the  quadratures  that  can  be  mea¬ 
sured  or  manipulated  by  a  microwave  IQ  mixer.  Signal  modulation  (on  timescales  long 
compared  with  27r/ws)  is  achieved  by  varying  X\  and  X2.  But  as  quantum  mechanical 


16We  use  the  electric  field  phase  convention  of  (B.10),  which  we  take  to  be  proportional  to  the 
voltage  operator.  Note  that  we  use  the  trajectories  of  as(t)  and  as(f)l  in  the  coherent  limit  and 
explicitly  include  the  time-dependence.  We  also  use  a  and  as  shorthands  for  as(0)  and  aj(0) 
respectively. 
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operators,  these  variables  do  not  commute.  Using  the  commutation  relations  for  a  and  a t 
we  straightforwardly  calculate 

[Xx,X2\=1-,  (2.82) 

so  that  by  the  Heisenberg  uncertainty  principle  (2.78)  we  have 

(AXf)(AXl)  >  i  (2.83) 

We  can  use  this  general  fact  about  quadratures  to  deduce  a  limit  on  the  performance  of 
linear  amplifiers  imposed  by  quantum  mechanics.  Let  us  write  the  quadrature  operators  of 
the  output  (i.e. ,  amplified)  signal  as 


Yi  =  VG~iX  i  +FU  Y2  =  s/G~2X2  +  F2 , 


(2.84) 


where  X\  and  X2  now  denote  the  input  quadratures,  G\  and  G2  denote  their  (potentially 
different)  power  gains,  and  F\  and  F2  are  operators  that  denote  quantum  fluctuations, 
which  are  necessary  for  the  output  operators  to  obey  the  canonical  quadrature  commutation 
relation  (2.82).  Specifically,  imposing  this  commutation  relation  on  the  output  and  input 
quadratures,  we  see17 

[W,y2]  =*-=>  VgiVg21-  +  \f^f2\  =  l- 

=>  [Fi,F2]  =  (l  -  y/Chy/G?)  ■  (2.85) 

The  uncertainty  principle  (2.78)  then  yields 

(AF12)(AF|)  >  ^  (l  -  V^V^)2  •  (2-86) 


Normalizing  this  to  the  gain  product  G\G2  yields  the  amplifier  uncertainty  principle  for 
added  noise,  as  described  in  1982  by  Carlton  M.  Caves  [51], 


A\A2  >  —  (  1  - 


16 


VG^VGl) 


(2.87) 


For  the  special  case  of  equal  gain  for  both  quadratures,  (y/Gi  =  VG  =  A\/G2,  A\  = 
A2  =  An/ 2  [the  factor  of  two  is  there  because  each  quadrature  contributes  to  half  of  the 
total  added  noise  when  both  quadratures  are  treated  on  equal  footing]),  this  reduces  to  the 

17We  assume  [A.;,  F3]  =  0,  meaning  that  the  input  and  internal  noise  fluctuation  modes  are 
independent. 
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fundamental  theorem  for  so-called  phase-insensitive  amplifiers, 


An  > 


1 

2 


(2.88) 


If  the  upper  sign  applies,  we  use  the  term  “phase-preserving,”  while  if  the  lower  sign  applies, 
we  use  the  term  “phase-conjugating.”  We  see  that  this  is  precisely  equation  (2.76)  in  the 
zero-temperature  limit,  which  shows  that  our  parametric  amplifier  when  operated  in  the 
phase-preserving  limit  is  quantum-limited  (we  would  obtain  the  analogous  result  for  the 
phase-conjugating  limit,  which  would  involve  considering  the  idler  frequency  —oj  as  the 
output  frequency).  In  particular,  this  means  that  increasing  the  gain  of  the  parametric 
amplifier  (within  the  range  in  which  it  functions  as  an  ideal  one)  allows  us  to  increase  the 
signal  to  noise  ratio  to  very  high  levels.  Of  course,  there  will  emerge  practical  issues  in 
achieving  ideal  parametric  amplifier  performance  in  any  real  device  -  in  particular,  in  our 
device  based  on  SQUID  loops,  we  should  be  cautious  of  magnetic  flux  noise,  which  can  lead 
to  a  fluctuating  resonance  frequency  and  therefore  a  fluctuating  gain  (we  will  mention  a  few 
precautions  taken  for  this  in  the  experimental  chapter  of  this  thesis). 

Above  we  have  in  fact  glossed  over  an  important  scenario  -  what  happens  when  the 
signal  is  exactly  resonant  with  the  pump?  In  this  case,  the  intermodulation  tone  in  (2.61) 
in  fact  has  the  same  frequency  as  the  signal  and  pump.  In  this  case,  the  gain  will  actually 
depend  on  the  phase  of  the  signal  relative  to  the  pump,  which  makes  sense  since  a  relative 
phase  can  be  defined  between  two  signals  of  the  same  frequency.  It  might  seem  strange  at 
first  that  there  is  only  a  phase  dependence  when  the  signal  is  exactly  on  resonance  with 
the  pump.  However,  real  signals  have  a  finite  (yet  potentially  very  small)  bandwidth.  In 
addition,  the  phase  dependence  can  in  fact  be  seen  at  other  frequencies  if  we  allow  for  a 
multimode  signal  and  define  so-called  generalized  quadratures  that  are  linear  combinations 
of  two  modes  symmetrically  placed  above  and  below  the  pump  frequency  (such  quadratures 
exist  essentially  because  they  consist  of  positive  and  negative  frequencies  in  the  rotating 
frame  of  the  pump.  For  a  brief  treatment  of  such  generalized  multimode  quadratures, 
see  [42].).  We  will  here  continue  to  consider  a  single  tone,  this  time  correctly  treating  the 
case  oj  =  0  on  resonance  with  the  pump,  leading  to  phase-dependent  gain  and  added  noise 
essentially  due  to  squeezing  of  the  output  mode. 

To  account  for  arbitrary  quadrature  angles,  we  generalize  (2.81)  to  the  (still  single-mode) 
quadrature  operators  for  a  mode  a 

ae~l 9  +  at  e*0 

Xq  = 


2 


(2.89) 
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so  that  X\  and  X2  in  (2.81)  correspond  to  the  special  cases  0  =  0  and  6  =  90°  respectively. 
We  wish  to  see  what  dependence  (if  any)  there  is  on  the  phase  of  the  signal  relative  to  the 
pump;  i.e.,  the  input  quadrature.  To  do  this,  let’s  decompose  the  output  (2.61)  into  two 
amplified  input  quadratures  that  are  90°  out  of  phase,  with  respective  quadrature  gains 
Aq  and  Be  [we  drop  the  frequency  arguments  here  so  that,  for  example,  M  =  M(0)  and 
cf  =  cin(0)]: 


c°ut  =  gcf  +  Mel 


inf 


=  Ae 


ere 


-id  +  c M0i$ 


+  Be 


Jne— i(0+7r/2)  _|_  ginfg i(6+ir/2) 


(2.90) 


We  can  then  solve  for  the  quadrature  gains  by  equating  the  coefficients  of  the  input  signal 
mode  and  its  adjoint,  leading  to  two  linear  equations  in  two  unknowns, 

Q  =  \e~ieAe  -  l- e~ieBe , 

M  =  \ewAe  +  '-eidBe  (2.91) 

with  the  solution 

Ae  =  Geie  +  Me~ie  (2.92) 

(since  6  is  arbitrary  there  is  no  need  to  solve  for  Be).  We  can  then  calculate  the  phase- 
dependent  gain 


Ge  =  |^6»|2 


(Geie  +  Me-ie)(G*e~i9  +  M*eie ) 

\G\2  +  \M\2  +  2 Re  \QM*e2ie 
2 Gs  —  1  +  2 \f  GS(GS  —  1)  cos  (2 [6  —  </>]), 


(2.93) 


where  we  have  used  the  definitions  (2.66)  and  (2.67),  the  identity  (2.74)  and  defined  the 
relative  phase  4>  between  the  “signal”  and  “image”  contributions  to  the  gain, 

GM*  =  \G\\M\e~2i^.  (2.94) 


We  therefore  see  that  there  exist  quadratures  90°  out  of  phase  which  maximize  and  minimize 
the  gain  (through  the  cosine  term).  In  the  high-gain  limit  Gs  3>  1,  we  obtain  the  asymptotic 
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expressions 


G*max  ^  4 Gs, 

Gmin  ~  2 Gs  —  1  —  2GS(1  —  ~  ^q2  4“  '  '  '  ) 


4  G„ 


(2.95) 


In  a  more  general  sense  (with  multimode  quadratures),  this  phase-dependence  arises  from 


(a)  Phase-Insensitive 

X, 


(b)  Phase-Sensitive 

X  min 


-  x2 


Figure  2.5  Illustration  of  phase-insensitive  and  phase-sensitive  modes  of  linear  amplifi¬ 
cation.  (a)  In  the  non-degenerate  case,  the  parametric  amplifier  acts  as  a  phase  insensitive 
amplifier,  both  amplifying  the  signal  noise  (blue  circle  mapped  to  red  circle)  and  adding 
at  least  one-half  quantum  of  noise  (brown  circle)  in  accordance  with  our  above  result,  (b) 
In  the  ideal  degenerate  case,  due  to  interference  between  Fourier  modes  symmetrically 
placed  about  the  pump  there  could  theoretically  no  added  noise  in  accordance  with  Caves’ 
theorem  (2.87).  In  this  case,  the  amplifier  must  also  amplify  one  quadrature  at  the  cost 
of  de-amplifying  the  other. 


interference  between  reflected  and  squeezed  output  tones  positively  and  negatively  detuned 
from  the  pump  [15,42].  In  the  ui  =  0  and  more  general  degenerate  multimode  case,  the  am¬ 
plifier  squeezes  the  output,  amplifying  one  quadrature  while  de-amplifying  (i.e.,  squeezing) 
the  other.  In  the  single-tone  case,  we  saw  earlier  that  the  amplifier  must  add  at  least  half 
of  a  quantum  of  noise  to  the  output  signal  via  the  image  channel.  However,  in  the  degen¬ 
erate  case,  by  choosing  the  correct  (generalized)  input  phase,  the  added  noise  can  be  made 
arbitrarily  small  in  the  desired  signal  quadrature.  To  learn  more  about  such  interesting 
squeezing  phenomena,  see  for  example  [15,40,42], 


Chapter  3 

Implementation  with 
Superconducting  Circuits:  Theory 


The  results  of  the  previous  chapter  are  very  promising,  but  they  assume  that  we  have 
a  nonlinear  quantum  oscillator  in  the  first  place.  In  this  chapter,  we  will  show  how  one 
can  construct  a  system  with  a  Kerr  Hamiltonian  out  of  a  A/4  superconducting  microwave 
coplanar  waveguide  resonator  terminated  with  a  series  of  Josephson  junction-based  super¬ 
conducting  SQUIDs.  An  implementation  with  Josephson  junctions  is  very  appealing  for 
several  reasons.  The  Josephson  junction,  which  is  based  on  superconductivity,  is  essentially 
the  only  known  nonlinear  circuit  element  that  is  practically  lossless  and  functions  at  the 
cryogenic  temperatures  from  which  the  signals  we  wish  to  pre-amplify  originate  in  the  first 
place.  In  addition,  Josephson  junctions  can  be  fabricated  on  a  chip  using  standard  clean 
room  fabrication  techniques  that  have  been  developed  over  the  past  several  decades.  This 
naturally  fits  in  with  the  circuit  QED  paradigm,  in  which  superconducting  qubits  are  sim¬ 
ilarly  fabricated.  Embedding  the  Josephson  junctions  in  a  coplanar  waveguide  resonator 
allows  us  to  control  the  coupling  (both  wanted  and  unwanted)  to  the  parametric  oscillator 
and  have  good  control  over  the  amplified  frequency  and  bandwidth.  In  the  implementation 
we  will  explore,  the  Josephson  junctions  modify  the  boundary  conditions  of  the  resonator 
and  through  their  nonlinear  inductance  give  rise  to  an  effective  Kerr  constant  for  the  sys¬ 
tem.  Capacitively  coupling  to  the  input  port  of  the  resonator  allows  one  to  input  pump 
and  signal  tones,  and  then  measure  the  amplified  signal  via  this  same  port.  We  note  that 
this  is  not  the  only  way  to  make  a  parametric  amplifier  with  Josephson  junctions  (for  ex¬ 
ample,  a  design  with  particularly  good  frequency  tunability  was  achieved  with  an  array  of 
many  SQUIDs  in  series,  forming  the  centerpin  of  a  A/2  resonator  with  separate  output  and 
input  ports  [22,42]),  but  it  is  one  that  requires  a  reasonable  pump  power  and  has  recently 
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enabled  the  demonstration  of  interesting  quantum  optics  experiments  characterizing  states 
of  microwave  radiation  at  the  single-photon  level  [15,24,31]. 

3.1  Superconductivity 

3.1.1  Background 

The  phenomenon  of  superconductivity  was  first  observed  in  1911  by  H.  Kamerlingh  O lines:, 
due  to  the  ability  to  induce  low  temperatures  thanks  to  his  1908  achievement  of  liquifying 
helium  for  the  first  time  [52],  Below  a  critical  temperature  Tc  and  applied  magnetic  field  Hc 
characteristic  to  the  material,  the  DC  resistance  of  many  metals  is  observed  to  abruptly  drop 
to  an  immeasurably  small  value  -  perpetual  supercurrents  have  been  observed  to  flow  for 
over  a  year  without  any  measurable  dissipation!  It  is  also  observed  that  the  superconducting 
metal  becomes  a  perfect  diamagnet  -  not  only  are  magnetic  field  lines  prohibited  from 
penetrating  into  the  metal  (which  might  have  been  explained  by  perfect  conductivity  alone), 
but  also  upon  the  transition  to  the  superconducting  state,  any  magnetic  flux  lines  originally 
penetrating  the  metal  are  ejected.  This  is  known  as  the  Meissner  effect. 

To  give  a  microscopic  derivation  of  the  intriguing  properties  of  the  types  of  superconduc¬ 
tors  that  are  currently  (approximately)  understood  is  a  very  formidable  task  and  certainly 
beyond  the  scope  of  this  thesis  (for  a  detailed  treatment  of  the  features  of  superconductors 
discussed  here,  see  for  example  [53]).  We  will  take  as  our  starting  point  the  results  of  the 
BCS  quantum  theory  of  superconductivity,  put  forward  by  Bardeen,  Cooper,  and  Schrieffer 
in  1957.  This  theory  describes  how  a  phonon-mediated  attractive  force  between  electrons 
allows  the  system  to  condense  into  a  nondegenerate  ground  state  separated  from  excited 
states  by  a  significant  energy  gap  2 A  [52],  unlike  in  an  ideal  Fermi  gas.  This  gap  is  of  order 
kTc  and  centered  at  the  Fermi  energy.  Put  very  broadly,  an  electron  deforms  the  atomic 
lattice,  allowing  a  second  electron  to  indirectly  interact  with  it  by  adjusting  to  the  defor¬ 
mation.  The  BCS  ground  state  involves  pairs  of  electrons  anticorrelated  in  momentum  and 
spin  forming  bosonic  Cooper  pairs  of  charge  — 2e  and  zero  spin.  The  BCS  theory  predicts 
that  it  takes  an  energy  of  2 A  to  break  a  Cooper  pair.  If  the  temperature  is  low  enough 
so  that  the  energy  gap  2A  of  this  Cooper  pair  fluid  is  larger  than  order  kT,  Cooper  pairs 
will  not  be  scattered  by  the  lattice.  For  aluminum,  Aai  ~  3.4  x  ICC4  eV  near  absolute 
zero  [52],  so  that  at  T  =  10  mK,  =Cf  is  a  few  hundred.  For  niobium,  A^b  ~  31  x  10-4 
eV  [52],  so  that  is  a  few  thousand.  Also,  compared  to  cavity  drive  photon  energies, 

which  are  on  the  order  of  hu<i=h  x  10  GHz,  ~  10  and  ~  100.  Shielding  the  sample 
and  filtering  input  lines  also  prevents  stray  radiation  from  breaking  Cooper  pairs.  It  is  thus 
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a  good  approximation  to  say  that  Cooper  pairs  will  remain  unbroken. 

If  we  wish  to  create  a  nonlinear  (i.e.,  anharmonic)  system  out  of  tiny  superconducting 
metallic  parts,  simple  LC  oscillators  alone  will  not  do.  Being  a  form  of  harmonic  oscillator, 
LC  oscillators  lead  to  equally  spaced  energy  spectra.  If  we  wish  to  create  a  qubit  out  of  tiny 
circuit  elements,  or  in  our  case  provide  an  effective  Kerr  constant,  we  will  need  a  non-linear 
element.  Luckily,  superconductivity  not  only  provides  us  with  low  dissipation  (allowing 
cQED  to  satisfy  (D3)),  but  also  provides  us  with  a  circuit  element  known  as  a  Josephson 
junction  which  unlike  other  known  non-linear  circuit  elements  is  dissipationless  and  consis¬ 
tently  retains  its  non-linear  properties  down  to  our  required  cryogenic  temperatures.  We 
will  more  specifically  see  that  the  Josephson  junction  acts  as  a  non-linear  inductor. 

The  key  feature  of  the  BCS  theory  of  superconductivity  that  we  will  use  is  the  fact  that 
since  the  Cooper  pairs  behave  like  bosons,  condensation  into  an  overlapping  common  ground 
state  is  possible.  We  can  describe  this  unique  ground  state  with  a  probability  amplitude 

ip(  r)  =  \fn( rje*<^r),  (3.1) 

where  n  =  |V’(r)|2  is  the  density  of  Cooper  pairs. 

3.1.2  Josephson  Junctions 

A  Josephson  junction  is  simply  a  thin  insulating  layer  (usually  a  layer  of  oxide)  sandwiched 
in  between  two  superconductors  (usually  two  evaporated  metal  films),  as  schematically 
shown  in  Figure  3.1(a).  The  two  superconducting  wavefunctions  decay  rapidly  in  the  in¬ 
sulator,  but  there  is  still  the  possibility  of  tunneling  as  usual  in  quantum  mechanics.  To 
start,  imagine  a  Josephson  junction  along  a  superconducting  wire,  with  the  wavefunctions 
of  each  (identical)  superconductor  given  by  ipi  =  sjnie1^1  and  V>2  =  ^/h^e*^2.  Suppose  also 
that  there  is  an  electric  potential  difference  V  applied  across  the  junction.  As  shown  in  [52], 
using  the  coupled  Schrodinger  equations  of  motion 

ifo-yy-  =  frTfa  -  eVipy  =  hTipi  +  eVip2  (3.2) 

ot  at 

(where  T  represents  the  coupling  rate  between  the  Cooper  pairs  of  different  superconduc¬ 
tors)  along  with  the  assumption  m  «  ri2  for  identical  superconductors,  one  can  arrive  at 
the  famous  Josephson  relations  predicted  by  Brian  Josephson  in  1962  for  the  current  I  and 
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(a)  Josephson  Junction 
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Figure  3.1  (a)  Schematic  of  a  Josephson  junction  consisting  of  two  superconducting 
electrodes  coupled  via  a  thin  insulating  layer.  In  our  implementation,  the  insulating  layer 
is  formed  between  two  evaporated  metal  films,  not  between  two  metal  blocks  as  drawn,  (b) 
Schematic  of  a  SQUID,  consisting  of  a  superconducting  loop  interrupted  by  two  Josephson 
junctions.  <f>  denotes  the  magnetic  flux  threading  this  loop.  Figure  modified  from  [54]. 


phase  difference1  cf>  =  —  4>i  across  the  junction: 


/  =  Ic  sin  4> 


(3.3) 


and 


<90  =  2e 
dt  h  ’ 


(3.4) 


where  Ic  is  the  critical  current ,  proportional  to  the  interaction  T.  This  describes  the  “AC 
Josephson  effect,”  which  gets  its  name  from  the  fact  that,  as  can  be  seen  by  combining 
the  two  equations,  a  constant  applied  voltage  leads  to  an  alternating  current  -  a  quite 
nonintuitive  feature  of  the  quantum  world  (similar  to  Bloch  oscillations  of  electrons  under 
the  influence  of  a  constant  electric  field  studied  in  conventional  solid-state  physics,  although 
these  are  almost  always  washed  out  by  scattering  of  electrons  due  to  imperfections  such  as 
impurities,  defects,  and  phonons  [52,55])!  We  can  also  see  how  the  Josephson  junction  acts 
as  a  non-linear  inductor.  Differentiating  the  first  Josephson  equation  with  respect  to  time 


^^The  superconducting  wavefunction  ip(r)  f°r  an  eigenstate  in  the  bulk  of  the  superconductor 
behaves  like  a  plane  wave,  with  the  phase  increasing  linearly  with  position  except  for  phase  slips 
at  material  interfaces.  In  the  lumped  circuit  element  approximation  of  a  Josephson  junction,  the 
relevant  wavelengths  are  usually  much  larger  than  the  size  of  the  circuit  elements,  so  that  we  can 
effectively  assign  a  phase  to  each  element  and  consider  only  phase  slips  at  junction  interfaces. 
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and  substituting  in  the  second  one  yields 


_  h  dl  dl 
2  elc  cos  0  dt  J  dt  ’ 


(3.5) 


indicating  an  effective  phase-dependent  (i.e.,  nonlinear )  inductance  Lj. 

We  would  like  to  have  a  Hamiltonian  describing  the  energy  associated  with  Josephson 
junctions  and  the  time  evolution  of  systems  involving  Josephson  junctions.  The  electrostatic 
part  of  the  classical  Hamiltonian  Hje\  can  be  easily  written  down.  As  in  Figure  3.1,  the 
two  superconducting  electrodes  are  capacitively  coupled  via  a  shunting  capacitance  Cj,  so 
that  if  n  =  _^e  denotes  the  number  of  excess  Cooper  pairs  on  one  of  the  electrodes,  the 
electrostatic  contribution  to  the  classical  Hamiltonian  becomes 


2„2 


H.J.el  = 


4  ezn 
~2Cj 


(3.6) 


Next  is  the  most  important  term  in  the  Hamiltonian,  the  Josephson  contribution.  The 
potential  energy  stored  in  the  junction  purely  due  to  Josephson  coupling  can  be  obtained 
by  inserting  the  Josephson  equations  into  the  usual  expression  for  the  energy  in  a  circuit 
element: 


Eii= So  lvit‘=L 

hh  r*® 

sm  (pacp  — >  — Aj  cos  0, 


2e 


(3.7) 


where  we  have  ignored  a  constant  offset  in  the  last  step  so  that 

„  ich 

Ej  =  — 

J  2e 

is  the  Josephson  energy  and  the  Josephson  contribution  to  the  Hamiltonian  is 


(3.8) 


Hj' j  =  —E j  cos  ( 


(3.9) 


3.1.3  SQUIDs 

With  a  single  Josephson  junction,  the  zero-0  energy  is  fixed  along  with  the  critical  current 
Ic.  However,  there  is  a  clever  and  simple  scheme  based  on  what  is  known  as  a  DC  SQUID 
(. superconducting  quantum  interference  device,  commonly  used  as  an  extremely  sensitive 
magnetometer)  which  can  act  as  an  effective  Josephson  junction  with  the  ability  to  tune  the 
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effective  Josephson  energy  via  applied  magnetic  flux.  If  we  have  two  Josephson  junctions 
in  a  superconducting  loop  as  in  Figure  3.1(b),  we  can  relate  the  difference  between  the 
phases  across  them  to  the  magnetic  flux  <F  through  the  loop,  due  to  the  phenomenon  of  flux 
quantization  [52],  To  arrive  at  this,  first  note  that  the  effect  of  gauge-transformations  on 
the  electromagnetic  Hamiltonian  (B.2),  as  discussed  in  Appendix  B.l,  leads  to  the  velocity 
in  semi-classical  Hamiltonian  mechanics  being 

v  =  —  (— ihV  —  qA) .  (3.10) 

m 

Plugging  in  a  wavefunction  of  the  form  (3.1)  yields  a  current  density  of 

j  =  =  —  (KVcj)  —  qA) .  (3-11) 

m 

Due  to  the  Meissner  effect,  B  vanishes  in  the  interior  so  that  the  current  density  also 
vanishes  there  by  Ampere’s  law.  Setting  (3.11)  to  zero  yields 

hXJd  =  qA.  (3.12) 


Now  imagine  traversing  a  closed  path  C  through  the  interior  of  the  superconducting  loop. 
If  we  require  the  wavefunction  to  be  single- valued,  we  need 

j)  V<^>  •  dl  =  27ts,  (3.13) 

where  s  is  an  integer.  But  by  Stokes’  theorem, 


C 


A  •  dl  =  /  /  (V  X  A)  •  der  =  /  /  B  •  der  =  4> 


D 


D 


(3.14) 


where  D  is  any  smooth  two-dimensional  surface  having  C  as  its  boundary.  Combining  this 
with  (3.12)  and  (3.13),  and  substituting  the  Cooper  pair  charge  q  =  2e,  we  obtain  the 
condition 


irhs 
4>  =  - 


e 

We  thus  see  that  flux  is  quantized  in  units  of  the  flux  quantum 


(3.15) 
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Combining  (3.13)  and  (3.15)  yields  a  relation  between  phase  and  flux, 

r  2-jt  4> 

<b  V0  •  dl  =  —  4*  =  2vr— .  (3.17) 

Jc  h  $  o 

If  0S) i  and  0S)2  are  the  phase  slips  across  junctions  1  and  2  respectively,  we  define  the 
variables 


0s  —  (0s,  1  T  0s,2)/2 

<5  =  0S,  1  -  0S, 2  (3.18) 

so  that  by  flux  quantization,  (5  =  27r<h/$o-  The  Josephson  part  of  the  classical  SQUID 
Hamiltonian  is  now  just  the  sum  of  the  two  individual  Josephson  junction  Hamiltonians: 


Us  =  — £ji  COS  0s,i  -  E.Y2  COS  0S) 2 


=  —  Ej\  cos 


5  .  „ 

+  2  )  -  En  cos 


6  5 

=  -E31  (  cos  (ps  cos  -  -  sin  <ps  Sin  - 


s  s 

—  Ej2  [  cos  0  cos  -  +  sin  6  sin  - 


_  ZTiSUm 

=  -Ej 


4> 


4>r 


cos  7r  — — —  cos  0s  +  d  sin  7 r— —  sin 


$ 


4>r 


(3.19) 


where  in  the  last  line  we  have  used  our  flux  quantization  condition,  introduced  Ejurn  = 
Eji  Ej2  and  d  =  (Eji  —  E J2) / Ef™ .  If  the  Josephson  junctions  are  made  to  be  symmetric2, 
this  effectively  acts  as  a  single  Josephson  junction  with  an  effective  Josephson  energy 


^iefF  _  ^isum 


(3.20) 


which  is  tunable  by  magnetic  flux  and  has  a  maximum  that  is  twice  as  large  as  the  single¬ 
junction  energy. 

2An  asymmetry  parameter  of  |rf|  <  1% — 10%  is  feasible  with  current  fabrication  techniques,  as  the 
energies  depend  critically  on  the  thickness  of  the  Josephson  junction  barrier  and  other  variations  in 
the  oxidation  process,  although  current  efforts  are  pushing  towards  mass  manufacturability  of  nearly 
identical  Josephson  junctions  for  quantum  integrated  circuits 
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3.2  Circuit  Quantization 

3.2.1  Linear  Part 

We  are  now  ready  to  describe  our  physical  implementation  of  a  quantum  Kerr  oscillator, 
which  is  schematically  illustrated  in  Figure  3.2.  We  first  find  the  linear  part  of  the  Hamil¬ 
tonian  of  the  isolated  resonator.  Our  starting  point  is  a  circuit  model  for  the  joint  coplanar 
waveguide  segment  +  SQUID  system  shown  in  Figure  3.2.  We  first  give  a  classical  treatment 


(a) 


(b) 


$ 


AT 


cAx\ 


- rM 

l  Ax 
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Ax 


x  =  0 


x=d 


Figure  3.2  (a)  Schematic  depiction  of  the  nonlinear  oscillator  system  coupled  to  a 
transmission  line,  (b)  Circuit  model  for  the  nonlinear  oscillator  system,  consisting  of 
a  segment  of  coplanar  waveguide  transmission  line  terminated  to  ground  by  one  or  more 
SQUIDs  corresponding  to  the  region  within  the  dashed  box  in  (a) .  The  coplanar  waveguide 
segment  is  modeled  by  an  infinitesimal  chain  of  capacitors  and  inductors  as  shown,  with  c 
and  l  being  the  capacitance  and  inductance  per  unit  length  of  the  waveguide  respectively. 


of  this  circuit  model,  and  then  quantize  the  Hamiltonian  to  obtain  a  quantum  mechanical 
model.  This  quantization  will  also  give  us  an  approximate  picture  of  the  spatial  mode  struc¬ 
ture  and  root  mean  square  voltages  of  the  transmission  line  segment.  It  may  seem  surprising 
that  we  should  even  be  able  to  notice  the  quantum-mechanical  behavior  of  a  transmission 
line  resonator  at  all,  but  at  cryogenic  temperatures,  these  tiny  superconducting  circuits 
reach  such  a  low  level  of  dissipation  that  discrete  energy  levels  become  observable.  It  is 
also  true  that  the  superconducting  gap  limits  the  degrees  of  freedom  of  the  system,  making 
them  even  simpler.  Not  only  is  the  energy  scale  of  thermal  fluctuations  small  compared  to 
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the  relevant  transition  frequencies  (kT  <C  Hujq)  and  also  to  the  superconducting  gap,  but 
the  resonator  has  a  high  enough  quality  factor  that  the  finite  widths  of  the  energy  levels 
are  less  than  the  separation  between  them.  This  also  requires  that  the  electronic  noise 
temperature  from  room  laboratory  equipment  coupled  to  the  device  be  low,  which  requires 
careful  filtering  of  signals. 

That  we  model  the  resonator  as  a  transmission  line  reflects  the  fact  that  we  are  treating 
it  as  an  essentially  one-dimensional  system.  This  is  valid  since  the  width  and  height  of 
the  resonator  are  several  orders  of  magnitude  smaller  than  its  length.  This  means  that 
many  modes  along  the  length  of  the  resonator  will  be  excited  before  any  non-fundamental 
lateral  mode  is  excited,  because  the  small  wavelengths  of  the  lateral  modes  are  associated 
with  correspondingly  much  higher  photon  energies.  Lumped  circuit  analysis  by  itself  is 
not  a  priori  accurate  for  our  resonator,  since  its  length  is  on  the  order  of  the  wavelength 
of  the  microwaves  propagating  within  it.  We  will  instead  be  justified  in  describing  the 
transmission  line  as  a  one-dimensional  series  of  infinitesimal  lumped  LC  circuits,  using  the 
ID  decomposition  of  the  transmission  line  from  Figure  3.2(b).  Our  approach  will  be  to 
consider  the  microwave  oscillations  in  each  of  the  infinitesimal  LC  circuits  (keeping  in  mind 
that  we  shouldn’t  consider  these  circuits  so  small  that  a  continuous  charge  distribution 
approximation  breaks  down).  Due  to  the  superconducting  nature  of  the  system,  we  will 
ignore  resistances  (which  would  require  LCR  circuits  rather  than  LC  circuits)  in  order 
to  simply  understand  the  ideal  quantum  behavior  of  these  degrees  of  freedom  using  the 
procedure  of  canonical  quantization. 

We  follow  the  treatment  given  by  Wallquist  et  al.  [56],  filling  in  the  gaps  of  the  derivation 
and  using  the  language  of  generalized  flux  for  the  transmission  line  to  obtain  some  intuition 
from  classical  circuits.  We  will  also  generalize  the  treatment  to  the  case  of  M  >  1  SQUIDs. 
We  will  first  calculate  the  Hamiltonian  in  the  linear  regime  (we  will  see  what  this  means  in 
terms  of  current  flowing  through  the  SQUIDs),  and  later  on  calculate  the  Kerr  constant  from 
the  leader  order  of  nonlinearity.  Let  us  divide  our  transmission  line  segment  of  length  d  into 
many  small  LC  circuits  of  length  Ax.  We  define  the  flux  variable  <f>n(i)  =  /[*  V(xn,t')  dt! 
for  each  node  between  the  nth  capacitor  and  inductor.  Then,  by  definition,  the  voltage  at  the 
nth  capacitor  node  will  be  V  (xn,  t )  =  <3?n(i).  The  voltage  drop  across  the  nth  inductor  is  then 
<hra(t)  —  4>n_i(t).  By  Faraday’s  law  of  induction,  the  current  through  the  nth  inductor  satisfies 
tE  =  Integrating  this  with  respect  to  time,  we  see  that  the  current  across  the  nth 

inductor  is  7^-($n-i(i)  —  <hn(i)).  Now,  using  the  familiar  formulas  from  freshman  physics 
[57]  for  the  energies  stored  in  a  capacitor  and  inductor,  we  can  write  down  a  Lagrangian 
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for  the  coplanar  waveguide  portion  of  system: 


-Ccpw  =  E 


cAx^2  m  _  (gnft)  ~  ^»-l(^))2 

2  nU  21  Ax 


(3.21) 


Next,  the  SQUID  array  will  have  a  contribution  to  the  system’s  Lagrangian.  We  first  con¬ 
sider  the  case  of  one  SQUID  [as  shown  in  Figure  3.2(b)]  for  simplicity,  and  later  give  a  simple 
argument  to  generalize  to  M  SQUIDs  in  series  [as  shown  in  Figure  3.2(a)]3.  There  will  be  a 
(small)  electrostatic  contribution  from  the  shunting  capacitance,  equal  to  ^T=1  2 
where  V).  is  the  voltage  drop  across  the  junction  in  the  two-junction  SQUID  loop  and  Cs 
is  the  (small)  shunting  capacitance  across  the  split  junctions  (assumed  to  consist  of  equal 
contributions  Cs/2  from  each  junction).  By  the  second  Josephson  equation  we  can  write 
this  as  2(2e)2  ('2  (<%  1  +  ^2)  =  9(2e)2  Cs<j)g,  where  </>s  is  the  effective  Josephson  junction  phase 
drop  (4>s,i  +  0s,2)/2  of  a  SQUID  loop  with  identical  junctions,  as  discussed  earlier4.  We  also 
have,  most  importantly  of  course,  the  Josephson  energy  contribution  —  Ej(&)  cos  (f>s  stored 
in  the  junctions,  which,  as  we  will  derive  in  the  next  section,  gives  a  nonlinearity  that  man¬ 
ifests  itself  as  a  Kerr  constant  for  the  joint  system.  To  make  the  dynamical  variables  of  the 
system  unitless  and  also  to  avoid  potential  confusion  between  the  <f>j  (generalized  flux  at  a 
node)  and  (magnetic  flux  threaded  through  a  SQUID  loop),  we  will  work  with  the  unitless 
superconducting  phases  cf>j  =  Finally,  using  this  convention,  there  is  one  more  term 

in  the  Lagrangian  -  the  phase  drop5  between  the  final  circuit  node  and  the  SQUID  node 
(j>s.  Adding  this  final  piece  of  bookkeeping  gives  us  the  total  system  Lagrangian6  in  terms 

3Keep  in  mind  that  the  SQUIDs  here  are  not  drawn  to  scale  -  in  reality  they  take  up  a  length 
that  is  small  compared  to  the  relevant  wavelength  of  the  resonator  so  that  they  can  be  lumped  into 
the  boundary  in  our  treatment.  We  will  also  of  course  show  why  more  than  one  SQUID  (rather  than 
simply  one  SQUID  with  a  larger  inductance)  would  be  desirable. 

4Really  we  have  (j>2  =  +  <fis,2)2  =  (j)2  1  +  2<^Sji ^Sj2  +  <j>2  2-  But  for  static  (or  adiabatically 

changing)  magnetic  flux  through  a  SQUID  loop,  we  have  <f>a  1  —  (f>Si 2  =  0,  which  when  squared  allows 
us  to  substitute  for  4>Sti<fis,2  above. 

5There  is  nothing  particularly  special  about  this  phase  drop  as  opposed  to  the  others,  but  it  is 
an  artifact  of  how  we  have  labeled  the  N  +  1  total  nodes.  Due  to  the  second  Josephson  relation, 
the  superconducting  phase  drop  across  the  SQUID  is  related  to  voltage  in  the  same  way  as  the 
generalized  phase  defined  for  the  circuit  nodes,  so  that  this  term  takes  the  same  form;  i.e. ,  the 
square  of  the  phase  difference. 

6It  may  seem  strange  to  have  factors  of  h  in  a  classical  Lagrangian,  but  this  simply  stems  from 
how  we  have  redefined  the  unitless  generalized  flux  to  have  units  of  phase,  and  the  fact  that  the 
Josephson  junction  can  be  modeled  as  a  nonlinear  inductor  only  if  it  obeys  the  Josephson  relations 
(which  involve  K)  to  begin  with. 
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of  the  unitless  dynamical  variables  (f>, 


r  _  fh  V  \cAx12u^  (<M*)  -  (j)n-l(t))2 

^ys  ~  Z^\o„)  9  9/A-r 


+  — ^Csj>2s+Emcos{<t>s)- 


h2{(j)s  -  4>n )2 
2(2e)2ZAx 


(3.22) 


We  now  make  two  approximations  limiting  us  to  a  linear  regime  for  now.  Namely,  we 
approximate  the  SQUID  as  a  harmonic  oscillator  by  Taylor  expanding  the  cosine  in  the 
Josephson  energy  to  only  second  order  and  writing  the  SQUID  part  of  the  Lagrangian  as 


/•harm 

^SQUID 


(3.23) 


This  is  valid  when  ^  <  1,  which  is  satisfied  when  the  SQUID  is  in  the  so-called  phase 
regime  (like  in  a  transmon  qubit),  Ej  S>  Eq.  Before  deriving  a  quantum  Hamiltonian 
under  these  approximations,  we  will  use  the  principle  of  stationary  action  from  classical 
mechanics  [34]  to  derive  boundary  conditions  for  the  system’s  modes  of  oscillation.  To  do 
this,  we  focus  on  the  dynamics  of  <j>s  at  the  boundary  of  the  resonator.  Replacing  the  cosine 
in  (3.22)  with  its  quadratic  approximation  and  using  the  Euler-Lagrange  equation  for  the 
dynamical  variable  (frs,  we  have 


d  \  dC 


dt  [  dcj), 


fi2Csz  ,  „  ,  2 &(<!>, -<f>N)  n 

7^**  +  Ej^°  +  2(2e)2/Ax  “  °’ 


(3.24) 


We  recognize  the  last  term  as  the  difference  quotient  for  the  first  spatial  derivative  of  <f)(x,  t) 
at  x  =  d,  so  that  taking  the  continuum  limit  and  recognizing  that 


<j>(d,t)  =  <j)s{t) 


(3.25) 


by  continuity  yields  the  boundary  condition 


h2d. 

J^<t>(d,  t )  +  Ej(<S>)d>(d,  t)  +  {2e)2Lcj\d,  t )  =  0, 


where  for  later  convenience  we  have  defined 


(3.26) 


Tcav  —  id,  CCav  —  cd. 


(3.27) 
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Next,  in  the  bulk  of  the  transmission  line  segment,  we  have  the  classical  Hamiltonian 


n/  1  d£bulk  r 

rt bulk  —  2_ ^  Vn  - -‘-bulk 


E 


94>  r 

~)2 

2eJ 


cAx  1  n\2  ,  (M)  ~  K-l{t))2 
+ - 2lKi. - 


(3.28) 

(3.29) 


We  now  introduce  the  conjugate  momentum  of  the  phase,  qn(t)  =  r'>f'?ulk  =  cAx<j)n(t), 

0<pn 

which  is  proportional  to  the  charge  on  the  nth  capacitor.  Hamilton’s  equation  of  motion 
then  gives 


2 

cAxcf>n  =  qn  =  - 


dU 

d(t>n 


/  \  2 

f  h  \  4*n+ 1  2(/>n  +  (j)n—\ 

\2  e)  lAx 


(3.30) 


Dividing  both  sides  of  this  equation  of  motion  by  cAx  and  recognizing  the  difference  quotient 
for  the  second  spatial  derivative,  we  obtain  in  the  continuum  limit  of  large  N  (i.e,  Ax  — >•  dx) 


d24>(x,t)  2  d24>(x,t) 

dt 2  V  dx2  ’ 


(3.31) 


where  v  =  1  /Vic  is  the  speed  of  propagation.  In  the  sinusoidal  steady-state,  we  can  take 
4>(x,t)  =  eiU}tg(x),  which  leads  to  an  eigenvalue  equation  for  g(x) 


d2g 

dx 2 


~—g  =  -k  g. 
vA 


(3.32) 


The  solutions  to  this  equation  depend  on  the  boundary  conditions  of  the  resonator.  The 
open  (or  capacitively  coupled)  boundary  condition  at  x  =  0  dictates  that  we  may  take  the 
current  to  vanish  at  the  ends,  so  that  the  spatial  eigenfunctions  for  (j>  have  a  vanishing 
derivative  at  the  boundary  (remember  our  expression  above  for  the  current  through  the 
nth  inductor)  and  so  are  cosines.  We  can  therefore  take  4>(x,t)  =  (j>o  sin(kft)  cos(kx). 
Substituting  this  bulk  solution  into  the  boundary  condition  (3.26)  and  dividing  through  by 
4> o  sin (kft)  then  yields 

7)2  h2kd 

-  J7T-C2  Cs(kf)2  cos  (kd)  -  '  '  sin  (kd)  +  Ej($)  cos  (kd)  =  0 

(2  e)2  (2e)-Lcav 

=>-  (kd)  tan(kd)  =  ^-LcavEj($)  -  LceLVCs(kf)2 

=>  (kd)  tan  (kd)  =  ~  ^r~( kd )2,  (3.33) 

Lj(®)  Ocav 

where  in  the  last  step  we  have  used  equation  (3.5)  for  the  effective  Josephson  inductance  Lj 
and  also  the  identity  /  =  d/\/LcavC'cav.  This  result  is  quite  interesting  and  very  important: 
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it  shows  explicitly  that  we  can  tune  the  wavelength  (and  hence  the  frequency)  of  the  allowed 
resonant  modes  of  the  system  by  tuning  the  Josephson  energy  of  the  SQUID  through  applied 
magnetic  flux  <f>.  For  realistic  resonators  and  junctions  used  in  cQED,  Cs  is  generally 
negligible  compared  to  Ccay  -  for  the  junction  areas  of  4  /mi2  we  will  use,  approximating  the 
oxide  thickness7  as  10  nm  and  its  dielectric  constant  as  1,  a  simple  parallel-plate  capacitor 
model  yields  Cs/2  ~  eo(4  /mi2)/(10  nm)  =  3  fF,  whereas  (using  parameters  we  will  calculate 
later)  our  resonator  will  have  Ccav  ~  (1.6  x  10-10  F/m)(4  mm)  =  .6  pF.  Figure  3.3  shows 


Figure  3.3  Plots  of  the  left-  and  right-hand  sides  of  equation  (3.33).  The  black  curve 
shows  tan(fcd),  while  the  green,  blue,  and  red  curves  denote  the  right-hand  side  of  the 

equation  for  values  ^-LcavEj  =  1,  8,  and  18  respectively.  The  small  ratio  Cs/Ccav  ~  .01 
is  as  in  our  crude  estimate  in  the  main  text,  and  leads  to  a  slight  bend  in  the  colored  curves. 


a  plot  of  the  left  and  right  sides  of  this  equation,  for  three  different  values  of  Ej  relative  to 
,Q  y,r — .  The  intersections  of  a  given  colored  curve  with  the  black  curve  yields  the  solution 
for  kd,  which  in  turn  determines  the  resonance  frequency  of  the  system  when  the  resonator 
length  d  is  specified8.  The  limit  Ej  — >  oo  can  be  understood  as  a  disappearing  junction, 
leading  to  a  short  of  the  centerpin  to  ground  at  x  =  d.  Along  with  the  open  boundary 
condition  at  x  =  0,  it  then  makes  sense  that  this  corresponds  to  a  A/4  resonator  according 
to  Figure  3.3.  As  Fj(4>)  is  lowered  (comparable  to  — )>  the  solution  for  k  starts  to 
drop  at  an  increasing  rate,  leading  to  a  resonance  frequency  that  tends  to  zero.  In  between 

7This  value  is  quoted  from  Sri,  as  I  have  not  measured  it  myself.  Other  junctions  in  the  literature 
for  similar  devices  also  have  thicknesses  of  a  few  nm  (see  for  example  [54,58]). 

®There  are  further  intersections  with  the  tangent  curve  for  kd  >  7t/2,  but  we’ll  only  consider  the 
low-frequency  solutions  in  the  range  0  <  kd  <  tt/2. 
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these  two  limits,  we  can  in  principle  interpolate  to  any  resonance  frequency  we  want  by 
applying  the  right  amount  of  magnetic  flux  <h.  However,  the  sensitivity  to  flux  variations 
changes  drastically  over  this  range.  For  small  Ej  it  is  both  the  case  that  the  solution  for  k 
becomes  very  sensitive  to  <f>,  and  is  also  the  case  that  Ej ( 4> )  itself  becomes  very  sensitive  to 
<h.  Also,  in  real  devices  there  is  a  finite  amount  of  asymmetry  between  the  two  Josephson 
junctions  in  a  SQUID  loop,  which  puts  a  limit  on  how  low  the  effective  energy  of  the  SQUID 
loop  can  be  tuned.  Also,  earlier  experiments  implementing  this  kind  of  tunable  resonator 
have  observed  a  decrease  in  the  internal  quality  factor  of  the  system  when  the  SQUIDs  are 
tuned  to  smaller  energies,  possibly  due  to  increased  power  dissipated  by  subgap  resistance9 
in  the  junctions  (not  just  from  inhomogeneous  broadening  due  to  flux  noise)  [59].  On  the 
other  hand,  if  Ej  is  too  large,  the  k  is  essentially  independent  of  Ej ,  with  kd  ~  ir/2.  For 
moderate  values  of  Ej  (say,  Lcav/Lj  around  3  to  12;  for  example,  the  blue  curve  in  Figure 
3.3),  there  is  in  some  sense  a  balance  between  tunability  and  stability.  In  the  experimental 
part  of  this  thesis  (the  aim  of  which  is  to  read  out  a  system  of  a  known,  fixed  frequency),  we 
will  aim  to  fabricate  SQUIDs  whose  maximum  energies  are  at  or  slightly  above  this  middle 
regime.  For  energies  in  this  regime  or  higher,  we  can  approximately  solve  the  transcendental 
equation  (3.33)  for  Cs  =  0  by  expanding  its  reciprocal  about  kd  =  n/2  [remembering  that 
cot(x)  about  7t/2  looks  like  —  tan(x)], 


cot (kd) 
kd 


\-kd  _  Lj($) 


kd 
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kv 
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=>  k  = 
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2d 


1  + 


1 

JaS 


1  + 


LjM 


(3.34) 


where  fo  =  ,  =  is  the  resonance  frequency  of  the  corresponding  bare  A/4  resonator. 

4\/ljcavC/cav 

Continuing  our  analysis  of  the  system’s  modes  of  oscillation,  we  focus  on  the  trajectory 
of  a  single  excited  mode  of  oscillation  for  fixed  k , 


4>(x,t)  =  4>k(t)  cos(kx).  (3.35) 

In  the  continuum  limit,  substituting  this  form  for  (j>  in  equation  (3.22)  and  integrating  over 
the  resonator  (using  the  harmonic  approximation  for  the  SQUID  energy  and  ignoring  Cs) 

9For  voltages  applied  across  a  junction  which  are  less  than  the  so-called  gap  voltage  (which  is 
able  to  break  Cooper  pairs),  the  thermal  energy  fc^T  can  create  a  small  density  of  quasiparticles 
(essentially  broken  Cooper  pairs)  whose  current  is  not  in  the  superconducting  state.  This  leads  to 
a  subgap  resistance  that  depends  on  the  number  of  excited  quasiparticles  [54],  As  the  Josephson 
energy  is  decreased,  the  SQUID  inductance  increases,  causing  more  of  the  current  to  flow  through 
the  subgap  channel,  which  acts  as  a  shunt  resistance. 
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yields  the  oscillator  Lagrangian10 
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(3.36) 


where  in  the  third  line  we  have  used  the  boundary  condition  (3.33)  for  i?j(4>),  valid  in 
the  Lagrangian  for  static  or  adiabatically  changing  flux.  We  then  see  that  in  the  linear 
approximation,  the  system  acts  like  an  LC  oscillator  with  an  effective  capacitance  and 
inductance 


Ck  = 


Cr. 


[1  +  sinc(2A;d)] ,  Lk  = 


2  Lc 


2  L  (kd)2  [1  +  sinc(2kd)} 

that  depend  on  k  and  therefore  on  <h.  This  is  consistent  with  (3.32),  since 
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Wfc  = 


kd 


(3.37) 


(3.38) 


s/LkCk  yfL  cav  cc; 

We  now  obtain  the  classical  Hamiltonian  of  the  oscillator  through  the  Legendre  transform, 
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(3.39) 


We  introduce  the  charge  number  variable  rik  =  \  dg£sc  =  so  that  hnk  is  the  canoni¬ 

cally  conjugate  momentum  to  (j>k,  allowing  us  to  express  the  Hamiltonian  in  the  very  simple 
form 

7~Losc  =  Ecknk  +  ELJk,  (3.40) 

where  we  have  defined  the  effective  capacitive  and  inductive  energies  of  the  oscillator 

(2e)2  „  h2 
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(3.41) 


10  When  we  factor  out  Aa;  — >  dx  from  the  sum  over  nodes,  the  second  term  in  the  summand 
becomes  a  squared  spatial  derivative  in  the  limit  Aa:  — >  0. 
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n  can  then  be  interpreted  as  analogous  to  the  excess  charge  (in  units  of  Cooper  pairs) 
on  the  effective  “capacitor”  necessary  to  give  the  oscillator’s  charging  energy.  Since  4>k 
and  hnk  obey  Hamilton’s  equations  of  motion  as  conjugate  variables  with  Poisson  bracket 
{4>k,  hnk}  =  1,  under  Dirac’s  method  of  canonical  quantization  [60],  they  can  be  promoted 
to  operators  obeying  the  canonical  commutation  relation 


[(/)k,hnk\  =  ih  =>•  [0fc,rafe]  =  i. 


(3.42) 


As  is  standard  [61]  for  the  treatment  of  quantum  harmonic  oscillators,  we  introduce  the 
analogous  dimensionless  lowering  and  raising  operators  for  our  oscillation  mode  (dropping 
the  subscript  k), 
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(3.43a) 

(3.43b) 


(3.44a) 

(3.44b) 


(3.45) 


(we  assume  k  is  fixed  and  drop  the  subscript  from  now  on).  Using  (3.38),  the  quantum 
Hamiltonian  for  our  oscillator  in  the  harmonic  approximation  then  takes  the  familiar  form 


H0 


A  +  . 


(3.46) 


3.2.2  Kerr  Nonlinearity 

In  the  previous  section,  we  derived  a  simple  harmonic  oscillator  Hamiltonian,  which  could 
have  been  done  without  any  SQUID  present  in  the  system.  In  the  harmonic  approximation 
(3.23),  the  only  role  of  the  SQUID  is  to  dictate  the  boundary  condition  of  the  system 
resonance,  and  to  provide  tunability  of  this  resonance  via  magnetic  flux.  In  a  sense,  all  we’ve 
done  so  far  is  generalize  a  linear  coplanar  waveguide  resonator.  For  making  a  parametric 
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amplifier,  we  of  course  need  a  source  of  nonlinearity.  We  will  now  show  how  an  effective 
Kerr  term  is  added  to  (3.46)  when  we  consider  the  nonlinear  part  of  the  SQUID  inductance 
(3.5). 

Doing  this  is  a  simple  extension  of  the  machinery  we  developed  in  the  previous  section. 
We  consider  the  next  term  in  the  expansion  of  the  cosine  in  the  Josephson  energy  term  of 
(3.22),  obtaining  a  correction  to  the  Lagrangian  of 


5C 


Ej(<S>) 

24 


cos4  {kd)4>k, 


(3.47) 


where  we  have  used  (3.25)  and  (3.35).  This  leads  to  the  same  change  in  the  Hamiltonian 

—  SC,  except  with  a  minus  sign.  Using  the  identification  (3.44a),  we  obtain  a  Kerr 

oq>k 

term  for  the  Hamiltonian, 
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(3.48) 


We  have  now  arrived  at  what  we  set  out  to  obtain,  a  quantum  nonlinear  oscillator  with 
a  Hamiltonian  of  the  form  (2.10),  with  a  Kerr  constant  of 


K  =  cos4  (kd)- 


(3.49) 


We  note  that  K  is  negative,  which  is  why  we  chose  a  negative  K  for  our  plots  in  the  previous 
chapter.  This  value  of  K  is  extremely  important  for  the  parametric  amplification  process, 
as  we  saw  in  the  previous  chapter.  For  moderate  to  large  values  of  Ej  (i.e.,  small  values  of 
Lj),  we  can  more  clearly  obtain  the  dependence  of  K  on  the  Josephson  energy  through  the 
first  line  of  (3.34), 
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(3.50) 


In  general,  a  higher  Josephson  critical  current  will  allow  for  the  use  of  higher  powers 
when  driving  the  resonator,  as  more  current  is  able  to  flow  through  the  junctions  in  ac¬ 
cordance  with  the  Josephson  relations.  More  specifically,  because  of  the  first  Josephson 
equation,  the  critical  current  determines  the  range  of  currents  that  are  considered  “small” 
in  that  the  phase  drop  across  the  junction  is  much  less  than  unity.  The  small-phase  limit 
was  one  of  the  important  approximations  we  made  in  our  derivations  above,  and  so  we 
must  ensure  that  it  is  satisfied  in  experiments.  From  (3.35)  and  (3.44a),  we  can  calculate  a 
“saturation”  pump  photon  number  lVgat  (not  to  be  confused  with  Writ)  for  which  we  have 
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the  highly  undesirable  equality  (f>s  =  1.  More  specifically,  we  can  calculate 
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which  we  would  like  to  keep  much  less  than  unity  because  we  would  like  to  drive  the  Kerr 
oscillator  near  its  critical  point  for  large  gains  without  worrying  about  violating  the  small 
phase  approximation.  We  could  keep  this  value  small  by  increasing  Lj ,  but  this  would  in 
fact  decrease  the  dynamic  range  of  the  amplifier  because  iVcrit  is  proportional  to  1/|K|, 
which  decreases  with  Lj.  And  since  Writ / fVsat  also  increases  with  71,  this  places  a  limit  on 
the  gain-bandwidth  product  of  the  resulting  parametric  amplifier. 

It  is  therefore  difficult  to  build  a  good  parametric  amplifier  with  a  single  SQUID.  But, 
as  hinted  earlier,  we  can  get  around  this  problem  by  using  an  array  of  SQUIDs  rather  than 
one  physical  SQUID.  First,  we  note  that  in  the  linear  regime,  inductances  add  in  series. 
Therefore,  if  we  use  a  series  array  of  M  SQUIDs  (as  depicted  in  Figure  3.2a),  each  with  a 
Josephson  energy  M- times  larger  than  what  would  have  been  used  with  a  single  SQUID, 
the  linear  Josephson  inductance  and  therefore  the  resonance  frequency  and  tunability  of  the 
resonator  [remember  equation  (3.34)]  will  remain  unchanged.  However,  the  nonlinear  part 
of  the  Josephson  inductance  will  in  fact  decrease,  essentially  because  we  are  diluting  it  with 
linear  inductance.  To  see  this,  assume  that  the  number  of  SQUIDs  is  small  enough  that 
the  spatial  region  taken  by  them  on  the  centerpin  is  still  small  compared  to  all  wavelengths 
used  to  probe  the  resonator.  When  this  is  true,  a  lumped  element  circuit  analysis  still 
holds,  and  we  can  in  addition  assume  that  the  phase  drop  across  each  junction  is  the  same 
by  symmetry11.  The  quadratic  term  of  the  harmonic  part  of  the  Lagrangian  (3.36)  then 
inherits  M  terms,  each  of  which  is  multiplied  by  M  due  to  the  assumed  M- fold  increase 
in  Ej  and  then  divided  by  M 2  due  to  the  assumption  of  equal  phase  drop  across  the  mth 
SQUID,  (j)S)m  =  4>n/M  (referring  to  Figure  3.2b),  which  reduces  to  the  same  Lagrangian 
as  for  one  SQUID.  However,  the  non-harmonic  quartic  term  (3.47)  is  proportional  to  the 
phase  drop  to  the  fourth  power,  leading  to  an  uncompensated  factor  of  1/M2  under  this 
same  analysis.  Therefore,  the  Kerr  constant  decreases  as  1/M2,  whereas  the  ratio  NCIit/Nsat 
depends  only  on  the  total  linear  Josephson  inductance.  In  summary,  using  multiple  SQUIDs 
allows  us  to  maintain  both  dynamic  range  and  tunability  in  the  small  phase  limit.  Note 
that  derivations  very  similar  to  the  ones  in  this  subsection  were  independently  given  in  [15]. 


nWe  are  implicitly  assuming  that  the  SQUIDs  are  identical,  including  identical  couplings  to 
external  magnetic  flux.  Later  we  will  see  how  we  can  experimentally  attempt  to  make  this  a  good 
approximation . 


Chapter  4 

Implementation  with 
Superconducting  Circuits: 
Experiment 


In  this  chapter,  we  describe  the  process  of  designing,  fabricating,  and  measuring  an  actual 
Josephson  parametric  amplifier  according  to  the  predictions  of  the  previous  two  chapters. 
The  most  substantial  (i.e.,  time-consuming)  part  of  this  thesis  was  the  process  of  learning, 
practicing,  and  optimizing  various  cleanroom  micro-  and  nanofabrication  techniques,  which 
would  have  taken  much  longer  were  it  not  for  the  gracious  help  of  others  (see  Acknowledg¬ 
ments)  . 

4.1  Design  and  Fabrication 

In  this  section,  we  describe  the  process  of  designing  and  physically  building  a  Josephson 
parametric  amplifier  on  a  chip.  We  will  discuss  the  logic  behind  our  chosen  design  pa¬ 
rameters  and  the  process  of  experimentally  realizing  these  parameters,  including  the  use  of 
important  tools  ranging  from  electromagnetic  simulation  to  electron  beam  lithography.  A 
large  part  of  this  was  engineering  not  only  designs  but  also  clean  room  fabrication  recipes. 
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4.1.1  Coplanar  Waveguide  Resonators 

Geometry 

The  type  of  resonator  chosen  for  our  implementation  is  a  segment  of  coplanar  waveguide 
(abbreviated  CPW,  incidentally  the  initials  of  its  1969  inventor,  C.  P.  Wen),  an  approx¬ 
imately  planar  structure  consisting  of  a  centerpin  surrounded  by  gaps  of  equal  width  on 
either  side  from  ground  planes  as  illustrated  in  Figure  4.1.  The  ground  planes  extend  far 
enough  out  relative  to  the  width  and  height  of  the  centerpin  and  gaps  that  they  can  rea¬ 
sonably  be  treated  as  semi-infinite  for  the  sake  of  resonator  parameter  calculations.  We 
first  consider  a  bare  A/4  resonator  without  any  SQUIDs  embedded  at  the  shorted  end.  We 
will  want  to  design  and  test  such  a  functioning  bare  resonator  before  we  attempt  to  embed 
SQUIDs.  A  schematic  top  and  cross-sectional  view  of  the  resonator  is  shown  in  Figure 


Figure  4.1  (a)  Bare  A/4  CPW  resonator.  Light  gray  denotes  a  Nb  film,  while  dark  gray 
denotes  the  sapphire  substrate  beneath.  Note  that  this  is  not  drawn  to  scale  -  in  reality 
lj  <  d.  Note  also  that  the  color  chosen  to  represent  the  Nb  him  is  somewhat  accurate 
with  respect  to  how  it  appears  to  the  naked  eye,  but  the  sapphire  substrate  is  not  darker 
and  is  in  fact  transparent  (pictures  of  real  devices  will  be  shown  below).  The  left  end  of 
the  CPW  is  coupled  to  an  output  line  via  a  coupling  hnger  capacitor  of  with  capacitance 
CK.  (b)  Cross-section  view  of  the  bulk  of  the  CPW.  Note  that  this  is  not  drawn  to  scale  - 
in  reality  IC/i.  er  denotes  the  dielectric  constant  (relative  permittivity)  of  the  substrate. 


4.1(a)  and  4.1(b)  respectively.  It  is  essentially  a  small  two-dimensional  version  of  a  coax 
cable,  with  the  gap  at  the  left  end  acting  like  a  mirror  that  can  reflect  or  transmit  pho¬ 
tons.  The  centerpin  height,  width,  and  gaps  ensure  that  the  resonator  can  be  treated  as  a 
one-dimensional  system,  as  was  done  in  the  previous  chapter  in  the  context  of  embedded 
SQUIDs.  The  choice  of  such  a  resonator,  which  supports  quasi-TEM  waves  (i.e. ,  essentially 
transverse  electric  and  magnetic  fields,  with  small  longitudinal  components  [16]),  is  a  good 
one  for  several  reasons.  Both  the  centerpin  and  ground  plane  can  be  easily  fabricated  on 
a  chip  with  standard  clean  room  techniques  in  a  single  layer,  without  using  any  dielectric 
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thin  films,  which  can  be  lossy.  The  CPW  has  a  characteristic  impedance  that  can  be  easily 
related  to  the  ratio  of  the  centerpin  width  to  gap.  This  means  that  we  can  connect  to  the 
CPW  from  larger  bonding  pads  while  impedance-matched  by  gradually  changing  the  width 
of  the  centerpin,  as  long  as  the  ratio  of  centerpin  width  to  gap  remains  constant  during  this 
change.  Because  of  this,  we  are  also  essentially  free  to  choose  the  width  of  the  centerpin 
so  that  it  is  both  easily  fabricated  with  photolithography  (described  later)  and  so  that  it 
is  of  a  convenient  size  for  the  embedding  of  SQUIDs.  The  use  of  niobium  on  sapphire  is 
known  to  consistently  yield  high  CPW  resonator  quality  factors  as  well.  And  of  course, 
the  CPW  allows  a  straightfoward  embedding  of  SQUIDs  for  inherited  nonlinearity,  as  seen 
in  the  previous  chapter,  due  to  the  presence  of  a  centerpin  close  to  the  ground  planes,  as 
opposed  to,  say,  a  slotline  waveguide. 


Frequency 


To  design  a  CPW,  we  need  to  understand  how  its  resonance  frequency,  characteristic 
impedance  (determined  by  capacitance  and  inductance  per  unit  length),  and  quality  factors 
are  determined  by  the  CPW  (plus  coupling  capacitor)  geometry  and  material  parameters. 
The  resonance  frequency  of  the  bare  A/4  resonator  is  simply  given  by  equation  (3.38)  [or 
equation  (3.34)]  in  the  limit  of  zero  SQUIDs  (i.e.,  zero  SQUID  inductance).  However,  un¬ 
like  in  Figure  3.2(b),  we  are  now  coupling  the  resonator  at  its  x  =  0  to  a  coax  line  via 
a  coupling  capacitor  (as  in  Figure  3.2(a)  and  Figure  4.1).  To  a  good  approximation,  this 
capacitive  loading  effectively  adds  to  Cc£cv,  yielding  a  downward-shifted  resonance  frequency 


of  approximately  [62] 


(4.1) 


As  we  will  show,  desirable  coupling  capacitances  for  our  paramp  are  small  enough  compared 
to  C/av  that  we  can  essentially  design  the  resonance  frequency  and  coupling  capacitance 
independently  (to  within  ~  1%).  If  we  wish  to  use  our  paramp  to  read  out  a  superconducting 
qubit  stored  in  a  cavity  of  frequency  6.85  GHz1,  in  order  to  reach  a  compromise  between 
stability  and  tunability  (as  discussed  in  the  context  of  Figure  3.3),  we  will  aim  for  a  bare 
quarter- wave  resonance  frequency  of  8  GHz.  This  also  conveniently  coincides  with  the 
nominal  maximum  frequency  specification  for  a  circulator  we  will  use  in  the  experimental 
setup.  To  determine  the  resonance  frequency,  we  therefore  need  to  know  the  inductance 
l  and  capacitance  c  per  unit  length  (or  the  characteristic  impedance2  Zc  =  \/Tfc)  of  a 


1This  is  the  cavity  resonance  frequency  of  a  cQED  sample  from  the  Houck  lab  that  we  will  use 
to  practice  operating  the  paramp  for  time-resolved  cQED  qubit  state  readout. 

2See  Appendix  C  for  a  brief  explanation  of  characteristic  impedance. 
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reasonable  CPW,  which  along  with  the  length  d  of  the  resonator  determines  its  resonance 
frequency. 


Characteristic  Impedance 

Defining  the  geometrical  parameters  [referring  to  the  dimensions  shown  in  Figure  4.1(b)] 


b  tanh(g) 


y=J 1^2  7>  _  gN 

^3-V1  *3.  K~  K(k)K(k'3)' 


(4.2) 


where  it'  is  the  complete  elliptic  integral  of  the  first  kind  (computable  with  the  Mathematica 
function  EllipticK  [] ),  it  can  be  shown  using  conformal  mapping  techniques  that  the 
characteristic  impedance  in  the  bulk  of  a  CPW  is  given  by  [7, 63] 


Zp/2  (  K(k)  Kjk^y1 
V^eff  \K(k')  K(kQJ 


(4.3) 


where  Zq  is  the  “impedance  of  free  space,”  Zq  =  \J go / eo  ~  1207T  D,  and 


K<1  +  e--)/2  <44> 

is  approximately  equal  to  the  average  of  the  free  space  and  substrate  dielectric  constants 
and  gives  the  speed  of  propagation  ueg  =  c/y/JDgee i  ~  c/y/eeff  for  nonmagnetic  substrates. 


Dielectric  Substrate 

To  evaluate  these  expressions,  we  need  to  know  the  dielectric  constant  of  our  substrate.  We 
chose  a  sapphire  (AI2O3)  substrate  partially  for  different  reasons.  Our  lab  (with  the  aid  of 
the  Houck  lab)  has  already  demonstrated  high-Q  A/2  resonators  based  on  niobium  (Nb)  on 
sapphire.  Our  lab  also  already  had  two  Nb-on-sapphire  wafers  but  little  remaining  Nb  on 
silicon  (Si).  Sapphire  is  also  known  to  allow  much  higher  internal  quality  factors  (i.e.,  smaller 
linear  loss  72  in  the  context  of  previous  chapters) .  The  choice  of  sapphire  leads  to  some  minor 
extra  work  in  the  fabrication  process  (it  is  harder  to  dice  into  chips  and  can  lead  to  charging 
problems  under  a  scanning  electron  microscope  [SEM]  and  during  e-beam  lithography) ,  but 
as  described  later  these  are  overcome  somewhat  easily.  The  substrate  originally  came  as  a 
53  x  53  nun2  square  wafer  of  h  =  .5  mm  thick  polished  c-plane  Sapphire  substrate  of  5N 
purity  from  CrysTec  Kristalltechnologie  (then,  t  =  200  nm  of  Nb  was  sputtered  on  top  by 
Star  Cryoelectronics  in  a  2  x  2  square  inch  region).  Sapphire  in  fact  has  an  anisotropic 
dielectric  constant  which  is  different  along  different  axes.  For  our  c-plane-oriented  crystal 
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substrate,  we  must  in  fact  substitute  an  effective  permittivity  given  by  [64] 

er,eff  =  y/ex,y  ■  ez  =  V9.3-11.5  «  10.34  (4.5) 

and  an  effective  height 

hes  =  ■  h  ss  0.90/i  (4.6) 

(the  first  of  these  modifications  is  much  more  important  than  the  second). 

We  would  like  the  centerpin  to  be  wide  enough  to  perform  accurate  photolithography 
while  thin  enough  that  superconducting  vortices  do  not  become  trapped  there3  (although, 
we  will  employ  sufficient  magnetic  shielding  so  that  the  latter  is  no  longer  a  limiting  factor 
on  our  centerpin  width).  We  also  want  the  gap  size  to  be  as  small  as  reasonable  with 
photolithography,  so  that  the  electromagnetic  field  is  confined  to  a  small  volume,  which 
should  help  avoid  certain  parasitic  modes  and  unwanted  decay  through  surrounding  lossy 
materials.  We’ve  seen  that  w  =  10  gin  width  worked  well  for  A/2  Nb  resonators  for  a 
different  experiment,  so  keeping  this  value,  along  with  a  gap  s  =  4.2  /mi  yields  Zc  ~ 
49.3  12,  matching  the  standard  50  Q  characteristic  impedance  of  microwave  equipment  (a 
close  match  is  not  important  since  the  CPW  is  somewhat  decoupled  from  the  laboratory 
transmission  line  through  a  coupling  capacitor,  although  as  we  will  shortly  see,  a  small 
unaccounted  amount  of  kinetic  inductance  should  actually  bring  Zc  closer  to  50  fi).  For  a 
short  discussion  of  what  characteristic  impedance  means,  what  impedance  matching  is,  and 
how  it  relates  to  the  velocity  of  wave  propagation  in  the  CPW,  see  (the  brief)  Appendix  C. 


Inductance 


Next,  we  need  to  know  l,  the  inductance  per  unit  length,  which  along  with  Zc  can  be  used 
to  calculate  c  and  therefore  the  A/4  resonance  frequency.  For  a  CPW,  the  expression  for  the 
familiar  magnetic  inductance  lm  per  unit  length  is,  quite  conveniently,  purely  geometrical, 
determined  (through  conformal  mapping  techniques)  by  the  equation  [7, 66] 


^  go  K(k') 
4  K(k)  ' 


(4.7) 


For  our  designed  ratio  of  w/b  ~  .543,  this  yelds  l  =  391  nH/rn.  We  can  now  indirectly 
calculate  the  capacitance  per  unit  length  via  the  relation  c  =  l  jZ‘l  ps  1.62  x  1CP10  F/m, 

3When  an  imperfect  strip  of  thin  superconducting  film  of  width  w  is  cooled  through  its  critical 
temperature,  the  Meissner  expulsion  effect  will  completely  occur  only  below  a  critical  field  Bm  ss 
4>o/ir2  [65].  A  width  of  roughly  5  /an  would  correspond  to  Earth’s  magnetic  field.  Otherwise, 
vortices  could  lead  for  example  to  flux  noise  in  the  SQUIDs. 
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justifying  the  comment  we  made  in  the  previous  chapter  about  the  relative  magnitudes  of 
the  cavity  and  SQUID  capacitances. 

However,  in  superconducting  resonators,  one  cannot  a  priori  neglect  kinetic  inductance 
associated  with  the  inertia  of  Cooper  pairs.  Classically,  kinetic  inductance  can  be  intuitively 
understood  by  comparing  the  equation  for  the  energy  stored  in  an  inductor  by  a  current  with 
the  kinetic  energy  associated  with  a  current  -  both  are  proportional  to  the  current  squared 
(the  latter  because  kinetic  energy  depends  on  velocity  squared).  Usually,  in  normal  metals 
we  don’t  have  to  worry  about  the  energy  stored  in  the  motion  of  charge  carriers,  because 
it  is  dissipated  into  the  metal  on  the  timescale  of  scattering.  However,  in  superconductors 
this  can  have  a  noticeable  effect.  It  can  be  shown  that  the  kinetic  energy  associated  with 
the  superconducting  current  in  our  centerpin  is  given  by  [7, 53] 


E  = 


1 

2 


nsmv2dV 


1  nsrn  2 

2  Her 


-Ho\2L 
2 h  L  A  ’ 


(4.8) 


where  A  =  wt,  is  the  cross-sectional  area,  ns  is  the  density  of  Cooper  pairs,  and  A l  = 
l/i J gm/nse2  is  the  so-called  London  penetration  depth,  the  distance  to  which  a  magnetic 
field  penetrates  into  the  superconductor  -  essentially  where  supercurrents  flow.  Identifying 
the  coefficient  of  |/2  with  the  kinetic  inductance  (after  conformal  mapping  taking  into 
account  the  CPW  geometry),  it  can  be  shown  that  the  kinetic  inductance  per  unit  length 
4  of  the  CPW  is  [66] 


h 


A2 

Mo -7  g(s,w,t), 
wt 


(4.9) 


where 


g(s,w,t ) 


1 

2  k2K{k)2 


2  (w  +  s) 

+  - -  In 


w  +  s 


3.1  (4.10) 


for  our  chosen  geometrical  parameters.  This  inductance  depends  somewhat  on  temperature 
and  magnetic  field  due  to  the  London  penetration  depth,  and  in  fact  this  formula  assumes 
that  t  <  Al,  but  it  can  give  us  an  order  of  magnitude  approximation.  Assuming  \l  ~  100 
nrn  in  our  Nb  film  [7],  we  estimate  lk/lm  ~  5%.  We  then  have  l  =  lm  +  Ik  ~  400  nH/m, 
pushing  Zc  to  just  over  50  ft.  Using  equation  (4.1),  we  see  that  choosing  a  resonator  length 
of  3.8  mm  will  yield  an  uncoupled  (i.e.,  CK  =  0)  resonance  frequency  of  8.1  GHz. 


Coupling  Capacitance 

The  last,  but  by  no  means  least,  design  parameter  for  the  CPW  is  the  coupling  capacitance 
and  its  associated  finger  capacitor  geometry.  The  capacitance  CK  determines  the  coupling 
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between  the  CPW  and  the  transmission  line  that  we  use  to  send  signals  into  and  out  of 
the  resonator,  as  depicted  in  Figures  3.2(a)  and  4.1(a).  More  precisely,  it  determines  the 
coupling  strength  71  (as  assumed  to  exist  in  the  context  of  quantum  optics  in  Figure  2.1) 
as  used  in  previous  chapters. 

How  can  we  extract  a  coupling  constant  from  a  capacitor?  To  do  this,  we  will  use  energy 
considerations  to  connect  the  languages  of  quantum  optics  and  electrical  engineering.  As 
originally  defined  in  the  quantum  Langevin  equation  (2.23)  with  a  single  port  7  =  71, 
(setting  input  and  output  signals  to  zero),  71  can  be  interpreted  as  the  rate  of  energy  decay 
out  of  the  resonator,  relative  to  the  energy  stored  in  the  resonator.  To  see  this,  treat  A 
as  a  coherent-limit  c-number,  so  that  \A\2  is  the  stored  energy  (normalized  to  the  energy 
of  a  single  photon).  Identifying  the  time-derivative  of  this  with  the  rate  of  energy  decay 
through  the  port  associated  with  71,  we  see  that  the  dissipation  rate  (in  angular  frequency) 
normalized  to  \A\2  is  27  (the  factor  of  2  comes  from  the  product  rule). 

We  can  perform  a  similar  analysis  for  the  CPW  circuit.  Suppose  we  have  initially  excited 
the  coupled  CPW  oscillator  and  allow  it  to  resonate  freely.  The  electric  energy  stored  in 
the  CPW  is 

Ee\eC  =  \c  [  (V2(x,t))dx,  (4.11) 

1  Jo 

where  (•)  denotes  a  time  average.  For  Ck  «  Cca,v,  the  normal  modes  of  the  CPW  resonator 
are  very  close  to  those  of  the  uncoupled  A/4  resonator,  so  that  for  the  fundamental  mode 
the  integration  takes  place  over  one  quarter  of  a  sinusoidal  envelope  which  has  its  maximum 
voltage  near  the  coupling  capacitor,  leading  to  Ec\ec  =  7)  •  \cd(V{ 0)2),  where  the  additional 
factor  of  1/2  comes  from  taking  a  spatial  average  over  a  quarter  wavelength.  As  usual,  the 
average  magnetic  energy  stored  is  equal  to  the  electric  energy,  yielding  the  total 

E=^Ccav{V(  0)2).  (4.12) 

Next,  we  calculate  the  rate  of  energy  loss  through  the  coupling  capacitor.  Assuming  that 
the  capacitor  dominates  the  CPW  impedance  (as  we  will  see,  we’ll  choose  l/uCK  to  be 
well  above  50  H),  this  is  simply  given  by  the  power  flowing  out  of  the  CPW  and  into 
the  laboratory  transmission  line  [the  red  transmission  line  in  Fig.  3.2(a)]  of  characteristic 
impedance  Zt  ~  50  H.  Because  an  infinite  (or  properly  terminated)  transmission  line  acts 
like  a  resistor  with  a  resistance  equal  to  its  characteristic  impedance  (as  in  Appendix  C), 
since  I  =  ijJoCKV(b)  is  the  associated  current  flowing  out  into  this  load  on  resonance  we  can 
write  the  dissipated  power  as 

P  =  (I2)Zt  =  (uoCK)2(V(0)2)Zt, 


(4.13) 
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so  that  we  obtain 


P 

271  =  E  = 


CrpiV 


—  —ZtZcC^ojQ, 
7 r 


(4.14) 


where  in  the  last  equality  we  have  used  Zc  =  \J Lcav/ Ccav  along  with  equation  (3.38),  which 
for  a  A/4  resonator  says  uq  =  n / (2\/ LcavCcajV) . 

We  would  like  71  to  be  small  enough  that  the  high-Q  approximation  holds  reasonably 
well,  while  at  the  same  time  large  enough  that  the  gain-bandwidth  product  (2.70)  of  our 
paramp  is  useful,  and  also  large  enough  that  the  system  can  respond  fast  enough  to  amplify 
pulsed  readout  tones  used  to  probe  cQED  systems  (usually  on  the  order  of  a  microsecond 
[or  less  with  the  aid  of  a  paramp]  in  duration  and  10  nanoseconds  in  risetime).  The  latter 
requirement  is  easily  fulfilled  for  71/27 r  >  5  MHz.  For  a  simultaneous  gain  of,  say,  25  dB 
and  a  bandwidth  of  2  MHz,  equation  (2.70)  says  that  271/(27 r)  =  30  MHz  is  necessary. 
Plugging  this  into  (4.14),  along  with  Zc  «  Zt  =  50  14  and  uo/(2ir)  =  7  GHz  yields  CK  ~ 
26.5  fF.  Plugging  this  value  into  equation  (4.1)  for  the  loaded  resonator  frequency  then 
yields  a  resonance  frequency  of  7.9  GHz.  To  obtain  such  a  capacitance,  we  must  come  up 


Figure  4.2  (a)  Screenshot  of  ANSYS  Maxwell  software  used  for  electrostatic  simulation 
of  capacitances,  (b)  3+3  finger  capacitors  design  used  for  paramps. 

with  the  necessary  physical  geometry  for  the  coupling  capacitor.  Using  the  commercial 
electromagnetic  simulation  software  ANSYS  Maxwell,  it  was  found  that  for  reasonable 
finger  capacitor  widths  and  separations,  a  3+3  finger  design  as  shown  in  Figure  4.2(b)  gave 
capacitances  in  the  20-30  fF  range.  Using  equal  finger  widths  and  separations  of  3.3  /xrn 
and  adjusting  the  lengths  of  the  fingers  (near  100  /xm)  allowed  the  capacitance  to  be  tuned 
to  the  desired  value. 

Photolithography 

The  technique  of  photolithography  was  used  to  make  the  resonators  (along  with  the  coupling 
capacitors)  by  patterning  the  pre-existing  Nb  on  our  chip.  Some  of  the  equipment  used  in 
this  process  is  shown  in  Figure  4.3.  The  first  step  in  the  process  was  to  pattern  a  “mask” 
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with  a  Heidelberg  DWL66  ultraviolet  laser  writer  at  PRISM  MNFL.  This  machine  takes 
a  design  file  and  various  parameters  as  input  and  exposes  the  desired  pattern  on  a  layer 
of  photoresist  that  comes  pre-deposited  on  top  of  a  chrome-coated  soda-lime  square  glass 
wafer.  A  4  mm  write  head  was  used,  which  can  nominally  achieve  a  minimal  feature  size 
of  ~  800  nm,  sufficient  for  our  resonators.  The  exposed  areas  of  photoresist  were  then 
dissolved  away  to  expose  the  chrome  below,  by  rinsing  in  a  known  amount  of  developer  in 
a  Laurell  EDC650Mz  Spin  Processor,  and  then  etching  the  exposed  chrome  with  a  similar 
tool.  The  remaining  resist  is  then  dissolved  away  with  a  resist-stripping  solution,  rinsed 
with  water,  and  blown  dry  with  a  nitrogen  gun.  Before  use  it  was  further  rinsed  with 
acetone  and  isopropyl  alcohol  (IPA)  and  blown  dry  with  a  nitrogen  gun. 


Figure  4.3  Photographs  of  equipment  used  for  photolithography  in  PRISM  MNFL 
cleanroom.  All  images  are  from  the  PRISM  website,  (a)  Heidelberg  DWL66  laser  mask 
writer,  (b)  Laurell  EDC650Mz  Spin  Processor,  (c)  Karl  Suss  MJB4  Mask  Aligner,  (d) 
Headway  spinners. 

The  next  step  is  to  transfer  this  pattern  onto  our  chip.  To  do  this,  the  wafer  was  first 
cleaned  by  sonicating  for  3  minutes  in  acetone  followed  by  3  minutes  in  IPA,  immediately 
followed  by  blow-drying  with  a  nitrogen  gun.  A  layer  of  S1818  photoresist  is  then  deposited 
on  the  wafer  by  vacuum-mounting  it  on  a  spinner,  depositing  enough  to  cover  a  good  portion 
of  the  chip,  and  then  spinning  the  resist  at  4500  rpm  for  60  seconds  (with  a  4500  rpm/s 
ramp-up),  which  corresponds  to  a  thickness  of  roughly  2  /mi.  The  resist  is  then  immediately 
baked  on  a  hotplate  at  115°  C  for  2  minutes.  We  note  that  all  photolithography  for  this 
thesis  was  performed  on  a  single  lxl  square  inch  quarter  of  the  2x2  square  inch  Nb-on- 
sapphire  wafer  mentioned  earlier,  including  the  paramp  resonators  with  gaps  for  SQUIDs 
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(a)  Clean  Wafer 


Nb 


Sapphire 


(b)  Spin  Resist 


Photoresist  (SI 81 8) 
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Sapphire 


(c)  Align  Mask 
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(d)  UV  Exposure 


Cr 

Resist 

Nb 

Sapphire 

(e)  Develop 
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(f)  Plasma  Etch 


(g)  After  Etch 
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(h)  Strip  Resist 
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Figure  4.4  Idealized  illustration  of  photolithography  steps.  Note:  thicknesses  are  not 
to  scale. 


(described  shortly)  and  chips  specifically  for  the  development  and  testing  of  SQUIDs  (also 
described  later).  Accordingly,  the  pattern  on  the  mask  was  a  3- by- 10  grid  of  multiple  designs 
spaced  apart  by  (Ax,  Ay)  =  (7.25  mm,  2.25  mm)  to  fill  the  wafer  (which  is  later  diced  into 
individual  chips).  This  has  the  advantage  of  diminishing  the  beading  of  thick  photoresist 
on  on  the  edges  of  the  chip,  which  could  be  more  of  a  problem  for  smaller  chips.  It  also 
of  course  has  the  advantage  of  only  needing  to  align  a  mask  and  perform  photolithography 
once.  The  photomask  is  described  in  Figure  4.5. 

The  resist-coated  wafer  is  exposed  to  UV  light  in  a  Karl  Suss  MJB4  mask  aligner,  which 
exposes  the  sample  through  the  photomask  made  earlier,  in  hard  contact  mode  (which  has 
a  nominal  best  resolution  of  1  /mi)  with  5  s  contact  time  and  5  s  exposure  time.  Before 
mounting  the  mask  in  the  aligner,  it  is  a  good  idea  to  quickly  blow  it  with  a  nitrogen  gun 
in  case  there  are  any  deposited  airborne  particles  (even  though  the  PRISM  lithography 
room  is  a  class  100  clean  room4).  The  exposed  areas  of  resist  are  then  dissolved  away 
by  dipping  the  wafer  in  MF-319  developer  for  90  seconds,  followed  by  a  quick  rinse  in  DI 
water  and  nitrogen  blow  dry.  The  wafer  was  then  observed  under  a  microscope  before  the 
final  niobium  etching  step.  It  was  noticed  that  resonator  finger  capacitors  came  out  better 
in  the  developed  resist  away  from  the  lower  corners  of  the  wafer,  with  fingers  on  chips 
near  the  corner  of  the  chip  having  a  distorted  width  and  edge  profile  (see  Figure  4.6  for 

4This  is  according  to  the  PRISM  website,  and  means  (according  to  a  now  outdated  standard) 
that  there  are  fewer  than  100  particles  with  greater  size  than  0.5  micron  per  cubic  foot  of  air. 
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Figure  4.5  (a)  Broad  overview  of  .gds  file  of  “Paramp_Mask_l”  photomask  made  for 
paramp  project,  (b)  Labels  indicating  the  contents  of  the  various  chip  designs  on  the  mask. 
The  letters  A  through  F  followed  by  the  number  of  “0”  have  3+3  finger  capacitors  (this 
column  appears  twice  because  it  is  our  target),  while  those  followed  by  the  number  “1” 
have  4+4  finger  capacitors,  for  an  approximate  2yi  «  30  MHz  and  50  MHz  respectively 
at  the  desired  operating  point  of  6.85  GHz.  “A”  denotes  resonators  with  no  gaps  at  the 
ends  for  SQUIDs  (i.e.,  bare  A/4  resonators),  “B”  denotes  a  gap  for  ~  1  SQUID,  “C”  for 
~  2  SQUIDs,  “D”  for  ~  3  SQUIDs,  “E”  for  ~  4  SQUIDs  (this  row  is  repeated  since  as 
discussed  later  it  is  our  target),  and  “F”  for  ~  5  SQUIDs.  “GO”  (repeated  six  times)  is  a 
chip  for  dose  testing  of  SQUID  e-beam  lithography  (as  described  in  the  next  section).  Each 
such  chip  contains  10  sites  for  a  4-SQUID  array  dose  test  and  10  sites  for  a  single  SQUID 
dose  test.  G1  (appearing  three  times)  is  a  chip  for  probing  the  normal  resistance  of  the 
junctions  via  a  four-probe  measurement  to  predict  the  Josephson  energies  (as  described  in 
the  next  section) .  (c)  Picture  of  the  photomask  (inverted  when  viewed  from  chrome  side) . 


photolithography  result  images).  This  is  possibly  due  to  the  visible  beading  of  resist  near 
the  edges  (especially  the  corners)  of  the  wafer,  which  may  have  led  to  uneven  contact  with 
the  mask  during  exposure. 

The  niobium  etch  was  then  performed  with  an  STS  ICP  Metal  Etcher,  which  uses  a 
beam  of  plasma  to  etch  metal.  This  has  the  advantage  of  being  directional  (anisotropic)  to 
produce  nearly  vertical  Nb  sidewalls.  The  process  is  also  largely  automated.  For  etching 
through  200  nm  of  Nb,  the  predefined  recipe  “Houck  Nb”  was  used.  After  stripping  the 
remaining  resist,  the  wafer  was  covered  in  a  protective  layer  of  thick  photoresist  (AZ  4330) 
and  then  diced  into  2x7  mm2  chips  with  an  ADT  proVectus  7100  Wafer  Saw.  As  sapphire 
is  a  relatively  hard  material  to  cut,  a  recipe  with  several  passes  at  different  cut  depths  had 
to  be  used. 
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Figure  4.6  Optical  images  of  resonator  photolithography  results,  (a)  Broad  view  of 
left  side  of  A/4  test  resonator  chip.  The  crosses  and  square  in  which  Nb  is  etched  form 
alignment  marks  for  e-beam  lithography  (as  described  later)  on  all  chips,  (b)  Finger 
capacitors  (chrome)  on  photomask,  (c)  Developed  resist  for  bare  A/4  resonator  test  chip 
near  lower-left  corner  of  wafer.  The  resulting  pattern  is  imperfect  (the  capacitor  in  the 
chip  right  at  the  corner  is  worse  and  was  not  used),  (d)  Developed  resist  for  paramp 
chip  near  center  of  wafer.  The  resulting  pattern  is  less  imperfect,  (e)  The  resulting  Nb 
capacitors  after  etching  and  resist  strip  of  (c).  (f)  The  corresponding  capacitors  for  (d). 


4.1.2  Fabrication  of  SQUIDs 

Because  the  feature  sizes  of  SQUIDs  are  generally  on  the  order  of  100  nm  or  smaller, 
photolithography  is  no  longer  sufficient  -  one  can  instead  use  electron  beam  lithography,  in 
which  a  focused  beam  of  electrons  is  used  to  beat  the  diffraction  limit  of  light5.  In  several 
ways  that  will  become  apparent  (especially  with  the  tools  used  in  this  thesis  work),  e-beam 
lithography  is  much  more  difficult  (and  more  time  consuming)  than  photolithography. 

5Electrons,  unlike  photons,  are  massive,  and  therefore  have  a  different  dispersion  relation  between 
their  energy  and  momentum  (and  therefore  their  deBroglie  wavelength,  which  characterizes  the 
length  scales  at  which  well-defined  features  can  be  exposed  via  bombardment  by  such  particles). 
Namely,  non-relativistic  massive  particles  have  kinetic  energy  E  =  H2k2 /2m,  whereas  photons  have 
E  =  Hck ,  where  k  oc  1/A  is  the  wavenumber  characterizing  its  momentum.  This  means  that  for 
an  energy  of,  for  example,  10  keV,  electrons  can  obtain  roughly  ten  times  a  smaller  wavelength 
than  a  photon  of  the  same  energy.  This  is  also  one  of  the  crucial  operating  principles  of  electron 
microscopes. 
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Electron  Beam  Lithography 

In  a  nutshell,  electron  beam  lithography  uses  a  different  kind  of  resist  (known  as  “e-beam 
resist”)  that  is  sensitive  to  bombardment  by  electrons  of  certain  energies,  rather  than  to  UV 
light.  The  Raith  eLiNE  e-beam  lithography  tool  used  in  this  thesis  is  essentially  a  scanning 
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Figure  4.7  Rough  illustration  of  the  Dolan  bridge  technique  used  for  paranrp  SQUID 
fabrication,  leaving  out  A1  charge  layer  step.  Figure  based  on  a  similar  one  in  [58].  (a) 
Spin  bottom  layer  of  e-beam  resist,  (b)  Spin  thin  top  layer  of  less  sensitive  resist,  (c) 
Expose  pattern  to  controlled  beam  of  electrons.  Faded  yellow  arrows  in  the  central  region 
denote  a  lower  “undercut”  dose  intended  to  only  significantly  expose  the  bottom  layer  of 
resist,  (d)  Develop  exposed  resist  in  MIBK.  (e)  Evaporate  first  layer  of  aluminum,  at  angle 
9\.  (f)  Expose  surface  to  oxygen  to  create  a  thin  layer  of  aluminum  oxide,  (g)  Evaporate 
second  layer  of  aluminum,  at  angle  02-  (h)  Strip  unwanted  resist  and  lift  off  unwanted 
aluminum  on  top  of  resist. 


electron  microscope  that  has  been  modified  and  hooked  up  to  control  electronics  to  allow 
the  beam  to  scan  over  a  chip  to  expose  user-defined  patterns  with  a  user-defined  charge 
dose.  The  most  common  method  of  making  Josephson  junctions  with  e-beam  lithography  is 
known  as  the  Dolan  bridge  technique  [67],  a  bilayer  process  that  is  illustrated  in  Figure  4.7. 
In  this  process,  a  suspended  bridge  of  thin  resist  blocks  evaporated  aluminum  from  being 
deposited  on  the  substrate  in  its  shadow.  By  evaporating  metal  at  two  different  angles  and 
exposing  the  metal  to  oxygen  between  these  two  evaporations,  one  can  create  Josephson 
junctions  out  of  aluminum  and  aluminum  oxide.  The  basic  procedure  is  as  follows.  First, 
the  photoresist  protecting  the  2x7  mm2  diced  chip  is  stripped  off  and  the  chip  is  sonicated 
for  8  minutes  in  acetone  and  then  8  minutes  in  IPA.  Next,  a  layer  of  copolymer  resist 
MMA  EL13  (methyl  methacrylate,  ethyl-lactate  13%)  from  MicroChem  is  spun  at  5000 
rpm  (with  a  1000  rpm/s  ramp-up)  for  60  s  to  a  thickness  of  about  600  nrn.  The  MMA 
layer  is  immediately  baked  on  a  hotplate  at  175  °C  for  two  minutes.  The  chip  is  then 
allowed  to  cool  on  the  spinner’s  vacuum  chuck  for  1  to  2  minutes  so  that  there  are  no 
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problems  depositing  the  second  layer  of  thin  resist.  This  second  layer  is  950PMMA  A3  [i.e. , 
poly  (methyl  methacrylate),  molecular  weght  950k,  3%  in  anisole]  from  MicroChem  spun  at 
4000  rpm  (with  a  1000  rprn/s  ramp-up)  for  60  s  to  a  thickness  of  ~  120  nrn.  The  chip  is 
then  baked  for  30  min  at  175  °C,  and  then  checked  under  an  optical  microscope  to  ensure 
that  no  damage  to  the  surface  has  accidentally  occurred. 

When  using  a  sapphire  substrate,  there  is  one  more  step  between  depositing  resist  and 
writing  a  pattern.  Because  sapphire  is  a  good  insulator,  it  is  susceptible  to  charging  under 
an  electron  beam.  This  is  detrimental  to  writing  precise  patterns  with  such  a  beam,  which 
can  move  around  under  the  influence  of  such  charge.  The  niobium  on  top  is  also  a  relatively 
poor  conductor  and  will  in  addition  have  poor  electrical  contact  with  the  metal  clip  holding 
the  chip  in  place  during  e-beam  lithography,  so  that  there  is  no  sufficient  conduction  path 
off  of  the  chip  for  the  excess  charge.  To  circumvent  this,  a  thin  (12  nm)  anti-charging 
layer  of  aluminum  is  evaporated  onto  the  surface  of  the  chip  before  writing  the  pattern. 
Both  thermal  and  e-beam  evaporation  (depending  on  the  availability  of  the  evaporators) 
were  used  successfully  for  this.  As  a  minor  note,  it  was  found  that  for  similar  evaporator 
vacuum  pressures,  ebearn  evaporation  yielded  a  less  granular  aluminum  layer  than  thermal 
evaporation  did  and  also  a  lesser  rate  of  cracking  under  the  Raith’s  30  kV,  20  nm  spot  size 
electron  beam  [see  Figure  4.8(c)  and  4.8(d)],  although  the  reason  for  this  difference  is  not 
known,  as  the  two  methods  were  performed  in  different  evaporators  and  only  a  few  times 
each.  Figure  4.8  illustrates  a  few  parts  of  the  Raith  process.  The  chip  is  first  mounted  on  a 
sample  holder  with  a  metallic  clip  gently  clamping  the  very  edge  of  the  chip,  and  loaded  into 
the  main  vacuum  chamber  via  a  load-lock  system.  For  SQUID  fabrication,  the  beam  was 
set  to  use  a  20  /jrn  aperture,  30  kV  energy  and  11  mm  working  distance.  The  SEM  is  then 
focused  on  ~  50  nm  gold  particles  deposited  on  the  corner  of  the  chip  in  order  to  optimize 
stigmation  and  aperture  alignment.  Then  various  steps  are  taken  to  align  the  system’s 
coordinates  with  the  chip  axes,  and  define  the  plane  of  the  chip  along  with  a  coordinate 
system  using  a  three-point  focusing  method.  Such  three-point  alignments  were  aided  by  the 
presence  of  square-  and  cross-shaped  alignment  marks  from  the  previous  photolithography 
step.  One  must  be  careful  not  to  accidentally  expose  critical  regions  of  the  chip.  Since 
the  aluminum  anti-charging  layer  tended  to  crack  under  a  concentrated  beam,  the  edges 
of  the  cracked  aluminum  were  used  to  focus  quite  nicely  on  the  surface  for  three-point 
alignment  (the  alternative,  which  was  found  to  be  quite  difficult,  was  to  “burn”  a  spot 
on  the  resist  and  focus  on  that).  For  thermally  evaporated  aluminum  charging  layers,  it 
was  also  possible  to  focus  on  the  relatively  granular  surface.  After  performing  a  write  field 
alignment  procedure  to  calibrate  the  coordinates  and  beam  deflection,  the  beam  current 
was  measured  via  a  Faraday  cup  on  the  sample  holder  in  order  to  calibrate  the  exposure 
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Figure  4.8  (a)  Picture  of  the  Raith  eLiNE  e-beam  lithography  system,  (b)  Picture 
of  the  inside  of  the  vacuum  chamber  in  which  the  chip  resides  during  the  e-beam  expo¬ 
sure.  (c)  Cracked  aluminum  on  one  of  the  small  alignment  marks,  conveniently  used  for 
focusing  during  three-point  alignment,  (d)  High-magnification  of  cracked  aluminum  edge, 
for  two  different  samples  (left:  thermally  evaporated  aluminum,  right:  e-beam-evaporated 
aluminum) . 


“dose.”  Then  the  desired  patterns  were  written  from  a  user-defined  .gds  file  and  user-defined 
dose.  At  a  given  energy,  the  dose  refers  to  the  amount  of  charge  deposited  per  area  within 
the  region  in  question.  To  create  PMMA  bridges,  the  region  around  the  bridges  needs  to 
be  fully  exposed  without  significantly  exposing  the  PMMA  bridge  region  [green  and  blue 
boxes  respectively  of  the  SQUID-array  design  shown  in  Figure  4.9(b)].  However,  it  is  also 
necessary  for  the  MMA  under  the  PMMA  bridge  to  be  exposed,  without  the  PMMA  on 
top  itself  being  exposed  above  its  threshold.  Since  the  MMA  has  a  relatively  low  threshold 
dose  of  around  60  gC/cm2  for  a  30  kV  beam,  an  “undercut”  box  of  slightly  over  this  dose 
was  used  at  the  bridges.  Several  rounds  of  dose  tests  were  needed  to  figure  out  the  optimal 
doses  for  the  various  exposure  regions.  An  idea  of  what  this  looks  like  is  given  by  Figure 
4.9  (c)  through  (e),  which  show  some  of  the  results  of  the  first  (most  coarse)  dose  test.  This 
and  all  future  dose  tests  used  a  constant  dose  of  70  gC/cm2  for  the  undercut  regions.  The 
first  dose  test  gave  the  entire  main  SQUID  exposure  a  single  base  dose  to  be  proximity- 
corrected.  Proximity  correction  is  a  built-in  feature  of  the  Raith  software  that  attempts  to 
correct  for  the  proximity  effect ,  in  which  primary  beam  electrons  accelerated  into  a  region 
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Figure  4.9  Results  of  first  (coarse)  dose  test,  (a)  Partial  image  of  one  of  the  dose  test 
chips  on  which  dose  tests  of  SQUID  fabrication  were  performed,  (b)  Proximity-corrected 
four-SQUID  design  displayed  in  LayoutEditor.  Smaller  and  denser  green  boxes  correspond 
to  higher  doses  in  the  main  exposure.  Blue  boxes  are  low-  and  constant-dose  undercut 
regions,  (c)  Underdeveloped  resist  (top)  and  resulting  SQUID  array  after  liftoff.  Base 
dose  used  in  proximity-corrected  exposure  indicated  on  right,  (d)  Best  dose  in  first  dose 
test  for  a  proximity-corrected  main  exposure,  (e)  Result  of  overexposure  (high  dose) .  The 
PMMA  bridges  have  visibly  collapsed. 


of  resist  can  scatter  and  expose  nearby  regions  of  resist.  As  can  be  seen  in  Figure  4.9,  with 
this  dose  strategy,  by  the  time  the  base  dose  becomes  high  enough  to  sufficiently  expose  the 
regions  directly  surrounding  the  SQUIDs  and  also  the  Nb-Al  overlap  region,  the  PMMA 
bridges  are  damaged  due  to  the  nearby  exposure.  In  a  second  round  of  dose  tests,  the 
Nb-Al  overlap  region  was  given  its  own  large  and  constant  dose  of  390  gC/cm2,  which  was 
sufficient  to  circumvent  any  significant  specks  of  undeveloped  resist.  In  this  second  dose 
test,  the  SQUID  region  was  still  proximity-corrected,  and  had  fewer  undeveloped  spots  of 
resist  near  the  SQUIDs.  The  first  two  dose  tests  were  performed  with  PMMA  bridges  of 
length  3.5  gm  and  ~  350  nrn  width.  These  10:1  aspect  ratio  bridges  would  occasionally 
collapse  with  a  dose  that  previously  produced  stable  bridges.  Reducing  the  bridge  length 
a  little  bit  made  them  more  robust  and  also  significantly  reduced  the  amount  of  specks  of 
undeveloped  resist  near  the  SQUID  loops.  Note  that  all  dose  tests  were  performed  with 
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the  exposure  done  in  a  region  that  emulates  the  actual  environment  used  in  a  real  paramp 
device  (i.e.,  a  centerpin  and  surrounding  Nb  ground  plane),  since  the  optimal  dose  could  be 
sensitive  to  a  significantly  different  immediate  surrounding  environment. 

After  exposure,  the  resist  must  be  developed.  The  A1  anti-charging  layer  is  first  removed 
by  dipping  the  chip  in  AZ  300  MIF  for  70  seconds,  followed  by  a  short  rinsing  dip  in  DI 
water  and  nitrogen  blow-dry.  The  resist  is  then  developed  by  gently  dipping  the  chip  in  1:3 
MIBK:IPA  for  50  seconds,  followed  by  a  10  second  IPA  dip  and  nitrogen  blow-dry.  Care  was 
taken  not  to  collapse  the  bridges  by  swirling  the  chip  or  using  too  much  pressure  on  the  N2 
gun.  Again,  sample  optical  images  of  exposed  resist  are  shown  in  Figures  4.9  (c)  through  (e). 
The  next  important  step  in  SQUID  fabrication  is  the  double-angle  evaporation  illustrated 
in  Figure  4.7  (e)  through  (g).  For  this,  the  sample  is  loaded  into  a  metal  evaporator  on  a 
stage  that  is  capable  of  tilting  at  user-specified  angles.  Before  the  evaporation,  a  moderately 
strong  plasma  etch  is  performed  to  both  clean  the  surface  (for  example  organic  residue  and 
oxide)  and  attempt  to  shave  the  edge  of  the  niobium  step  over  which  the  evaporated  A1 
must  transition  from  the  substrate  to  the  niobium  in  a  continuous  manner.  This  extra  etch 
was  found  not  to  damage  the  PMMA  bridges,  but  it  is  not  clear  whether  or  not  it  smoothed 
the  niobium  edges  [an  SEM  image  of  one  of  the  resulting  overlaps,  which  were  found  to  be 
satisfactory  for  evaporation  angles  of  ±45°,  is  shown  in  Figure  4.10(c)].  The  Al-Nb  overlap 
was  chosen  to  be  large  (~  600  /.mi2)  so  that  at  paramp  operating  frequencies,  the  capacitance 
essentially  acts  as  an  ac  short  even  if  there  is  little  dc  contact.  Large  evaporation  angles 
were  used  for  the  sake  of  Al-Nb  contact  and  also  to  obtain  a  large  SQUID  area  for  desired 
Josephson  energy,  as  will  be  discussed  shortly.  Figure  4.10  shows  environmental  SEM  (in 
environmental  mode,  to  avoided  charging)  images  of  the  result  of  one  of  the  dose  tests. 
After  the  double-angle  evaporation,  the  entire  chip  is  covered  with  aluminum.  We  can  “lift 
off”  the  unwanted  aluminum  (which  is  sitting  on  top  of  a  bilayer  of  e-beam  resist)  by  letting 
the  chip  sit  in  acetone.  The  junction  aluminum  deposited  directly  on  the  sapphire  substrate 
and  niobium  will  remain  attached,  whereas  when  the  resist  is  stripped  by  the  acetone  the 
aluminum  on  top  of  it  will  be  lifted  off.  Liftoff  was  performed  for  3  hours  at  70  °C.  Care 
must  be  taken  to  make  sure  the  acetone  does  not  evaporate  and  expose  the  chip  to  air,  as 
liftoff  can  be  difficult  or  impossible  if  that  happens.  After  3  hours,  a  pipette  is  used  to  squirt 
the  chip  50-100  times  so  that  the  aluminum  visibly  strips  off,  ideally  in  one  piece  (generally, 
it  was  found  that  smaller  pieces  would  stick  to  the  edges  longer  than  the  bulk).  The  chip 
is  then  sonicated  for  5  seconds,  dipped  in  IPA  for  5  minutes,  sonicated  for  2  seconds,  and 
then  blown  dry  with  a  nitrogen  gun.  After  optical  inspection,  the  sample  is  now  ready  for 
resistance  tests  or  packaging,  depending  on  the  purpose  of  the  chip. 
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Figure  4.10  Environmental  SEM  images  of  a  SQUID  array  from  the  second  round  of 
dose  tests,  (a)  Overview  of  SQUID  array  and  overlap  with  Nb  centerpin  and  ground  plane 
(the  e-beam  alignment  was  a  little  off  on  this  run) .  Small  artifacts  presumably  from  specks 
of  undeveloped  resist  are  visible,  (b)  Side  view  of  one  of  the  SQUID  loops.  The  thin  film 
seen  is  possibly  the  “veil  of  death”  (seen  and  colloquially  named  so  by  other  groups  [7]), 
potentially  due  to  aluminum  scattered  by  residual  outgassing  of  resist.  The  rough  edges 
of  the  SQUID  loops  are  possibly  due  to  imperfect  liftoff  or  the  brief  sonication  used  after 
liftoff,  (c)  Side  view  of  Al-Nb  overlap,  (d)  Side  view  of  one  of  the  Josephson  junctions. 


SQUID  Characterization 

With  optimal  doses  in  hand,  we  can  start  to  think  about  tuning  the  SQUID  parameters  to 
obtain  the  desired  Josephson  energies.  In  accordance  with  Figure  3.3  and  equation  (3.34), 
an  effective  SQUID  energy  of  1000  GHz  (1  THz)  would  allow  us  to  reach  a  compromise 
between  tunability  and  stability,  reaching  the  desired  operating  frequency  of  6.85  GHz  at 
around  three-quarters  of  the  maximum  SQUID  energy.  To  have  a  reasonable  critical  power 
in  accordance  with  equations  (2.47)  and  our  discussion  in  section  3.2.2,  an  array  with  four 
SQUIDs  (therefore  each  of  energy  4  THz)  was  found  to  be  sufficient  to  obtain  a  reasonably 
high  critical  input  power  of  —100  dBm  at  an  operating  frequency  of  6.85  GHz. 
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How  can  we  predict  the  SQUID  Josephson  energies?  Remarkably,  it  can  be  shown  [53,68] 
that  the  critical  current  of  a  Josephson  junction  obeys  the  Ambegaokar-Baratoff  relation 


Ic 


A  (T) 
2  eRn 


(4.15) 


where  A (T)  is  the  (temperature-dependent)  superconducting  gap  and  Rn  is  the  normal 
state  resistance,  which  is  measurable  at  room  temperature!  This  resistance  is  very  sensitive 
to  the  junction  (i.e. ,  oxide)  thickness,  and  should  also  scale  inversely  with  the  area  of  the 
junction.  To  get  the  desired  resistance  (73  D  per  junction  in  a  four-SQUID  series  array),  we 
can  therefore  fabricate  various  SQUID  sites  while  varying  junction  geometry  and/or  oxida¬ 
tion  time.  Based  on  the  Houck  lab’s  standard  recipe  for  much  smaller  junctions  (resistances 
of  kiloohms),  a  naive  extrapolation  by  junction  area  would  tell  us  that  a  junction  area  of 
roughly  4  p2  is  desired,  which  is  why  relatively  long  bridges  and  large  evaporation  angles 
were  chosen.  This  neglects  edge  effects  in  the  geometry,  but  is  a  reasonable  starting  point. 
This  estimate  turned  out  to  perform  better  than  expected.  Two  series  of  resistance  test 
chips  were  iterated,  in  which  a  four-probe  lock-in  amplifier  measurement  (at  80  Hz)  was 
used  to  measure  the  resistance  of  SQUID  arrays  of  varying  bridge  lengths  (and  hence  junc¬ 
tion  areas).  A  four-probe  measurement  was  used  because  the  desired  SQUID  resistances 
are  relatively  low  (less  than  100  fl),  and  we  want  to  isolate  only  their  resistance  by  measur¬ 
ing  the  voltage  drop  across  them  with  a  known  amount  of  current  flowing  through  them. 
Furthermore,  niobium  is  known  to  have  a  very  high  contact  resistance  when  probed  with 
the  small  electrodes  of  our  Lakeshore  probe  station,  so  aluminum  electrodes  were  fabricated 
over  the  initially  deposited  niobium  electrodes  [see  Figure  4.11(a)].  For  the  electrode  e-beam 
lithography  step,  an  aperture  of  120  /mi  was  used  for  the  sake  of  time.  Care  was  taken  to 
remain  electrically  grounded  (through  a  wrist-strap  and  BNC  tees  that  could  be  shorted  at 
one  end  when  contacting  the  probes)  so  that  the  SQUIDs  do  not  become  damaged  from  a 
large  electrical  discharge6.  Figure  4.11(d)  shows  the  measured  four-probe  resistances.  The 
slope  of  the  data  with  respect  to  bridge  length  makes  sense,  with  the  resistance  roughly 
halving  as  the  junction  area  doubles.  These  resistances  are  in  the  desired  range  (~  150  D 
for  a  series  of  four  SQUIDs). 


6Unless  one  is  in  the  mood  for  fried  calamari  and  bad  puns. 
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PMMA  Bridge  Length  (pm) 


Figure  4.11  (a)  Optical  image  of  a  four-probe  resistance  test  site,  (b)  Closeup  of  elec¬ 
trode  overlaps  near  SQUID  array,  (c)  Picture  (of  a  picture)  of  a  four-probe  measurement 
being  performed  with  the  Lakeshore  probe  station.  The  probes  are  manually  lowered  onto 
the  bond  pads  for  good  electrical  contact,  (d)  Results  of  resistance  tests.  Note  that  while 
the  dependence  on  bridge  length  makes  sense,  the  dependence  on  oxidation  time  was  op¬ 
posite  from  what  would  be  expected.  However,  the  variation  was  relatively  small,  and 
it  seems  that  both  oxidation  times  already  exceed  the  saturation  point  of  the  oxidation 
process,  so  that  the  variation  is  due  completely  to  unknown  (or  uncontrollable)  variations 
between  the  two  fabrications  (it  was  checked  that  junction  aging  was  not  the  reason.) 

4.2  Amplifier  Implementation  and  Characteriza¬ 
tion 

After  dialing  in  on  a  coupling  capacitor,  e-beam  lithography  dose,  and  SQUID-array  resis¬ 
tance,  we  are  ready  to  try  and  make  a  Josephson  parametric  amplifier.  After  the  junctions 
are  fabricated,  the  chip  must  be  packaged  in  a  way  that  allows  us  to  connect  the  resonator 
input  port  to  our  laboratory  equipment.  This  is  done  by  mounting  the  chip  in  a  notch 
of  matching  size  on  a  custom  printed  circuit  board  (PCB)  [see  Figure  4.13(a)]  with  silver 
epoxy'.  The  launch  pad  leading  to  the  centerpin  of  the  resonator  must  be  connected  to 

7For  the  epoxy  to  set,  it  was  necessary  to  heat  the  PCB  (along  with  the  sample)  to  120  °C  in  a 
vacuum  oven.  To  make  sure  this  wouldn’t  change  the  junction  resistances,  this  was  tested  on  one  of 
the  resistance  test  chips,  and  it  was  found  that  the  process  increased  the  resistances  by  only  up  to 
around  10  percent. 
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the  centerpin  of  our  coax,  which  is  done  via  three  wire  bonds.  A  dense  array  of  wire  bonds 
was  also  used  to  connect  the  ground  plane  of  the  resonator  to  the  ground  of  the  PCB,  in 
order  to  avoid  unwanted  resonances.  In  general,  shorter  bonds  are  desirable  (though  harder 
to  achieve)  because  they  add  less  inductance  to  the  circuit  and  have  better  high-frequency 
properties.  For  the  paramp  measured  in  this  thesis,  extra  wire  bonds  across  the  middle 
of  the  chip  connecting  the  two  sides  of  the  ground  plane  across  the  centerpin  were  used 
to  further  reduce  the  chance  of  spurious  modes.  A  picture  of  the  wire  bonding  process  is 
shown  in  Figure  4.12(a)  and  images  of  the  resulting  bonded  chip  are  shown  in  Figure  4.13(a) 
through  (c). 


Figure  4.12  (a)  Picture  of  the  wirebonder  in  use.  (b)  From  left  to  right:  octobox,  shim, 
and  lid.  (c)  Illustration  of  how  sample  is  mounted  via  surface  mount  SMP  (subminiature 
coax)  connectors  and  SMP  bullets  in  the  fridge  via  the  octobox. 

The  PCB  and  sample  are  then  covered  with  a  shim  [similar  to  the  one  shown  in  Figure 
4.12(b)]  and  mounted  to  the  cover  piece  of  an  “octobox,”  which  is  capable  of  providing  up 
to  8  SMP  connections  (hence  the  name).  This  reduces  the  volume  that  could  contribute  to 
unwanted  parasitic  modes.  The  cover  piece  is  then  mounted  on  the  octobox  in  a  manner 
that  allows  for  the  sample  to  connect  to  metal  coax  cables  in  the  fridge  on  a  rod  that 
will  hold  it  in  place  in  a  magnetic  shield  [see  Figures  4.12(c)  and  4.13(g)].  Bright  field 
optical  micrographs  of  the  paramp  device  measured  in  this  thesis  are  shown  in  Figure 
4.13  (d)  through  (f).  For  the  paramp,  a  PMMA  bridge  length  of  2  /jrn  was  used.  The 
reader  may  have  noticed  that  the  region  of  niobium  ground  plane  on  the  right  hand  side 
of  the  resonator  has  been  removed,  except  for  a  square  region  right  near  the  SQUID  array 
termination.  This  was  an  attempt  to  minimize  inhomogeneous  SQUID-flux  coupling  that 
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Figure  4.13  (a)  PCB  with  connectors  and  two  paramp  samples  mounted.  The  PCB 
copper  directly  surrounding  the  chip  is  shinier  because  it  was  polished  with  a  fiberglass  pen 
(and  subsequently  sonicated  in  acetone  and  IPA)  in  order  to  ensure  good  wire  bonds,  (b) 
Closeup  of  connection  from  SMP  to  chip,  (c)  Closeup  of  chip  with  individual  wire  bonds 
visible,  (d)  Bright  field  optical  micrograph  of  paramp  finger  capacitors,  (e)  Bright  field 
optical  micrograph  of  paramp  SQUID  termination  region,  (f)  Close-up  bright  field  optical 
micrograph  of  paramp  SQUIDs.  There  are  a  few  small  specks  due  to  undeveloped  resist 
near  the  SQUID  loops,  but  they  are  away  from  the  junctions,  (g)  Picture  of  experimental 
setup  inside  the  dilution  refrigerator;  includes  banana  for  scale. 


had  been  observed  in  [15],  and  to  prevent  over-coupling  to  magnetic  field  via  the  flux 
focusing  effect,  whereby  the  effective  magnetic  field  in  the  gap  of  the  resonator  could  be 
two  orders  of  magnitude  larger  than  the  applied  magnetic  field  [7].  In  the  interest  of  time, 
an  in-shield  superconducting  magnet  was  not  obtainable,  and  for  the  first  paramp  tests  we 
instead  used  an  in-shield  copper  (non-superconducting)  coil  magnet,  which  turned  out  to 
provide  insufficient  magnetic  flux  through  the  SQUID  loops  (for  our  12  gm2  loop  area,  a 
flux  quantum  corresponds  to  .17  millitesla  at  the  SQUID)  before  significantly  heating  up 
the  fridge,  despite  thermal  anchoring  of  the  flux  bias  line  and  of  the  magnet  itself.  In  a 
second  fridge  cooldown,  we  used  a  superconducting  magnet  without  a  shield,  demonstrating 
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frequency  tunability  (but  with  fluctuations  due  to  lack  of  shield,  as  discussed  later).  The 
paramp  gain  measurements  in  this  thesis  are  from  the  first  cooldown  [with  a  shield],  so  we 
operated  the  paramp  at  a  fixed,  stable  value  of  ujq. 


(a)  (W  xio4 


Figure  4.14  Low-power  (—125  dBm  at  sample)  reflection  spectroscopy  of  paramp. 
(a)  Simplified  schematic  of  experimental  setup.  Attenuation  is  added  on  the  input  line  to 
attenuate  otherwise  significant  thermal  Johnson  noise  coming  from  high  temperatures,  (b) 
Unwrapped  raw  phase  from  PNA  network  analyzer,  (c)  The  same  data  with  the  phase- 
accumulation  background  slope  subtracted  off  (and  negated,  as  a  choice  of  convention  for 
the  phase  of  an  outgoing  vs.  incoming  wave). 


After  cooling  the  dry  dilution  refrigerator  down  to  base  temperature,  we  are  ready  to  test 
our  attempt  at  a  paramp!  A  picture  of  the  physical  setup  is  shown  in  Figure  4.13(g)  and  a 
schematic  of  the  measurement  circuitry  is  shown  in  Figure  4.14(a).  The  first  measurement 
performed  was  reflection  spectroscopy  with  a  network  analyzer.  A  single  tone  was  sent 
through  a  circulator  to  the  input  port  of  the  resonator  and  reflected  back  through  the 
circulator  to  an  output  port  that  was  connected  to  the  second  port  of  the  network  analyzer 
at  room  temperature.  As  expected,  our  niobium  resonator  has  an  internal  quality  factor 
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much  higher  than  its  coupling  quality  factor,  meaning  in  accordance  with  equation  (2.38) 
essentially  all  of  the  information  is  contained  in  the  reflected  phase.  As  shown  in  Figure 
4.14,  the  raw  phase  measured  by  the  network  analyzer  has  a  large  slope  as  a  function  of 
frequency  simply  due  to  the  natural  change  in  optical  path  length  through  the  various  coax 
cables  as  the  frequency  of  propagation  changes.  Luckily,  we  can  approximate  the  dispersion 
relation  of  our  coax  as  linear  and  subtract  off  this  large  background,  allowing  us  to  observe 
a  clear  resonance  at  7.26  GHz  as  expected  with  the  zero  applied  magnetic  flux  through 
the  SQUIDs  (the  predicted  value  of  7.2  GHz  at  maximum  SQUID  energy  was  surprisingly 
close).  By  fitting  the  phase  response  to  theory,  we  can  extract  a  quality  factor  of  198, 
meaning  271/(27 r)  ~  37  MHz.  This  is  surprisingly  close  to  the  predicted  value  of  35  MHz 
at  this  operating  frequency  (remember,  in  accordance  with  equation  (4.14),  71  depends  on 
the  resonance  frequency  in  addition  to  the  coupling  capacitance). 

4.2.1  Nonlinearity 

To  probe  the  Kerr  nonlinearity  of  our  resonator,  we  can  crank  up  the  probe  power  relative  to 
the  predicted  critical  input  power.  Figure  4.15  shows  reflection  spectroscopy  measurements 
as  a  function  of  probe  power  and  probe  frequency.  In  the  phase  response,  we  see  an  increase 
in  the  slope  of  the  phase  jump,  reaching  and  going  through  the  critical  point  as  described 
in  the  caption.  In  the  amplitude  response,  for  high  powers  we  can  start  to  see  a  small  dip 
on  resonance,  presumably  corresponding  to  the  onset  of  a  small  nonlinear  loss  term.  This 
term  is  presumed  to  be  small  because  of  the  small  dip  it  yields  along  with  the  fact  that  it 
was  not  needed  to  fit  the  functional  form  of  the  resonance  frequency  with  respect  to  probe 
power.  Namely,  in  both  the  phase  and  amplitude  response,  we  can  observe  a  small  dip  in 
the  system’s  resonance  frequency  as  a  function  of  increasing  probe  power  that  agrees  very 
well  with  the  quadratic  dependence  of  the  resonance  frequency  on  the  input  field  discussed 
in  Section  2.2.2,  as  shown  by  the  fit  to  theory  in  the  inset  of  Figure  4.15(b). 

The  probe  powers  shown  in  Figure  4.15  are  based  on  a  calibration  of  the  coax  lines  (and 
a  previous  calibration  of  the  lines  in  the  fridge).  From  the  onset  of  the  critical  point,  we 
deduce  a  critical  power  of  approximately  —93.5  dBm,  which  is  close  to  the  predicted  value 
[based  on  equations  (3.49)  and  (2.47)]  of  —89  dBm.  This  is  consistent  with  our  uncertainty 
in  the  SQUID  energy  and  in  the  performance  of  the  circulator  during  the  room-temperature 
calibration  versus  its  performance  while  cold. 
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Figure  4.15  Reflection  spectroscopy  of  the  paramp  with  a  single  tone,  (a)  Reflected 
phase  as  a  function  of  tone  power  (at  the  paramp)  and  frequency,  (b)  Vertical  cuts  of  (a) 
with  varying  powers  ranging  from  the  linear  regime  (magenta)  through  the  critical  point 
(green)  to  the  bistable  regime  (blue).  Note  that  the  resonance  frequency  decreases  with 
power,  as  expected.  The  extracted  resonance  frequency  along  with  a  fit  to  theory  is  shown 
in  an  inset,  (c)  Magnitude  of  the  reflection  coefficient  |T|,  relative  to  the  value  in  the 
linear  regime  (which  was  subtracted  from  each  of  the  columns).  Note  the  appearance  of 
loss  at  high  powers,  suggesting  a  (small)  amount  of  nonlinear  loss.  The  fit  of  frequency 
to  theory  also  suggests  that  the  nonlinear  loss  is  negligible,  (d)  Vertical  cuts  of  (a)  with 
various  powers. 


4.2.2  Parametric  Gain 

With  the  critical  power  known,  we  can  now  test  the  device’s  performance  as  a  parametric 
amplifier.  The  setup  for  this  measurement  is  also  shown  in  4.14,  using  two  microwave  signal 
generators  combined  at  room  temperature  and  a  spectrum  analyzer,  rather  than  a  network 
analyzer.  For  all  measurements,  the  spectrum  analyzer  (an  Agilent  E4405B  lent  to  us  by 
the  Houck  lab)  was  operated  with  its  minimum  bandwidth  resolution  (i.e.,  a  bin  for  power 
collection)  of  1  KHz.  Figure  4.16  shows  the  result  of  biasing  the  pump  near  its  critical 
power  and  frequency,  while  also  adding  a  much  smaller  test  signal  tone  detuned  10  KHz 
above  the  pump.  If  all  of  the  assumptions  we  made  in  Section  2.2  are  satisfied,  we  could 
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Figure  4.16  (a)  Illustration  of  input  and  output  spectrum  for  a  non-degenerate  paramp 
operated  in  the  ideal  linear  regime.  Figure  adapted  from  [15].  (b)  Observed  spectral 
densities  for  a  signal  power  60  dB  less  than  the  pump  power,  with  the  paramp  operating 
at  a  gain  of  about  20  dB.  (c)  Signal  input  50  dB  weaker  than  pump,  (d)  Signal  input  40 
dB  less  than  pump. 


expect  to  see  amplification  of  the  signal  and  the  creation  of  an  idler  tone  on  the  opposite 
side  of  the  pump.  This  is  precisely  what  is  observed.  For  a  signal  tone  approximately  60 
dB  below  the  pump  critical  power  (and  thus  the  pump)  [Figure  4.16(b)],  when  the  pump  is 
off  we  observe  a  very  small  power  spectral  density  just  above  the  room-temperature  noise 
floor  corresponding  to  the  small  signal.  When  the  pump  is  on,  we  observe  a  drastic  increase 
in  the  signal  power  and  also  the  appearance  of  an  intermodulation  tone  on  the  other  side  of 
the  pump  as  predicted  by  theory.  However,  if  we  increase  the  power  of  the  signal  relative 
to  the  pump  [Figure  4.16(c)  and  (d)],  other  peaks  in  the  spectrum  appear.  This  is  perhaps 
due  to  the  violation  of  the  stiff  pump  approximation,  given  that  the  signal  is  amplified  by 
20  dB,  although  the  pump  does  not  appear  to  be  visibly  depleted.  Further  experimental 
exploration  of  this  phenomenon  in  addition  to  a  theoretical  understanding  of  the  potential 
importance  of  terms  beyond  the  linearized  signal  approximation  are  worth  further  pursuit8. 

8A  hand-wavey  and  perhaps  incorrrect  intuition  behind  the  appearance  of  these  equally-spaced 
peaks  is  that  the  signal  starts  to  act  as  a  new  pump  for  the  original  pump,  causing  a  small  amount 
of  the  original  pump  to  appear  as  an  idler  relative  to  it  and  the  signal,  and  so  on  with  further 
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Figure  4.17  Maximum  gains  as  a  function  of  pump  bias,  (a)  Maximum  signal  gain 
Gs(10  KHz)  «  Gs(0)  as  a  function  of  pump  power  and  pump  frequency.  Panels  (a)  and 

(c)  share  the  same  colorscale.  (b)  Cuts  across  the  data  in  panel  (a),  (c)  Maximum  inter¬ 
modulation  gain  G,(10  KHz)  ss  Gi(0)  as  a  function  of  pump  power  and  pump  frequency. 

(d)  Cuts  accross  the  data  in  panel  (c). 


To  systematically  test  the  system’s  performance  as  an  amplifier,  we  can  measure  the 
maximum  signal  gain  as  a  function  of  pump  power  and  pump  frequency.  The  results  of  such 
a  measurement  are  shown  in  Figure  4.17.  In  this  experiment,  the  signal  is  kept  at  a  constant 
power  50  dB  below  the  pump  and  at  a  detuning  of  10  KHz  from  the  pump  frequency.  Since 
this  detuning  is  much  less  than  71 ,  measuring  the  signal  gain  will  give  a  good  approximation 
to  the  maximum  gains  Gs(0)  and  Gj( 0).  It  is  necessary  to  have  some  finite  detuning  u  7^  0, 
because  otherwise  we  would  be  operating  the  amplifier  in  degenerate  mode  and  would 
also  need  to  worry  about  the  relative  phase  of  the  signal  and  pump  (this  is  explored  in  a 
later  measurement).  We  see  a  maximum  gain  for  a  pump  power  of  approximately  —93.1 
dBm,  consistent  with  our  (somewhat  more  ambiguous)  spectroscopy  measurements.  We 
note  that  the  intermodulation  deamplification  in  Figure  4.17(d)  would  continue  to  decrease 
below  —10  dB  for  low  pump  frequencies  if  it  weren’t  for  the  noise  floor  becoming  important 


harmonics.  This  could  be  related  to  classical  intermodulation  distortion  [16]. 
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at  the  spectrum  analyzer. 

After  mapping  out  the  pump  bias  landscape,  we  can  explore  the  dependence  of  gain 
on  signal  frequency  with  a  fixed  pump;  i.e.,  the  bandwidth  of  the  amplifier.  We  do  this 
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Figure  4.18  (a)  and  (b):  Signal  and  intermodulation  gains  as  a  function  of  signal 
frequency,  for  four  different  pump  biases  ( fp,Pp ):  —96.3  dBm  at  7237.8  MHz  (magenta), 
—95  dBm  at  7230.9  MHz  (red),  —94  dBm  at  7225.9  MHz  (green),  and  —93.4  dBm  at  7222.5 
MHz  (blue).  Equivalently,  the  ratios  Pp/Pc rit.  in  increasing  order  are  {.48,  .65,  .81,  .94}. 
Each  trace  also  shows  a  fit  to  theory  (black  curves),  (c)  Exploration  of  the  dynamic  range 
of  the  amplifier  at  a  particular  bias  point  [near  the  bias  point  of  the  green  curves  in  (a) 
and  (b)]  such  that  the  signal  gain  is  20  dB.  (d)  Phase-dependence  of  gain  in  degenerate 
mode  (signal  frequency  the  same  as  pump  frequency,  originating  from  the  same  microwave 
source) . 


by  fixing  the  pump  at  one  of  the  four  points  of  maximum  gain  found  at  each  of  the  four 
pump  powers  explored  in  Figure  4.17.  Namely,  the  ratios  Pp/PCnt  in  increasing  order  are 
{.48,  .65,  .81,  .94}.  The  resulting  gain  curves,  shown  in  Figure  4.18(a)  and  (b),  are  symmetric 
about  the  pump  as  expected,  and  fit  to  theory  (black  curves)  quite  well.  We  can  clearly 
see  that  the  bandwidth  of  the  amplifier  decreases  with  the  maximum  gain.  From  the  fitted 
curves,  we  extract  signal  gain-bandwidth  products  of  ^  \J G{C)B^I (27 r)  ss  {19,20,22,8} 
MHz.  The  first  three  of  these  products  agree  surprisingly  well  with  our  high-gain  prediction 
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(2.70),  which  says  this  should  simply  equal  7i/(27r)  ~  18.5  MHz.  The  discrepancy  for  the 
highest  gain  could  possibly  be  due  to  the  onset  of  pump  depletion.  Either  way,  we  have 
attained  our  designed  performance  of  2  MHz  of  bandwidth  at  25  dB  gain. 

Figure  4.18(c)  shows  the  power  spectral  density  (at  the  spectrum  analyzer)  of  the  am¬ 
plified  signal  as  a  function  of  input  signal  power,  with  the  paramp  operated  with  a  gain  of 
20  dB.  We  see  that  the  paramp  performs  well  as  a  linear  amplifier  for  circuit  QED  exper¬ 
iments,  having  a  linear-regime  dynamic  range  of  at  least  30  dB  (the  noise  floor  limits  the 
probing  of  signals  below  the  lower  limit  of  this  range).  This  range  is  limited  by  saturation 
of  the  amplifier  at  high  signal  input  powers  due  to  the  onset  of  pump  depletion  (in  the  plot 
shown,  this  corresponds  to  the  input  signal  being  about  30  dB  less  than  the  pump  [and 
then  amplified  by  20  dB  to  10  dB  less  than  the  pump]). 

Regarding  the  noise  floor  at  the  spectrum  analyzer,  it  was  usually  found  to  be  —73  to 
—70  dBm  (although  on  the  first  day  of  measurements  it  was  only  —80  dBm  for  unknown 
reasons).  This  drops  to  —95  dBm  when  the  output  of  the  room  temperature  amplifiers  is 
disconnected  from  the  spectrum  analyzer.  If  we  instead  leave  the  analyzer  connected  to  the 
output  of  the  room  temperature  amplifiers  and  instead  unplug  the  input  to  the  first  room 
temperature  amplifier,  the  noise  floor  is  still  —89  dBm.  To  extend  paramp  measurements  to 
be  able  to  observe  lower  deamplfications,  one  could  try  and  bypass  the  room  temperature 
amplifiers  altogether.  But  with  the  paramp  operated  at  25  dB  gain,  we  can  then  see  the 
amplified  signal  on  the  analyzer  but  not  the  bare  signal  without  the  pump  turned  on.  This 
means  that  the  signal  itself  must  have  been  below  the  noise  floor. 

Finally,  we  test  our  prediction  for  phase-dependent  gain  when  the  paramp  is  operated  in 
degenerate  mode.  To  do  this,  we  operate  the  pump  near  the  same  bias  point  (19  dBm  gain), 
but  provide  the  signal  through  the  same  source  as  the  pump  via  a  splitter.  The  phase  of  the 
pump  relative  to  the  signal  is  then  changed  by  routing  the  pump  path  through  a  Colby  PDL- 
10A  programmable  delay  line,  which  can  provide  delays  (i.e.,  path  length  differences)  for 
microwaves.  Because  the  “signal”  and  “pump”  are  at  the  same  frequency,  in  single- frequency 
degenerate  mode  we  must  define  the  signal  gain  as  the  gain  in  sideband  modulation  power 
when  the  signal  is  amplitude-modulated  at  a  certain  frequency.  Namely,  we  amplitude- 
modulate  the  signal  path  at  20  KHz  via  a  Marki  mixer,  and  measure  the  sideband  power 
spectral  density  with  the  pump  on  and  off.  The  degenerate  gain  is  then  defined  by  these 
two  relative  powers.  Figure  4.18(d)  shows  the  results  of  this  experiment.  We  can  clearly  see 
a  strong  periodic  dependence  on  the  PDL  delay,  which  is  ideally  proportional  to  the  path 
length  difference  (and  hence  phase  difference)  of  the  pump  relative  to  the  signal.  However, 
there  are  some  imperfections  due  to  the  imperfect  nature  of  the  delay  line  and  the  fact  that 
changing  the  delay  also  sometimes  slightly  changes  the  pump  power.  It  is  quite  interesting 


4.2  Amplifier  Implementation  and  Characterization 


82 


to  note  that  the  maximum  gain  is  greater  than  the  nominal  G(10  KHz)  ~  G(0)  =  19  dB, 
and  reaches  up  to  25  dB.  This  is  just  as  predicted  by  the  quadrature  squeezing  relations 
(2.95).  The  de-amplified  quadrature  should  have  a  correspondingly  lower  gain,  but  we  are 
unable  to  oberve  the  full  dip  due  to  both  the  finite  resolution  of  the  delay  line  and  also  the 
noise  floor  at  the  spectrum  analyzer. 

In  a  second  fridge  cooldown,  we  removed  the  magnetic  shield  from  the  sample  and 
instead  used  a  larger  external  superconducting  magnet  below  the  base  plate  of  the  fridge. 
This  allowed  us  to  thread  enough  flux  through  the  SQUIDs  without  heating  the  fridge,  so 
that  we  could  observe  the  frequency-tunability  of  the  paramp.  The  results  of  this  experiment 
are  shown  in  Figure  4.19.  The  resonance  frequency  as  a  function  if  applied  magnetic  field 


Phase  (rad.) 


0  2  4  6  8  10 


Applied  Magnetic  Field  (piT) 

Figure  4.19  Reflected  phase  of  a  single  tone  in  the  linear  regime,  as  a  function  of 
applied  magnetic  field.  The  quasi-periodic  resonance  indicates  clear  tunability  of  u> o . 

is  not  strictly  periodic  [as  would  be  predicted  for  a  single  SQUID  by  equation  (3.20)], 
because  there  are  in  fact  four  SQUIDs  that  are  not  exactly  identical.  The  quasi-periodic 
pattern  is  therefore  to  be  expected.  We  note  that,  based  on  a  previous  calibration  of  the 
superconducting  magnet,  the  applied  magnetic  field  needed  to  tune  through  a  period  is 
roughly  50  times  less  than  what  would  be  expected  without  considering  the  flux  focusing 
effect.  This  means  that  our  attempts  to  reduce  flux  focusing  by  carving  out  the  Nb  ground 
plane  near  the  SQUIDs  had  at  most  a  factor  of  ~  2.  Note  that  the  zero- field  value  here 
leads  to  a  different  resonance  frequency  than  the  previous  cooldown  with  the  magnetic 
shield  and  no  applied  field.  This  is  because  during  the  second  cooldown,  magnetic  flux 
was  trapped  in  the  superconductor  due  to  lack  of  shielding,  leading  to  a  flux  offset.  The 
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maximum  resonance  frequency  observed  agrees  with  what  we  designed,  and  the  desired 
operating  point  of  6.85  GHz  is  at  a  point  with  a  good  tradeoff  between  tunability  and 
stability.  We  note  that  because  there  was  no  internal  magnetic  shield  at  the  sample  during 
the  cooldown,  the  resonance  frequencies  drifted  slightly,  due  to  either  external  or  internal 
(i.e.,  superconducting  vortices  from  trapped  flux)  magnetic  noise.  Namely,  for  a  resonance 
frequency  of  7.1  GHz,  the  resonance  drifted  by  a  about  5  MHz  on  the  timescale  of  a  minute, 
whereas  at  6.2  GHz,  the  resonance  drifted  by  about  100  MHz  on  this  same  timescale.  It 
will  therefore  be  necessary  to  implement  an  in-shield  superconducting  magnet  to  obtain  a 
tunable  paramp  that  is  also  stable. 

4.3  Future  Directions 

This  thesis  has  performed  a  theoretical  analysis  of  parametric  amplification  and  has  pre¬ 
sented  preliminary  data  demonstrating  the  successful  implementation  of  a  paramp  for  circuit 
QED  based  on  a  superconducting  resonator  inheriting  a  Kerr  nonlinearity  from  an  array  of 
SQUIDs.  The  paramp  characterization  measurements  performed  in  section  4.2  were  only 
started  shortly  before  the  completion  of  this  thesis,  so  there  are  still  further  measurements 
that  would  be  worth  pursuing,  including  a  more  detailed  analysis  of  the  data  shown  and  a 
verification  that  the  second  paramp  device  also  works.  It  would  also  be  desirable  to  see  how 
quantum-limited  the  added  noise  of  the  amplifier  is.  The  added  noise  could  be  probed  for 
example  through  a  so-called  Y-factor  measurement  in  which  two  calibrated  noise  sources 
are  input  to  the  paramp  [16,42],  Implementing  a  superconducting  magnet  (or  perhaps  an 
on-chip  flux  bias  line)  within  the  magnetic  shield  is  the  next  design  step  if  we  wish  to  have 
a  tunable  and  stable  amplifier  to  use  in  actual  circuit  QED  experiments.  The  most  ele¬ 
gant  and  relatively  localized  solution  would  probably  be  modifying  an  octobox  by  shaving 
an  indent  around  its  circumference,  allowing  many  layers  of  thin  superconducting  wire  to 
be  wrapped  around  it.  Demonstrating  that  the  paramp  can  be  used  to  read  out  a  super¬ 
conducting  qubit  with  a  decreased  number  of  averages  would  then  be  a  proof-of-principle 
demonstration  of  paramp-based  readout.  A  simple  and  effective  readout  scheme  to  use  was 
recently  demonstrated  by  Riste  et  al,  whereby  a  small  dispersive  readout  signal  leads  to  a 
large  change  in  the  phase  of  the  reflected  pump  when  the  paramp  is  operated  in  single-tone 
degenerate  mode  near  its  sensitive  critical  point  [69]. 

After  achieving  readout  of  a  superconducting  transmon  qubit  with  the  paramp  (as  was 
done  in  [69]),  we  would  like  to  develop  paramps  to  be  used  in  novel  experiments  involving 
quantum  dots.  Other  spin  qubit  projects  in  our  lab  would  also  benefit  from  such  a  Josephson 
parametric  amplifier;  this  will  require  lower  ranges  of  operating  frequencies  (on  the  order 
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of  1  GHz)  and  would  be  an  interesting  engineering  challenge  (although  presumably  not  as 
challenging  as  conversely  going  to  higher  frequencies).  Realistic  progress  towards  useful 
paramps  in  our  lab  is  now  in  the  foreseeable  future. 


Appendix  A 

Quantization  of  the 
Electromagnetic  Field  in  a 
Perfectly  Conducting  Cavity 


We  briefly  perform  a  variant  of  a  standard  method  of  field  quantization  based  on  Maxwell’s 
equations  (see  for  example  [39])  that  is  accessible  to  someone  with  an  understanding  of 
undergraduate  quantum  mechanics  and  electromagnetism.  In  the  main  text  we  will  give 
a  rather  different  and  specific  treatment  of  cavity  quantization  by  analyzing  the  modes  of 
the  particular  type  of  resonator  used  in  making  a  cQED  paramp.  However,  the  derivation 
in  this  section  is  quite  general  and  will  help  us  to  understand  where  the  Jaynes- Cummings 
Hamiltonian  comes  from  in  the  next  section.  We  start  by  examining  the  electromagnetic 
vector  potential  satisfying 


B(r,  t)  =  V  x  A(r,  t)  (A. la) 

E(r,f)  =  ~dAgt,tS>  ~  V(Kr^)>  (A. lb) 

where  4>  is  the  scalar  potential.  As  studied  in  an  undergraduate  course  in  electromagnetism 
[70],  choosing  the  Coulomb  gauge 

V  •  A  =  0,  (A. 2a) 

V2</>  =  --  (A. 2b) 

eo 
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and  combining  with  Maxwell’s  equations  yields  that  the  vector  potential  in  vacuum  obeys 
the  wave  equation 


V2A(r,  t)  = 


1  <92A(r,  t) 


(A.3) 


c2  dt2 

The  theory  of  harmonic  analysis  tells  us  that  it  is  possible  to  express  any  solution  of  (A.3) 
as  a  linear  combination  of  solutions  of  the  form  A(r,  t)  =  u(r )T(t).  Substituting  such  a 
separated  solution  into  (A.3)  yields  for  each  Cartesian  component  of  u(r) 


VV*}(r)  m 

«W(r)  c2T(t ) 

Since  the  left-hand  side  depends  only  on  r  and  the  right-hand  side  only  on  t,  each  side  must 
be  equal  to  a  mutual  constant.  For  an  approximately  non-leaky  cavity,  we  expect  that  this 
constant  should  be  a  negative  real  number,  say  —k2.  Defining  ujk  =  c\k\ ,  this  implies  that 
Tk(t)  =  e±iuJkt  and  that 

(v2  +  ufc(r)  =  0.  (A. 5) 

For  a  field  confined  to  well-defined  cavity  of  finite  volume  V,  the  spectrum  will  be 
infinite  but  discrete,  and  the  expression  for  A  takes  the  form  of  a  unique  discrete  sum  over 
a  complete  set  of  orthonormal  modes.  The  Coulomb  gauge  condition  (A. 2a)  takes  away  a 
degree  of  freedom,  so  that  we  expect  the  spatial  vector  eigenfunctions  of  (A. 5)  to  have  two 
independent  orthogonal  mode  polarizations,  indexed  by  //  =  1,  2.  We  define  these  functions 
to  be  unitless  and  normalized  so  that 


1 

V 


ufc,M(r)  '  ufc',M'(r)  dr 


fikk'dfi/i'  ■ 


(A.6) 


The  exact  nature  of  these  functions  depends  on  the  nature  of  the  cavity  boundary.  For 
example,  in  a  rectangular  box  with  perfectly  conducting  walls,  they  take  the  familiar  Fourier 
form  efc)/ie*k'r,  where  is  a  unit  polarization  vector  [necessarily  perpendicular  to  the 
wavevector  k  by  (A. 2a)].  More  generally,  (A. 2a)  dictates  that 


V  •  ufcj/1(r)  =  0.  (A. 7) 

From  now  on,  we  will  suppress  the  double  index  on  Uk,,n  in  favor  of  the  more  compact 
notation  u.m  (m  stands  for  a  single  mode).  An  arbitrary  real  solution  for  A(r,  t)  can  now 
be  expressed  as 

A(M)  =  [am(*)um(r)  +  <4(f)u^(r)] , 


m 


(A.8) 
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where  the  am  are  complex  mode  amplitudes  with  time  dependence  absorbed  into  them  and 
with  SI  units  of  V-s-m-1.  Substituting  (A. 8)  into  (A.l)  yields 


E(r) f)  =  J2iuJm  [am(f)um(r)  -  a*m{t)u*m(r)} 
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B(r ,t)  =  ^2  [am(i)V  x  um(r)  +  a*n(t)V  x  u^(r)]  (A.9) 

m 

The  classical  held  Hamiltonian  is  the  electromagnetic  energy 

U  =  |  f  (E2  +  c2B2)  dr  =  Ve0  ^  u2m{a*mam  +  ama*J,  (A.10) 

where  we  have  used  (A. 6)  and  (A. 7),  and  in  anticipation  have  withheld  from  simplifying 
the  expression  using  the  usual  commutative  property  of  complex  multiplication1.  We  can 
notice  a  similarity  between  the  form  of  (A.  10)  and  an  infinite  set  of  uncoupled  harmonic 
oscillators  of  the  kind  studied  in  elementary  quantum  mechanics  [50].  We  postulate  that 
the  creation  and  annihilation  operators  should  be  defined  by  the  promotion 


h 


2uJmV  €q 


°m(*) 


]/  2umVe/m 


and  obey  the  canonical  bosonic  commutation  relation 


(A.ll) 


6 


mm'  • 


(A.12) 


The  coefficients  in  (A.ll)  are  chosen  so  that  using  (A.12),  the  quantum  version  of  (A.  10) 
can  be  written  in  the  standard  form 


/ 1  cavity  —  ^  hwr i 


Z _ / 

in 


dmdm  +  2 


(A.13) 


Each  harmonic  oscillator  excitation  corresponds  to  a  photon,  and  as  usual  ( dlndm )  is  the  ex¬ 
pected  number  of  photons  of  a  given  mode  occupying  the  resonator.  The  energy  eigenstates 
are  the  number  states  (also  known  as  Fock  states)  with  a  definite  number  of  excitations  nm 

1More  specifically,  for  the  magnetic  part  we  have  used  the  expansion  (V  x  u*)  •  (V  x  u)  = 
V  •  [u*  x  (V  x  u)]  +  u*  •  [V  x  (V  x  u)]  and  then  Stokes’  theorem  to  convert  the  integral  of  the  former 
term  into  a  boundary  term  that  vanishes  for  a  perfectly  conducting  cavity  wall  at  which  u  x  dS  =  0. 
The  integral  of  the  latter  term  becomes  fv  u*  •  [V  x  (V  x  u)]dr  =  (using  (A. 7))  fv  u*  •  —  V2u  = 
(using  (A. 5)  and  (A. 6))  w^/c2. 


in  each  mode,  leading  to  energies 


E  —  ^  ^  +  1/2). 

m 


(A.14) 
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The  Jaynes- Cummings 
Hamiltonian 


B.l  Derivation 


The  problem  we  now  want  to  address  is  this:  how  does  a  single  atom  interact  with  photons 
in  a  resonator?  Under  the  approximations  we  are  about  to  make,  there  turns  out  to  be  a 
rather  simple  Hamiltonian  describing  the  evolution  of  the  joint  atom-cavity  system. 

Any  atom  is  a  many-electron  system  bound  by  an  effective  atomic  potential  I4,tm(r)- 
The  electronic  wavefunction  of  an  isolated  atom  will  have  a  countable  number  of  bound 
energy  eigenstates  |j)  with  energy  Ej  forming  a  complete  orthonormal  basis  of  its  Hilbert 
space,  so  that  the  Hamiltonian  takes  the  form  ")T  Ej\j)(j\.  If  we  truncate  the  atom’s  energy 
levels  to  the  lowest  two  in  order  to  obtain  a  two-level  system  (valid  when  the  spectrum  is 
anharmonic  enough  that  the  desired  transition  can  be  isolated),  we  obtain,  after  shifting 
the  zero  of  energy  to  the  average  energy,  the  simple  atomic  Hamiltonian 


H&tm  —  2 


1  0 
0  -1 


(B.l) 


in  the  { 1 1) ,  1 0) }  basis,  where  huja  is  the  energy  difference  between  the  two  energy  eigenstates. 
We  now  have  the  atomic  Hamiltonian  and  the  cavity  Hamiltonian  (A.  13),  and  would  like  a 
term  Hin t  describing  atom-photon  interactions. 

The  classical  minimal  coupling  Hamiltonian  describing  the  interaction  of  an  electron 
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with  the  electromagnetic  field  is  [61] 

"Hem  =  (P  +  eA(r,  t))2  -  e(f>(r,  t),  (B.2) 

where  — e  is  the  electron  charge  and  A  is  the  vector  potential.  It  is  easily  verified  that 
Hamilton’s  equations  of  motion  with  Hem  yield  the  correct  classical  Lorentz  force  expres¬ 
sions.  We  now  use  the  the  semi-classical  version  of  (B.2),  under  which  the  electron  is  treated 
quantum- mechanically  (neglecting  spin)  but  the  electromagnetic  field  is  not,  to  write  down 
the  Hamiltonian  for  an  atomic  electron 

Hc\  =  (  V  -  —  A(r,t)J  -  e0(r,t)  +  I4tm(r)  (B.3) 

(the  electric  field  due  to  the  atom  is  contained  in  I4tm)-  We  know  from  classical  electro¬ 
magnetism  [or  infer  from  (A.l)]  that  the  electric  and  magnetic  fields  are  invariant  quantities 
under  the  gauge  transformation 


A  ->  A'  =  A  +  VA 

</>-></>'  =  *--&  (B-4) 


for  some  scalar  A.  However,  it  is  not  obvious  that  such  a  transformation  is  consistent  with 
Schrodinger’s  equation.  Directly  applying  such  a  transformation  yields  a  very  different 
looking  time  evolution,  but  we  can  try  to  restore  the  original  form  by  performing  a  local 
phase  transformation  on  the  position  space  wave  function  of  the  form 


(B.5) 


Plugging  this  into  (B.3)  shows  us  that  the  Schrodinger  equation  retains  its  form  provided  we 
choose  x(r,  t)  =  §A(r,  t).  Observables  will  also  remain  unchanged  provided  we  transform  a 
given  Hermitian  operator  A  to  UAU t,  with  U  =  eix^r,t\ 

With  the  semi-classical  picture  in  hand,  we  choose  again  to  start  in  the  Coulomb  gauge 
(A. 2).  We  consider  the  interaction  with  photons  in  vacuum  with  no  charge  density  in  the 
vicinity  other  than  that  due  to  the  atom,  so  we  may  take  the  instantaneous  potential  </>(r,  t) 
due  to  outside  sources  to  be  0  as  we  did  in  our  field  quantization.  The  Schrodinger  equation 
becomes 

=  ~2 nT  (V  _  +Vatm(r)  (B-6) 

where  we  have  made  the  dipole  approximation  [38]  A(r,  t)  ~  A(ro,  t),  where  ro  is  the  center 
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of  the  atom.  This  approximation  is  valid  if  the  wavelength  of  the  mode  is  large  compared 
with  the  linear  size  of  the  bulk  of  the  atomic  wavefunction.  Note  that  we  can  eliminate  the 
A(ro,f)  term  in  (B.6)  by  making  the  gauge  transformation  (B.4)  with  A  =  —  A(ro,t)  •  r, 
along  with  the  required  phase  transformation 

Y>(r, t)  — >  i//(r,t)  =  e_«" A(ro,t'l'rV;(r, t).  (B.7) 


Substituting  these  transformations  into  (B.6)  yields,  after  rearrangement, 


ih 


W(r,t) 

dt 


2m 


V2  +  14tm(r) 


er  •  E(r0,  t) 


(B.8) 


where  we  have  used  (A.  la)  for  the  electric  field.  Summing  over  all  of  the  atomic  electrons 
(indexed  by  a),  we  arrive  at  the  semi-classical  interaction  Hamiltonian  in  this  gauge, 


Hint  =  •  E(i-0,f).  (B.9) 

a 

We  can  rewrite  e£ara  as  =  e|*)(*|  £ara|j>01  =  £ijdy|*>01  =  EijVli-  where 

d ij  =  e(i|£ara|j)  =  £a  fff  ^>*(ra)era^fj(ra)d3ra  is  the  transition  dipole  moment.  Using 
(A. 9)  and  (A.  11),  we  obtain  after  shifting  each  am  by  an  appropriate  phase  factor  (this  is 
just  a  convention  and  leaves  the  cavity  Hamiltonian  (A.  13)  unchanged)  the  electric  field 
operator 

e  =  (B-10) 

Combining  these  results,  we  therefore  expect  a  fully  quantized  (non-relativistic)  atom- 
photon  interaction  Hamiltonian  of 


H'mt  —  ti  ^  ^  ^  ®  (®m  “b  Q-m);  (B.ll) 

i,j  in 

where  (jihm  =  —  dy  are  dipole  coupling  factors.  On  timescales  small  compared  to 

the  transit  time  Ttransit  of  the  atom  through  the  cavity,  we  can  treat  the  position-dependent 
coupling  factors  as  constants.  This  is  not  always  valid  in  CQED,  but  is  true  for  the  sta¬ 
tionary  artificial  atoms  of  cQED. 

In  general,  the  coupling  strengths  increase  as  the  volume  decreases,  as  could  be  seen 
directly  for  a  rectangular  box  cavity  by  plugging  in  the  aforementioned  Fourier  form  of 
um.  More  intuitively,  this  is  due  to  the  fact  that  at  a  given  photon  frequency,  a  smaller 
volume  leads  to  a  higher  energy  density  because  the  total  energy  contained  in  the  cavity 
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must  take  on  the  same  quantized  value.  Circuit  QED  is  usually  concerned  with  systems  in 
the  strong  coupling  regime  for  which  g  3>  k,  7j_.  For  simplicity,  we  take  d,j  to  be  real,  so 
that  gij:m  =  gji.m.  =  9m  when  there  are  only  two  atomic  levels.  Thus,  for  a  two-level  atom 
(qubit),  we  have 

ll i 1 1 1  —  ft  ^  ^  9m(&+  T  )(®m  "I-  (H.12) 

m 

where  we  have  defined  the  atomic  raising  and  lowering  operators  <r+  =  jl)(0|  and  cm  =  |0)(1| 
and  dropped  the  explicit  tensor  product  symbol.  If  electromagnetic  signals  are  sent  with 
frequencies  and  polarizations  intended  to  excite  only  one  photon  mode,  say  one  of  the  the 
lowest  ones  of  frequency  ojr,  we  only  need  to  consider  one  mode  in  the  sum  (B.12).  In 
the  microwave  cavities  of  cQED,  the  one-mode  approximation  is  certainly  possible  for  low 
electronic  noise  and  cryogenic  temperatures:  at  T  =  20  mK  and  ojr  =  2ir  x  7  GHz,  the 
Bose-Einstein  distribution  yields  a  thermal  population  of 


nr  = 


ehuJr/kT  _  2 


10“ 


(B.13) 


for  the  lowest  mode  and  even  smaller  numbers  for  higher  modes.  Thus  we  may  take  m  =  1 
and  drop  all  other  terms  (along  with  the  subscript).  Finally,  we  make  the  rotating  wave 
approximation  [38],  which  keeps  only  the  terms  proportional  to  a+a  and  <T_cd  that  conserve 
excitation  number  (and  also  energy,  when  the  atom  transition  frequency  matches  that  of 
the  photon).  For  example,  cr+a  represents  the  excitation  of  the  atom  with  an  annihilation 
of  a  photon.  The  name  given  to  this  approximation  stems  from  the  fact  that  in  the  in¬ 
teraction  picture  Hj  =  eiHca-vityt^hH(]fe~iHcavityt/h,  the  non-energy-conserving  terms  oscillate 
rapidly  and  average  out  to  zero.  Combining  our  interaction  Hamiltonian  with  the  atom 
and  cavity  Hamiltonians,  we  arrive  at  the  well-known  Jaynes-Cummings  Hamiltonian  of 
quantum  optics  for  an  atom  interacting  with  a  single  mode  of  the  electromagnetic  field: 


Hjc 


hujr  +  0  +  l>hua<jz  +  hg(a+a  +  <7_at'), 


(B.14) 


where  we  have  dropped  the  hats  on  operators  and  the  identity  operator  is  implied  for  the 
photon  field  in  the  purely  atomic  term  and  vice  versa.  We  have  arrived  at  this  result  using 
many  approximations,  but  its  predictions,  some  of  which  we  will  now  explore,  have  been 
extremely  well-verified  by  many  CQED  and  cQED  experiments. 
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B.2  Coupling  Phenomena 


Analytic  Eigenstates  and  Eigenenergies 

We  now  discuss  the  eigenstates  and  eigenergies  of  the  coupled  atom-cavity  system  under 
the  Jaynes- Cummings  Hamiltonian.  Still  treating  our  atom  as  a  two-level  system,  we  will 
use  the  basis  states  \g)\n)  =  \g,n)  and  |e)|n)  =  |e,n),  which  denote  atom  in  the  ground  and 
excited  state  (respectively)  with  n  photons  in  the  cavity.  We  should  suspect  that  these  basis 
states  will  no  longer  be  eigenstates  of  the  whole  Hamiltonian  due  to  the  interaction  term. 
We  can  easily  solve  for  the  eigenstates  by  first  writing  the  Jaynes- Cummings  Hamiltonian 
in  matrix  form.  Since  iJjc  only  couples  neighboring  atom  and  photon  levels  (i.e,  swaps 
excitations  between  \g)\n)  and  |e)|ra—  1)),  the  matrix  takes  a  block  diagonal  form.  Defining 
the  detuning  A  =  uia  —  ujr  and  plugging  into  Hjc  (we  suppress  the  hats  on  Hamiltonians 
from  now  on),  we  see  that  the  nth  block  (n  >  1)  is 


Hin )  =  h 


nur  —  %  g\fn 
g\/n  nur  +  ^ 


(B.15) 


in  the  n-excitation  subspace  basis  {|g)|n),  |e)|n—  1)}.  We  can  diagonalize  the  entire  Hamil¬ 
tonian  by  diagonalizing  each  of  the  two-dimensional  H^n\  This  simple  calculation  yields 
eigenenergies 

E±>n  =  nhujr  =t  —  y/ 4  g2n  +  A2.  (B.16) 

For  n  =  0,  we  have  the  ground  state  energy  of  —  ^  corresponding  to  the  zero  excita¬ 
tion  eigenstate  |g}|0)  (which  wasn’t  formally  included  in  our  blocks).  Since  these  are  two- 
dimensional  subspaces,  we  may  find  an  angle  9n  such  that  the  orthonormal  eigenstates 
are 


+,n)  =  sin0n|g)|n)  +  cos9n\e)\n  -  1) 

—  ,n)  =  cos9n\g)\n)  —  sin0„|e)|n  —  1).  (B.17) 


To  find  the  form  of  9n ,  note  that  we  require  H^\+,n)  =  E+y1l\  +  ,  n).  Equating  the  compo¬ 
nents  of  \g)\n)  then  yields 


tan  9n 


ZgVh 

A  +  \J  4g2n  +  A2 


1 

=  -  arctan 
2 


2  tan  6 
1— tan2  0 ' 


(B.18) 


where  we  have  used  the  identity  tan  29 
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The  Resonant  Regime 

At  zero  detuning  (the  so-called  resonant  regime),  the  eigenstates  are  equal  superpositions 
of  atom  and  photon  states, 


H’±)  =  (\g)\n)±\e)\n-1))/V2  (B.19) 

Looking  at  (B.16),  in  this  regime  the  once  degenerate  subspace  of  \g)\n)  and  |e)|n  —  1)  is 
split  into  two  energies  differing  by  2 hg^/n,  a  phenomenon  known  as  vacuum  Rabi  splitting. 
Spectroscopically  measuring  the  lowest  such  splitting  is  a  common  way  of  experimentally 
obtaining  the  value  of  the  coupling  strength  g.  We  can  also  now  predict  the  fundamental 


ln> 

13} 


12) 

ID 

10) 


lg> 


Figure  B.l  Jaynes-Cummings  energy  level  diagram  (in  angular  frequencies)  in  the 
resonant  regime  A  =  0,  where  we  observe  vacuum  Rabi  splitting  (uncoupled  degenerate 
blue  and  red  lines  become  split  purple  lines).  |n)  refers  to  photon  number  and  |c/)  and  |e) 
correspond  to  the  atom  in  the  ground  or  excited  state.  Figure  adapted  from  [7]. 


CQED  phenomenon  of  vacuum  Rabi  oscillations.  This  is  obvious  at  zero  detuning  from 
the  nature  of  our  eigenstates;  starting  in  =  |e,  0),  the  state  will  oscillate  to  and  from 
|g,  1),  since  these  states  are  the  standard  |+)  and  |— )  states  on  the  Bloch  sphere.  More 
generally  at  any  detuning,  if  we  start  in  the  state  |e)|n  —  1),  the  standard  time  evolution  of 
the  separated  Schrodinger  equation  applied  to  the  n-excitation  subspace  yields  (after  using 
(B.17)  to  express  |e)|n  —  1)  in  terms  of  the  eigenstates  |  —  ,n)  and  |+,n)) 

ip(t)  =  cos9ne~lE+’n^h\ +,n)  —  sin  9ne~lE~'n^h\  — ,  n).  (B.20) 


Plugging  in  the  eigenenergies  (B.16)  and  taking  the  squared  magnitude  of  the  | g,n)  com¬ 
ponent  yields  the  probability  Pg(t)  of  the  excitation  having  transferred  to  a  cavity  photon: 


Pg{t)  =  | cos  9n  sin  9ne  lQt  —  sin0ncos0„e*nt| 


2 


(B.21) 
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where  Q  = 
sin  arctan  x 


\J g2n  +  (^)2.  Using  the  trigonometric  identities  cos  6  sin  6 
=  -7==,  we  obtain  the  expression 


^  sin  20  and 


Pg(*)  =  ^^2  Sil1'2  Qt  (B'22) 

Thus,  the  excitation  is  swapped  between  atom  and  cavity  at  the  vacuum  Rabi  frequency 
2fL  It  is  also  interesting  to  note  that  the  coefficient  in  front  of  the  oscillation  term  is  a 
Lorentzian  as  a  function  of  the  detuning,  which  is  a  very  special  instance  of  commonly 
measured  Lorentzian  absorbtion  spectra.  Unlike  spontaneous  emission  into  a  continuum 
of  vacuum  modes  (which  is  essentially  the  result  of  the  destructive  interference  of  many 
modes),  this  is  a  coherent,  reversible  phenomenon  that  occurs  when  the  atom  interacts 
with  a  single  quantized  mode  of  the  electromagnetic  field  in  a  well-defined  cavity.  As  such, 
it  has  no  obvious  classical  counterpart.  When  measuring  the  vacuum  Rabi  splitting,  the 
finite  lifetimes  of  the  excitations  due  to  decoherence  manifest  themselves  as  a  finite  width 
in  these  energy  levels  in  accordance  with  the  energy-time  uncertainty  principle.  Since  in 
the  resonant  regime  excitations  are  half  atomic  excitations  and  half  photons,  we  might 
guess  that  the  decay  rate  of  these  states  should  be  T  «  7x2+K .  If  we  are  in  the  strong 
coupling  regime,  the  splitting  is  well-defined  and  has  sharp  peaks,  with  many  vacuum  Rabi 
oscillations  occurring  before  the  decay  of  the  excitation. 


The  Dispersive  Regime 

A  now  paradigmatic  way  of  utilizing  cQED  for  quantum  computation  is  to  operate  in  the 
dispersive  regime ,  where  the  detuning  is  large  compared  to  g  (but  not  so  incredibly  large 
that  our  rotating  wave  approximation  breaks  down)  so  that  photons  are  not  absorbed  by  the 
atom.  This  regime  provides  us  with  a  means  of  state  readout  (D5)  by  providing  an  effective 
qubit  state-dependent  shift  in  the  cavity  resonance  frequency,  which  can  be  detected  by 
probing  the  cavity  with  microwaves,  and  also  provides  a  way  to  couple  distant  qubits  via 
the  cavity.  In  this  regime,  Hmt  will  be  weak,  as  even  the  terms  cr+a  and  a _cd  don’t  conserve 
energy.  Rather  than  simply  Taylor  expanding  (B.16)  and  (B.17),  it  will  be  useful  to  find  a 
unitary  change  of  basis  U  such  that  the  transformed  Hamiltonian  H'K,  =  UH^clJi  eliminates 
the  coupling  term  to  first  order.  We  will  now  motivate  the  transformation  commonly  stated 
in  the  cQED  literature  [8].  We  know  that  in  this  regime  g/ A  <C  1,  so  that  if  we  express 
U  =  e^x ,  where  X  is  antihermitian,  we  can  write  Hj(,  =  e%x Hjce~^x  and  Taylor  expand 
for  small  g/ A  (recognizing  and  factoring  out  a  commutator  each  time  we  take  a  derivative) 
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so  that 

H'j c  =  HjC  +  ^[X,HJC\  +  ^  (|)2  [X,  [X,HjC]}  +■■■•  (B.23) 

We  want  the  second  term  to  eliminate  Hm t  from  the  first  term.  Since  Hm t  is  first  order  in  g, 
we  thus  require  an  X  which  is  zeroth  order  in  g  and  satisfies  [X.  Hcav ;ty  +  Haim\  =  —Hmt. 
Dividing  through  by  ^  yields 

ur[X,  a^a]  +  -ua[X,  az\  =  ay((T+a  +  cr_a^)  —  ay(<r+a  +  cr_a^).  (B.24) 


Using  the  commutation  relations  [a, a)  a]  =  a,  [a) ,0) a]  =  —  a*",  [a+,az]  =  — 2<r+,  and 
[(j_,  az ]  =  2<r_,  we  see  by  inspection  that  this  equation  is  satisfied  by  X  =  a+a  —  cr_a^,  so 


that 


U  =  exp 


'  9_ 
.A 


(cr+a 


—  a_a^) 


(B.25) 


Remembering  that  the  conjugate  transpose  distributes  over  tensor  products,  we  see  that 
this  X  is  indeed  antihermitian.  We  compute 


[X,  Hint]  =  2%[<r+a,  O-a))  =  2hg(a+a-aa ^  —  a~a+a^a) 

=  2%[|e)(e|(l  +  a) a)  -  |c/)(sf|ata]  =  2hgaz(a^a  +  ^),  (B.26) 


where  we  have  subtracted  a  constant  offset  of  hg.  Using  this  to  calculate  (B.23)  to  second 
order  in  (^)2  and  factoring  yields  the  dispersive  Hamiltonian 


Hd 


isp 


h 


Ur  +  A 


1  ^  +  ^  +  \hwa(Tz. 


(B.27) 


Remarkably,  we  see  that  this  is  effectively  a  Hamiltonian  for  an  “uncoupled”  cavity  and 
atom,  but  with  the  cavity’s  effective  resonance  frequency  ui'r  =  ay  ±g2/ A  dependent  on  the 
state  of  the  qubit!  This  is  what  causes  the  shift  in  transmission  resonance  which  we  can 
detect  to  infer  the  state  of  the  qubit  without  disturbing  it. 

We  can  also  regroup  the  dispersive  Hamiltonian  terms  as 


Hd 


isp 


hujr 


a'a+- 


h 
+  2 


2g2  ^  g2 

coa  +  -^-a*a  +  ^ 


°z- 


(B.28) 


In  this  form  we  can  interpret  the  effect  of  the  dispersive  interaction  as  creating  a  new 
effective  atomic  transition  frequency  u'a  =  ay  +  ^-a^a  +  ^  (a  photon-induced  Stark  shift 
of  2g2n/A  and  a  vacuum  fluctuation-induced  Lamb  shift  of  g2/ A  [7]). 

The  dispersive  eigenstates  |+,n)  «  |e)|n  —  1)  and  |  —  ,n)  ~  |g)|n)  can  be  found  by 
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expanding  (B.17)  to  first  order  in  g/ A: 

I +,n)  «  ^^\g)\n)  +  |e)|n  —  1) 

\-,n)^\g)\n)-^^\e)\n-l).  (B.29) 

The  n  =  1  state  only  has  one  state  (the  ground  state)  to  decay  into,  and  so  we  can  expect 
by  the  simple  argument  of  how  much  time  |+,  1)  spends  as  a  photon  (remember  vacuum 
Rabi  oscillations)  that  the  rate  of  decay  of  the  atom  from  its  excited  state  is 

r  ~  (£)  k  +  7±  =  +  7±.  (B.30) 

We  thus  see  that  if  n  >  J±,  the  decay  rate  of  the  atom  can  be  enhanced  an  amount  by 
the  presence  of  the  cavity.  This  is  a  special  instance  of  the  Purcell  effect  [71,72],  and  can 
be  used  in  single  photon  sources  for  quantum  optics  experiments  by  encouraging  a  single 
atomic  excitation  to  decay  through  a  large  n  cavity  as  a  photon  [7].  This  effect  can  be 
the  dominant  source  of  relaxation  for  qubits  in  the  vicinity  of  or  above  the  lowest  cavity 
mode  [72],  although  at  lower  bias  points,  dielectric  loss  is  currently  believed  to  be  dominant 
in  Josephson  junction-based  superconducting  qubits  [13,14,73]. 


Appendix  C 

Transmission  Line  Circuits 
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Figure  C.l  Distributed  circuit  model  of  transmission  line.  Figure  adapted  from  [7]. 


In  what  follows,  x  denotes  the  length  coordinate  along  the  line,  V (x,  t )  is  the  potential  on 
the  center  conductor  relative  to  the  ground  planes,  and  I(x,  t )  denotes  the  current.  As  shown 
symbolically  in  Fig.  C.l,  suppose  the  CPW  (or  some  other  one-dimensional  waveguide 
or  transmission  line)  has  a  uniform  capacitance  to  ground  Ci ,  total  self-inductance  Li, 
resistance  Ri,  and  shunt  conductance  to  ground  Gi  per  unit  length.  Following  Pozar  [16], 
Kirchhoff’s  laws  in  infinitesimal  form  then  read 


dV{x,t) 

dx 

dl(x,t ) 
dx 


- Ril(x,t )  -  Li 
- GLV(x,t )  -  Ci 


dl(x ,  t ) 

at 

dV(x,t) 


at 


(C.la) 

(C.lb) 


For  sinusoidal  steady-state  signals  of  frequency  oj,  these  then  become  (when  using  the 
complex  phasor  representation  Ael^elut  for  V  and  I) 

=  -{Rt  +  iuLi)I{x)  (C.2a) 

dx 

- ®  =  ~{Gi  +  iuCi)  V{x).  (C.2b) 
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These  transmission  line  equations  combined  yield 

x )  \  d2I(x )  2 

dx2  “  7  ’  ~~d^  “  7  ^  (C‘3) 

where  7(0;)  =  (ify  +  iuiLi)(Gi  +  iuiCi )  =  a(oj)  +  i(3(uj)  is  the  complex  propagation  con¬ 
stant.  The  traveling  wave  solutions  of  (C.3)  are 

V(x)  =  V0+e~^x  +  V0~e^x  (C.4) 


and  the  analogous  result  for  for  I(x).  From  this  we  see  that  a 
j 3  =  2n/X  =  w/nphase  =  WyTelr/c  contains  wave  parameters, 
yields 

I(x)  =  T(C+e-^  -  V0-e^), 

Zq 


represents  dissipation  while 
Plugging  (C.4)  into  (C.2a) 

(C.5) 


where  Zq  =  y  the  characteristic  impedance.  Let  us  consider  the  termination  of 

the  TL  by  a  load  impedance  Zl  (x  =  0  at  termination).  Using  (C.5), 


„  -V(0)  U0+  +  Vq  ^ 

L~  W)  v+-v0-  °' 

(C.6) 

from  which  we  obtain  the  reflection  coefficient 

r  _  Vo~  _  ZL  ~  Zq 

U0+  Zl  +  Zq 

(C.7) 

Note  that  no  reflection  occurs  when  Zl  =  Zq,  which  is  known  as 
and  (C.5)  then  give 

impedance  matching.  (C.4) 

V{x)  =  U0+  (e_73;  +  Te7a;) 

(C.8a) 

u+ 

I(x)  =  -Te7*) 

Zq 

(C.8b) 

The  impedance  at  a  point  x  =  —l,  which  can  be  interpreted  as  the  impedance  of  a  TL  of 
finite  length  l,  is  then 


V{—1)  _  _  z  ZL  +  Zotanhy/ 

/(— V)  °  Zq  +  Zl  tanh  7 1 


(C.9) 


after  a  few  lines  of  algebraic  manipulation. 
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